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CAPITULO 0 


Introduccion y Revision de 
Electricidad y Magnetismo 


I. Revision de elementos matematicos 


En este apartado se introducen los elementos basicos del calculo vectorial 
que seran utiles en el desarrollo de los contenidos de los apuntes. Tambien 
se hace un repaso breve del problema de las magnitudes y unidades en las 
leyes fisicas. 


0.1 DEFINICIONES BASICAS 


0 . 1.1 _ 

Definiciones de campo 

Definimos un campo escalar ip como: 

ip : D ClZ 3 — >n 

funcion escalar y de punto que asocia a cada punto x en su dominio D un 
escalar ip. 

Definimos un campo vectorial A como: 

A: D ClZ 3 —> H 3 

funcion vectorial y de punto que asocia a cada punto x en su dominio D un 
vector A. 

Las versiones complejas de los campos escalares y vectoriales se definen de 
manera analoga. 

0 . 1.2 _ 

Vectores unitarios y cosenos directores 

a= Vector unitario en la direction de a: 

a = a/\a\ 
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Si los vectores estan definidos en el cuerpo complejo se debe emplear: 

a = a/Va- a* 

Cuando empleemos el sistema de coordenadas cartesiano, podemos expresar 
cualquier vector unitario f en funcion de los cosenos de los angulos (a, /3, 7 ) 
que forma con los ejes x, y, z: 

r = cos ax + cos (3y + cos 7 5 

conocidos como cosenos directores. 


0.1.3 _ 

Lfneas de campo 

Las lineas de campo son lineas tangentes al campo vectorial en cada punto. 
Indican direction y sentido e intensidad, para lo cual se dibujan de forma 
que hay mayor densidad de lineas en las zonas donde la intensidad de campo 
es mayor. 

En dinamica de fluidos, una linea de campo es precisamente una trayectoria 
recorrida por una particula de fluido. 

0.1.4 _ 

Flujo de A a traves de una superficie S 

Dado un campo vectorial A decimos que su flujo a traves de una superficie 
S viene dado por: 



En dinamica de fluidos, donde A representa un campo de velocidades, el 
diferencial de flujo mide la cantidad de fluido que pasa por el paralelogramo 
tangente por unidad de tiempo. 

0.1.5 _ 

Radian y estereorradian 

El radian es la medida de un angulo piano y se define como el angulo piano 
cuyo vertice, en el centro de una circunferencia de radio r, subtiende un 
arco cuya longitud es r. 

La medida del angulo solido es el estereorradian , que corresponde al angulo 
solido cuyo vertice, en el centro de una esfera de radio r, subtiende una 
superficie esferica cuyo area tiene como valor r 2 . 

Dado que la superficie de una esfera es S = 4 - 7 rr 2 , hay 47 r estereorradianes 
en una esfera cerrada. 

El elemento infinitesimal de area ds sobre la superficie de la esfera de radio 
r viene dado por: 

ds = r 2 sin Odddf) (to 2 ) 

por lo que el elemento diferencial de angulo solido sera: 

ds 

dfl = — = sin 9d6d(j) (sr) 



0.2. Sistemas de coordenadas 9 


y vemos que se cumple, efectivamente que: 


dfi = 4n 


0 . 2.1 _ 

Definition 


0.2 SISTEMAS DE 
COORDENADAS 


Dependiendo del problema y de las simetrfas que presente, es conveniente 
elegir el sistema de coordenadas. 

■ Coordenadas rectangulares ( x,y,z ). 

■ Coordenadas cilmdricas ( p,<j>,z ). 

■ Coordenadas esfericas (r, 9, <f>). Notar que el angulo 9 es de revolution 
con respecto al eje z y varta entre 0 y 7r mientras que el <f> se mide 
siempre en el piano xy y varia en un margen de 0 a 27 t. 


0 . 2.2 _ 

Cambios de coordenadas rectangulares-cilindricas 

Las coordenadas se relacionan segun: 

x = p cos <j) 
y = p sin (f> 
z = z 

y los vectores unitarios (que dependen de </>): 


A p 


cos 4> sin 4> 0 


A x 

Afj, 

= 

— sin (f> cos 4> 0 


Ay 

. ^ . 


0 0 1 


A z 

A x 


cos <f> — sin (f> 0 


A P 

Ay 

= 

sin (f> cos 4> 0 


A,/, 

A z 


0 0 1 _ 


A z 


0.2.3 _ 

Cambios de coordenadas cilindricas-esfericas 

p = r sin 9 
z = r cos 9 

y los vectores unitarios (que dependen de 9): 


A r 


sin 9 

0 cos 9 


' A p - 

Ag 

= 

cos 9 

0 — sin 9 


Ap 

Ap 


0 

1 0 


. ^ . 

Ap 


sin 9 

cos 9 0 


A r 

Ag 

= 

0 

0 1 


Ag 

A z 


cos 9 

— sin 9 0 


Ap 
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0.2.4 _ 

Cambios de coordenadas rectangulares-esfericas 

Las coordenadas se relacionan segun: 

x = r sin 9 cos (f> 
y = r sin 9 sin </> 
z = r cos 9 

y los vectores unitarios (que dependen de 9 y <f>): 


A r 


sin 9 cos (p sin 9 sin <f> cos 9 


A x 

Aq 

= 

cos 9 cos </> cos 9 sin <f> — sin 9 


Ay 

A^ 


— sin (j> cos (f> 0 


. . 

A x 


sin 9 cos (f> cos 9 cos </> — sin (f> 


A r 

Ay 

= 

sin 9 sin 4> cos 9 sin (f> cos </> 


A& 

. A, . 


cos 9 — sin 9 0 


A<p 


0.3 CALCULO DIFERENCIAL 
VECTORIAL 


0.3.1 _ 

Gradiente de una funcion escalar de punto ip 

Sea ip(x, y, z) una funcion escalar y de punto. Se define el gradiente de ip en 
coordenadas rectangulares como: 

dip(x, y, z) „ dip(x, y, z) , dip(x,y,z) „ 

Vip(x,y,z) = ---xH--- y-\ --- 2 : 

ox ay dz 


Interpretation geometrica del gradiente 

Se demuestra que si Vip 7 ^ 0, entonces \7ip apunta en la direction en la cual 
ip esta creciendo mas rapidamente. 


Propiedad 

Sea S la superficie que consta de aquellos ( x,y,z ) tales que f(x,y,z ) = 
k= cte. El piano tangente a S en un punto ( x 0 , y 0 , z 0 ) de S esta definido por 
la ecuacion: 


< V^(x G , y Q , z 0 ), (x - x Ql y - y Q , z - z Q ) >= 0 
si Vi/’(x 0 , y 0 , z 0 ) 7 ^ 0, siendo < • > el producto escalar. 

0.3.2 _ 

Laplaciana de una funcion escalar de punto if) 

Sea i/j(x, y, z) una funcion escalar y de punto. Se define la laplaciana de ip 
en coordenadas rectangulares como: 

d 2 ip{x,y,z) d 2 ip(x,y,z) ff 2 ip{x,y,z ) 

dx 2 dy 2 dz 2 


A ip(x,y,z) = X/ 2 ip(x,y,z) 
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0.3.3 _ 

Divergencia de una funcion vectorial de punto A 

Sea A(x, y, z) = A x (x, y , z)x + A y (x, y, z)y + A z (x, y, z)z una funcion vecto¬ 
rial y de punto. Se define la divergencia de A en coordenadas rectangulares 
como: 

V A = 9Ax I dAy I dAz 
dx dy dz 


0.3.4 _ 

Rotacional de una funcion vectorial de punto A 


Sea A = A x x + A v y + A z z una funcion escalar y de punto. Se define el 
rotacional de A en coordenadas rectangulares como: 


V x A = 


dA h 


f d A 

\ dy dz 
que simbolicamente puede expresarse como: 


(dA x 

I dz 


dA z 

dx 


Vxi = 


x y z 
d/dx d/dy d/dz 


A x 


A,, 


A, 


(dAy 

V dx 


dA x 

dy 


0.3.5 _ 

Laplaciana de una funcion vectorial de punto A 

Sea A(x, y, z) = A x {x , y, z)x + A y (x, y, z)y + A z (x, y , z)z una funcion vecto¬ 
rial y de punto. Se define la laplaciana de A en coordenadas rectangulares 
como: 

/SA = A A x x + A A y y + A A z z 


0.3.6 _ 

Teorema de la divergencia de Gauss 


Sea V una region y S la superficie cerrada orientada que acota V. Sea A un 
campo vectorial definido en V. Entonces: 


V • Adv = (ft) A ■ ds 


El significado fisico de la divergencia es que en un punto P, V • A es la razon 
del flujo neto que sale por P por unidad de volumen. 


0.3.7 _ 

Teorema de Stokes 


Sea S una superficie orientada clefinida por una funcion C 2 , z = f(x,y), 
(x, y) € D y sea A un campo vectorial C 1 en S. Entonces, si L denota una 
curva frontera orientada de S, tenemos: 


V x A ■ ds = (b A ■ dl 
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El significado fi'sico del rotacional en un punto P, V x A es la circulation 
de A por unidad de area en una superficie perpendicular a ds. 


0.3.8 _ 

Operadores en coordenadas cilindricas 


Vip 


dip , 1 dip - 

dp p+ ~ p ^ 



_ 1 d(pAp) 1 9A 0 dA z 
P dp P d<p dz 


1 9x1, 

dA^ 

P + 

\ dA p 

dA z - 

0 + 

'1 9(pA 0 ) 

1 dA p ' 

p dcp 

dz 

dz 

dp 

p dp 

p dcp 


Aip = \7 2 ip 


1 d ' dip 
p dp P dp 


1 d 2 ip d 2 ip 

pi Qfjyl Q z 2 


AA = 


AA P 


2 dAt 
p 2 dcp 



A A^ + 


2 dA p 
p 2 dcp 



<p + A A z z 


0.3.9 _ 

Operadores en coordenadas esfericas 


Vip 


dip , l dip* 1 dip* 
dr r d6 r sin 0 dcp 


V -A = 


1 d{r 2 A r ) 
r 2 dr 


1 d 
r sin 0 dO 


(sin OAg) 


1 dA ^ 

r sin 0 dcp 


Vxl = 


1 

r sin 6 


d_ 

dO 


(A^ sin 0) 


9A e l „ 1 
- rH— 

dcp J r 


1 9A r 
sin0 dcp 


d{rA 0 ) 
dr 


n , 1 [ 5 Me) 

^- 

r or 


9A r 

“90~ 


= VV 


1 9 

' 2 dip~ 

1 9 

I 

’ . ,9V>' 

r 2 9r 

dr 

r 2 sin 0 90 

smU do_ 


1 d' 2 ip 
r 2 sin 2 0 9<^ 2 


AA = 


Ad r 


AA e - 



" ^ ( Ar + + CSC0 ^ + °w) 

^ (csc 2 0Ae - + 2cotOcscO°^j 

(esc 2 OAg, — 2csc0-—2— 2cot0csc0 —-/■ 
\ dcp dcp 


f+ 

0+ 

1 * 
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0.3.10 _ 

Identidades vectoriales 

Suma y multiplication 

A- A= |,4| 2 

A - A* = |A | 2 

A + B =B+A 

A ■ B = B ■ A 

Ax B = -B x A 

(.A + B)-C = A-C + B-C 

(A + B) x C = A x C + B x C 

A-BxC = B-CxA = C-AxB 

Ax B xC = (A- C)B - (A-B)x B 

(Ax B) ■ (C x D) = A - B x (C x D) = A(B ■ DC - B ■ CD) 

(Ax B) x (C x D) = (Ax B ■ D)C -(Ax B- C)D 


Diferen ciacion 


V • (V x A) = 0 

V x V'l' = 0 
V($ + 't) = V$ + 

V-(A + B) = V-A + V-B 
Vx(A + B)=VxI+VxB 
V-(^l) = A ■ V\t + • A 

V x (1/.4) = V4> x A + x A 

V(A-B) = (A-\7)B+(B -\7)A + Ax (VxB)+Bx(Vxl) 

V-(AxB) = B- VxA — A-VxB 

Vx(lx5) = 4V • B - B\7 ■ A + (B ■ V)A- (A-V)B 

V x V x A = V(V • A) - AA 


Para expresar una magnitud en un cierto espacio de referenda, se toma 
como patron la unidad U de forma que una magnitud X puede expresarse 
de la forma: 

X = Ux = U'x' 


0.4 UNIDADES, SISTEMAS 
DE UNIDADES Y 
ECUACIONES DE 
DIMENSIONES 


siendo x la medida de esta magnitud, que variara segun la unidad elegida. 

Para comparar cantidades de la misma magnitud es necesario efectuar sus 
medidas con la misma unidad. 


Las leyes fisicas relacionan entre si medidas de magnitudes de distinta na- 
turaleza. Las leyes fisicas, en general, se expresan mediante relaciones de la 
forma: 

x = Kx^ 1 ■ x^... • x°C ( 1 ) 

donde K es una constante y ai, a 2 , ..., a„ son numeros racionales positivos 
o negativos, de manera que: 


Las magnitudes relacionadas mediante una ley se dice que son cohe 
rentes entre si cuando la constante K es la unidad. 
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■ Las magnitudes fisicas se dicen fundamentals cuando son indepen- 
dientes entre si, es decir, no hay ninguna ley ffsica que las relacione. 
Las unidades que tomemos para la medida de estas, las denominare- 
mos unidades fundamentales. 

■ Toda magnitud ffsica relacionada con las fundamentales mediante una 
ley ffsica se denomina magnitud derivada. Sus unidades se denomi- 
naran unidades derivadas y se obtendran mediante una relation del 
tipo ^ entre las unidades de las magnitudes fundamentales. 

Como magnitudes fundamentales comunes a todos los campos de la ffsica, 
se han tornado la longitud y el tiempo; ellas son suficientes para desarrollar 
la cinematica; en los demas campos de la mecanica se hace necesario utilizar 
una tercera magnitud fundamental, que puede ser la masa o la fuerza. En los 
restantes campos de la ffsica es necesaria una cuarta magnitud fundamen¬ 
tal, tomandose en termodinamica la temperatura; en optica la intensidad 
luminosa y en electricidad el amperio. 

Un sistema de unidades queda definido cuando se han elegido sus magnitu¬ 
des fundamentales con sus respectivas unidades de medida. A partir de estas 
y mediante las correspondientes leyes ffsicas se deduciran sus magnitudes y 
unidades derivadas. 

M.K.S.A.= metro, kilogramo, segundo, 
amperio 

El sistema C.G.S. es una variante del an¬ 
terior que utiliza el centi'metro, el gramo 
y el segundo. 

■ L para la longitud 

■ T para el tiempo 

■ M para la masa 

Si en la expresion de una ley ffsica sustituimos cada magnitud fundamental 
por su sfmbolo, obtendremos la llamada ecuacion de dimensiones de la 
magnitud derivada de que se trate. 

La ecuacion de dimensiones nos permite establecer la homogeneidad de las 
formulas ffsicas, diciendose que una expresion es homogenea cuando ambos 
miembros tienen la misma ecuacion de dimensiones. Ademas facilita el esta- 
blecimiento de la coherencia entre unidades ffsicas, diciendose que una serie 
de unidades son coherentes entre sf, cuando la ecuacion de dimensiones de la 
ley ffsica que las relaciona no contiene factor numerico multiplicative salvo 
la unidad. 

un sistema de unidades U , U\, C/ 2 , ■ U n para el cual es X = Ux , X\ = 
U 1 X 1 , X 2 = U 2 X 2 , X n = U n x n se dice que es acorde con una ley ffsica 
como la ([l]) cuando verifica: 

x = Kx“ 1 -x 2 2 ... •<" (2) 

Un nuevo sistema de unidades U', U[, U 2 , ..., U' n para el cual es X = U'x', 
X\ = U[x' ll X 2 = U 2 X 2 , ..., X n = U' n x' n se dice que es no acorde con una 
ley ffsica cuando cambia la cosntante de la misma, es decir, verifica: 


En lo que sigue utilizaremos el sistema M.K.S.A que ha pasado a formal' 
parte del Sistema Internacional (S.I.), adoptado oficialmente en Espaha 
desde 1967. 

Cada magnitud fundamental se representa por un sfmbolo que en el Sistema 
Internacional y para las magnitudes mecanicas es: 


x' = K'x'l 1 • x ,a 2 2 ... ■ x ,a n 


(3) 


0.4. Unidades, sistemas de unidades y ecuaciones de dimensiones 


El cociente entre ®)y ©: 


x 

x' 


K 

~K' 





Q n 


en funcion de las unidades de los respectivos sistemas se expresa: 

u = K(u 1 y i (u^y* (Un\ a " 
v K'\u'J [uj -\uj 


expresion que nos permite cambiar de sistema de unidades en la expresion 
de una ley fisica, al darnos el coeficiente k' de la expresion correspondiente 
al nuevo sistema de unidades: 






(Un 

W, 


a n 


(4) 


Una misma ley fisica, tal como la ([I]), en un cierto sistema de unidades, se 
expresa mediante una relacion como la ([2]) , mientras que en otro sistema de 
unidades se expresa mediante la <§, estando dado el coeficiente K' por la 
ecuacion 0- 
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0.5 ELECTROSTATICA 


Electricidad deriva de qXenrpou que sig- 
nifica ambar en griego 


La unidad de carga en el sistema MKSA 
(SI) es el Coulombio (C). La carga de un 
electron es de —1.602 x 10 —19 C 


II. Revision de Electricidad y Magnetismo 


En este apartado se revisan los conceptos mas importantes de la electricidad 
y magnetismo clasicos hasta llegar al planteamiento de las ecuaciones de 
Maxwell. 


0.5.1 _ 

La carga 


Introducimos la funcion p(r) conocida como densidad de carga 


p{r) = 


dq 

dv 


La carga esta cuantizada, de forma que la carga total de un cuerpo es 
un multiplo entero de la carga de un electron. Las distribuciones de carga 
que emplearemos son continuas. Por ello, consideraremos elementos de vo- 
lumen suficientemente grandes como para englobar un numero elevado de 
parti'culas cargadas, de forma que deje de manifestarse su caracter discreto 
y veamos efectos de promedio. 

Se entiende por carga puntual situada en r*i la definida como: 

p(r) = q 5(r- fq) 

donde £(•) es una funcion delta de Dirac. Una carga puntual tiene sentido 
al considerar distribuciones de carga a distancias muy grandes respecto a 
su tamano. 


0.5.2 _ 

Ley de Coulomb y principio de superposicion 


Hay dos leyes fundamentals en la electrostatica: la Ley de Coulomb y el 
principio de superposicion. Supongamos dos cargas puntuales de valor q y 
q' separadas una clistancia r. La Ley de Coulomb establece que la fuerza 
entre ellas responde a la expresion: 


donde C es una constante que depende del medio. En un sistema racionali- 
zado la ecuacion anterior resulta: 


f = 


qq' , 

- T 

Airer 2 


siendo e una constante conocida como permitividad o constante dielectrica 
del medio. 

Cuando las cargas se encuentran en el vacio, clecimos que e = e D . El valor 
de esta constante en MKSA es (notar que tiene dimensiones): 


e o = 9 (Faradios/m .) = 8.854 x 10 12 (F/m) 

3o7r 
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Para medios lineales, homogeneos e isotropos distintos al vacfo, la permiti- 
vidad se expresa como: 

e = e r e Q 


siendo e r la permitividad relativa del medio, que cumple siempre e r > 1. Mas 
adelante comentaremos brevemente como se introducen los medios materia- 
les en la description del fenomeno electromagnetico. En cuanto al principio 
de superposition, establece que la fuerza que experimenta una carga dada 
debida a una distribution de cargas es la suma vectorial de las fuerzas que 
producirfan cada una de las cargas de la distribution si estuvieran aisladas 
del resto. 


Un medio es homogeneo cuando sus ca- 
racterfsticas no dependen de las coorde- 
nadas, e isotropo, cuando las propieda- 
des del mismo vari'an igual en todas las 
direcciones en tor no a un pun to dado. 


0.5.3 _ 

Definition de £ 


Aunque las dos leyes anteriores son suficientes para resolver cualquier pro- 
blema de la electrostatica, es convieniente por razones practicas y concep- 
tuales introducir dos ideas secundarias que seran de gran utilidad: el campo 
electrico y el potencial electrostatico. 

Se define el campo electrico £ creado por q' en la position de la carga q 
como la fuerza que experimenta la unidad de carga: 


de donde se sigue que: 



£ = 


q’ - 

- r 

47rer 2 


(5) 


Observamos que se trata de un campo de fuerzas centrales. Sus unidades 
mas habituales en MKSA son N/C 6 Voltio/m = V/m. 

El principio de superposition se expresa como: 


N 




i—1 


47 t erf 


0.5.4 _ 

Definition de V 


Resulta de interes introducir un nuevo vector mas relacionado con las cargas 
que £. Definimos: 

V = e£ (6) 

asumiendo una relation lineal entre £ y V. Este nuevo campo se denomina 
vector desplazamiento electrico o vector de densidad electrica de flujo. Para 
cargas puntuales, este vector tiene la propiedad de que es radial e indepen- 
diente del material en cuyo seno se encuentran las cargas. Ademas su flujo 
a traves de una superficie cerrada que englobe a la carga es igual a la carga 
total: 

(jj) V ■ ds = T,q t 

y esto es independiente de la superficie en cuestion. 


Cuando se trabaja con medios materiales reales (no en el vacio), la presen- 
cia de un campo electrico en dichos medios puede provocar que las cargas 
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Recordamos que una distribution con un 
momento dipolar V e crea un campo ex- 
terno equivalents al que producirfa una 
distribution volumetrica de carga p p = 
—V • V e y una distribution superficial 

ps = Pen 


de las moleculas (neutra en un principio) se desplacen de forma que dichas 
moleculas se polarizan. Tambien puede ocurrir que la molecula inicialmen- 
te se encuentre polarizada (o lo que es lo mismo, que tenga un momento 
dipolar permanente), por lo que el campo hara que se oriente de determi- 
nada forma. Grosso modo, el campo final sera la suma del campo inicial 
mas las contribuciones de los distintos dipolos y la distribution de campo 
en el interior del material sera diffcil de determinar y presentara grandes 
variaciones a nivel microscopico. De acuerdo con la evidencia experimental 
se verifica: 

V = ej + V e (7) 

El vector V e es el vector adicional de polarization electrica. 

En medios lineales, homogeneos e isotropos 0 se convierte en: 

V = e 0 £ + V e = + CoXe£ = £o(l + Xe)£ = 

donde Xe es una constante positiva que depende del material y se denomina 
susceptibilidad electrica 


0.5.5 _ 

Propiedades integrales y diferenciales de £ y V 


Segun el teorema de Helmholtz, un cam¬ 
po vectorial queda um'vocamente deter- 
minado si se conocen su divergencia y su 
rotational en todos los puntos del espa- 
cio 


Estudiamos las propiedades integrales y diferenciales de £ a partir de la 
expresion ([5| utilizando resultados conocidos del calculo vectorial. 

Se demuestra que se verifica la ecuacion: 


V ■ ds = 


pdv 


donde S es la superficie cerrada que encierra el volumen V. 

Esta expresion se conoce como Teorema de Gauss. Esto significa que el flujo 
neto a traves de una superficie cerrada que no encierra ninguna carga es 
cero. 

La version diferencial de este teorema es: 


V V= p 


Por otra parte, es facil demostrar que 

Vxf = 0 


y aplicando el teorema de Stokes a esta integral, resulta: 


Un campo que verifica 0 se denomina 
conservative. 


E-dl = 0 


( 8 ) 


0.5.6 _ 

Potential electrostatico 

Esta propiedad permite una simplification. Sabemos que un campo que 
cumpla esta propiedad se puede calcular como: 



(9) 
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siendo $ una funcion escalar y de punto denominada potencial electrostdtico 
que tiene una propiedad muy util: la diferencia entre dos puntos es igual al 
trabajo que hay que realizar para trasladar una carga unidad entre ellos. 

En efecto, el trabajo que hay que realizar para transportar una carga del 
punto 1 al punto 2 contra el campo £ sera: 


W- l2 = ~ f £-dl = <f > 2 -$i 


Asi, la diferencia de potencial entre dos puntos es igual al trabajo que hay 
que realizar para trasladar una carga unidad entre los dos puntos. La ecua¬ 
cion <f>(r) = tie. define una superficie equipotencial, siendo el campo per¬ 
pendicular a estas superficies en todo punto. 

Si ahora tomamos la divergencia de ([9]), resulta: 



se obtiene la conocida como Ecuacion de Poisson, que para puntos del 
espacio donde p = 0 se convierte en: 


Ad> = 0 


conocida como Ecuacion de Laplace. 

Teniendo en cuenta que el potencial creado por una carga puntual es: 


9 47rer 

para un elemento de carga pdv' sera d<f> = pdv'/\r— r*j por lo que el potencial 
total sera, por superposition: 






dv' 


que corresponde a la version integral de la ecuacion de Poisson. Esta ecua¬ 
cion plantea un problema de suma en contraste con los problemas de con- 
diciones de contorno, bastante mas complejos, en los que p(r) no se conoce 
de antemano. 


0.6 MAGNETOSTATICA 


0 . 6.1 


Corriente y conceptos relacionados 

Supongamos una distribution de carga que se mueve con velocidad v cons- 
tante. Definimos un vector densidad superficial de corriente electrica J en 
cada punto como: 

J = pv (10) 


La palabra magnetismo deriva de Mag¬ 
nesia, region de Grecia donde se encon- 
traron las primeras piedras iman con es¬ 
tas propiedades 


esto es, como la cantidad de carga que por unidad de tiempo atraviesa una 
unidad de area (A/m 2 en MKSA). 

Supongamos ahora una densidad de corriente dependiente del tiempo J(r , t) 
y una superficie infinitesimal ds situada en el seno de la corriente y sea 
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dl la carga que atraviesa la superficie por uiiidad de tiempo. Claramente 
dl = J ■ ds por lo que se verificara: 


I se denomina corriente a traves de la superficie S. 

Supongamos que la superficie S es cerrada y que admitimos el principio de 
que la carga no se crea ni se destruye. Entonces, la carga total en el volumen 
V sera f v pdv y su tasa de variation con el tiempo d/dt f v pdv. Por tanto, 
el principio de conservation de la carga se expresa como: 


J ■ ds+ — / pdv = 0 
s ot Jv 


Si aplicamos el teorema de la divergencia al primer sumando podemos susti- 
tuirlo por f v V • J ds y dado que la expresion anterior debe verificarse para 
cualquier volumen V, se sigue que: 




( 11 ) 


expresion matematica de la Ley de conservacion de la carga o ecuacion de 
continuidad de la carga. 

Si la densidad de carga p es una funcion solo de r y no varfa con el tiempo, 
la ecuacion anterior se convierte en: 


V-J = 0 


que implica que las lfneas de flujo de corriente se cierran sobre si mismas, 
esto es, forman lazos cerrados. 

Cuando la corriente se produce dentro de un conductor y dentro del conocido 
como margen lineal , se verifica experimentalmente la conocida Ley de Ohm 
(generalizada): 

J = a £ 

siendo cr la conductividad del medio en MKSA). El valor <j per- 

mite clasificar a su vez a los medios como conductores, semiconductores y 
aislantes, si bien esta division es relativa, como veremos mas adelante. 

A partir de la anterior se deriva la Ley de Joule: 

^ = J-£ = a\£\ 2 

dv 


0 . 6.2 _ 

Fuerzas entre conductores 

Supongamos dos lazos por los que discurren las corrientes I\ y I 2 respec- 
tivamente. Los experimentos de Ampere mostraron que la fuerza total F\ 
sobre el circuito 1 debida a su interaction con el circuito 2 es proporcional a 
las dos corrientes I\ y I 2 y a una integral que depende solo de la geometrfa 
(ver figura [l]) : 

A = CI J2 ( f ) 

Jc 1 

La constante de proporcionalidad C depende del medio. En un sistema 
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Figura 1: Ley de Ampere 


racionalizado, la ecuacion anterior se escribe: 


* = 


dh x ( dl 2 x fi 2 ) 


Cl Jc 2 


12 


siendo /i una constante conocida como permeabilidad magnetica del medio, 
cuyas unidades son ( Henrio/m)=H/m en MKSA. La permeabilidad del 
vacio en este sistema es 


/io = 47 t x 10 7 ( H/m ) 

Para cualquier otro medio lineal, homogeneo e isotropo, podemos escribir: 


p P’rl^o 


0.6.3 _ 

Definicion del vector de intensidad de campo magnetico Tt 


Partimos de una situation en la que tenemos corrientes electricas estaciona- 
rias, esto es dj/dt = 0. Tal y como hicimos en electrostatica, es conveniente 
definir un nuevo campo, imaginando que una de las corrientes produce un 
campo que actiia sobre la otra. 


Definimos el vector intensidad de campo magnetico Tt producido por / 2 
como: 



dl 2 X fi2 


r 


3 

12 


Entonces, postulamos que este campo Tt ejerce una fuerza T\ sobre la 
corriente I\ dada por: 

T\ = pl\ (j) dli xTt 

siendo las unidades de Tt A/m en MKSA. 

Es conveniente extender la definicion de TL asumiendo que cada elemento 
dl de corriente produce un campo dTt dado por: 


dTL = 


/ dl x r 
4-7T r * * 3 


siendo r el vector que une el elemento y el punto en el que se calcula el 
campo. Esta expresion recibe el nombre de Ley de Biot y Savart. Notar que 
implicitamente estamos utilizando el principio de superposicion de fuerzas 
magnetic as. 
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0.6.4 _ 

Definition del vector de induccion magnetica B 


De la misma forma que hicimos con el campo electrico en electrostatica, es 
necesario introducir el efecto de lo microscopico en nuestro planteamiento 
macroscopico. A efectos de calcular el campo magnetico, los medios materia- 
les se caracterizan mediante momentos dipolares magneticos. Introducimos 
un vector vector de induccion magnetica B definido como: 

B = nH (12) 

que puede expresarse como 

B = Hoii + Vm (13) 


siendo V m el vector adicional de magnetization. En medios lineales, ho- 
mogeneos e isotropos (131 se convierte en: 


B = /I 0 H +V m = Ho'H + HoXm'H 


donde Xm es una constante positiva o negativa que depende del material y 
se denomina susceptibilidad magnetica 


0.6.5 _ 

Propiedades diferenciales e integrales de B y H 

Operando a partir de las ecuaciones anteriores, se puede demostrar que: 

V xH = J (14) 

V B = 0 (15) 

esto es, el campo B es rotational puro. Empleando los teoremas de la diver- 
Notar que e! campo B no es conservative: gencia y Stokes llegamos a las versiones integrales: 


H ■ dl = I 


ic 


conocida como Ley de Ampere, junto con 

B ■ ds = 0 


0 . 6.6 _ 

El Potential vector magnetico 

Dado que V x B es distinto de cero, no puede existir un potencial real 
para calcular el campo mediante una diferenciacion. No obstante, podemos 
introducir un potencial generalizado de tipo vectorial que si cumple esta 
propiedad. Ademas el resultado sera valido tambien para el interior de las 
corrientes. 

Introducimos un vector A al que llamaremos potencial vector de forma que: 

B = Vxd 
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por lo que se cumplira automaticamente que: 

V B = 0 


Entonces, la ecuacion (14) resulta: 


V x (V x A) = pj 


que corresponde a un conjunto de ecuaciones diferenciales acopladas y que 
se puede expresar como: 


V(V ■ A) - AA = nJ 


Por el teorema de Helmholtz, A queda determinado solo cuando definamos 
V • A. Como tenemos libertad para hacerlo como queramos, a la vista de la 
ecuacion anterior, interesa imponer V • A = 0 de forma que: 



Dado que esta ecuacion tiene la misma forma que la ecuacion de Poisson, 
su solution sera analoga, de forma que podemos decir que: 


A?) = 


4-7T 




dv' 


'V' 


0.7 ENERGIA EN CAMPOS 
ESTATICOS 


0.7.1 _ 

Energfa electrostatica 


Suponemos una distribution arbitraria de cargas. Se clemuestra que la energia 
potential del sistema conocida como energia electrostatica viene dada por: 


W E 



P(r)p(r') 

| f- T* | 


dvdv' 


Una variante de esta expresion es: 

We = lj p(r)$(r)dv 

Aunque esta integral se extiende solo a los puntos donde existe densidad de 
carga, puede ser extendida a todo el volumen de espacio Vt de la siguiente 
forma: 

W E =\\ £-Vdv 

z Jvt 


Esta formula expresa la energia linicamente en terminos del campo elec- 
trostatico, concepto basico Maxwelliano (la energia reside en el campo no 
en las fuentes). 
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0.7.2 _ 

Energfa magnetostatica 

Se demuestra que la energia necesaria para formar una distribution de 
corriente viene dada por: 


Wm = ^ J A- Jdv 

Aunque la energia parece residir en las corrientes, esta expresion puede 
modificarse para que la energia se exprese unicamente en funcion del campo: 


W M =■= \ B- Hdv 
1 Jvt 


0.8 ELECTRODINAMICA 


0 . 8.1 _ 

El problema de la variacion temporal 

La electrostatica trata distribuciones de cargas estaticas, esto es, inmoviles 
en promedio. Por tanto, los campos electricos que se producen son indepen- 
dientes del tiempo. La magnetostatica estudia los campos creados por cargas 
en movimiento, con la particularidad de que estas corrientes son estaciona- 
rias, esto es, independientes del tiempo. Los campos electricos y magneticos 
son, por tanto, independientes del tiempo. Ahora vamos a calcular el efecto 
del movimiento de las cargas de manera arbitraria, por los que los campos 
seran variables con el tiempo. 

En espacio libre, en presencia unicamente de cargas y distribuciones de 
corriente pero no de materia, hay cuatro ecuaciones basicas para los cam¬ 
pos estaticos, dos para definir la divergencia y el rotacional de £ y otras 
dos para B. Para el caso electrodinamico, las ecuaciones que definen los 
rotacionales deben completarse con terminos que dependen de la variacion 
temporal. Uno de ellos fue encontrado experimentalmente por Faraday. El 
otro, fue postulado por Maxwell. Este conjunto de 4 ecuaciones completas 
se denomina ecuaciones de Maxwell para espacio libre. 


0 . 8.2 


La Ley de Faraday 

Faraday encontro, empleando corrientes en lazos, que la integral curvilfnea 
del campo electrico alrededor de un lazo es proporcional a la variacion tem¬ 
poral del flujo magnetico a traves de ese lazo: 


£■ dl = - [ ^-ds 

Js dt 


(16) 


La evidencia empirica indica que esta relation entre un campo magnetico 
variable y el campo electrico, es una propiedad de los campos en el espacio: 
los cables y las corrientes que fluyen en ellos sirve unicamente para que se 


manifieste (16 ). 
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Empleando el teorema de Stokes podemos escribir: 


V x £ = - 


dB 

~dt 


que es la generalization de la ley V x £ = 0 de la electrostatica. 

Si generalizamos lo anterior a lazos de forma variable se obtiene la conocida 
Ley de Lenz: 

, 9 $ n 7\ 

9 ~dt ( 17 ) 

siendo <fi el flujo a traves de la superficie que define el lazo y fem la fuerza 
electromotriz generada. 


0.8.3 _ 

La corriente de desplazamiento y las ecuaciones de Maxwell 

Discutimos ahora una generalizacion de las ecuaciones anteriores debida 
a James C. Maxwell. Esta generalizacion es un postulado que se entiende 
mejor si miramos al contexto de la epoca en que se introdujo. 

Consideremos la ecuacion de la magnetostatica Vx6 = pij. Si tomamos 
la divergencia de esta ecuacion obtenemos: 

V -J = 0 

Sabemos que el hecho de que la divergencia de un vector sea cero implica 
que las liiieas de campo se cierran sobre sf mismas. Por lo tanto, la ecuacion 
anterior implica que la corriente fluye en lazos cerrados. Esto es lo que 
ocurre con campos estaticos. Sin embargo ahora, nos encontramos ante una 
situation diferente. 

Si p. ej. consideramos un condensador que se carga conectando las dos pla- 
cas a una bateria, podemos decir que existe corriente en el cable pero no 
habra corriente entre las placas del condensador. La conception de este 
proceso por parte de Maxwell es diferente. Para entenderla acudimos a los 
descubrimientos de Faraday sobre el fenomeno de la polarization o lo que 
es lo mismo, a la separation de cargas en un medio material debido a la 
presencia de un campo electrico. Las cargas deben moverse para separarse 
y, dado que las cargas en movimiento constituyen una corriente, el proceso 
de polarization implica la existencia de una corriente que puede llamarse 
corriente de polarizacion. En aquel momento se supom'a que las acciones 
electromagneticas no se transmitian a traves del vacfo sino a traves de un 
eter que no era distinto en nada a la materia dielectrica. Por lo tanto, no 
era raro que Maxwell extendiera el concepto de corriente de polarizacion al 
eter. 

Maxwell considero que, en el proceso de carga de un condensador, los cables 
del condensador y el eter entre dichas placas formaban un circuito material 
continuo. Para ello, modified las ecuaciones anteriores como se indica a 
continuation. 

A partir de V • V = p se sigue: 
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Si introducimos esto en la ecuacion de continuidad de carga (|Tl| obtenemos: 

J+f =0 (18) 


El segundo sumando dT>/dt tiene las dimensiones de una corriente y fue 
denominado por Maxwell corriente de desplazamiento. Si consideramos que 
la corriente total consiste en una densidad de corriente mas una corriente de 
desplazamiento, vemos que la corriente total siempre fluye en lazos cerrados. 
Sin embargo, debemos recordar que la corriente de desplazamiento no es una 
corriente en el sentido ordinario, ya que no esta asociada a un flujo de carga. 


A partir del resultado (181, Maxwell generalizo la Ley de Ampere diferencial 
V x 7 ~t = J asumiendo que para los campos variables con el tiempo, la 
corriente en esta formula deberia ser la corriente total incluyendo la corriente 
de desplazamiento, de forma que esta ecuacion resulta: 


V x H = J + 


dV 

~dt 


( 19 ) 


Esta ecuacion es un postulado que ha quedado contrastado ampliamente 
por la experiencia. Entonces, las ecuaciones del modelo de Maxwell son: 


V x H = 3 


dV 

ht 


_ f - dB 
Vx£ = -3lT 

y si imponemos como postulado la ley de conservation de carga: 

X7J+^ = 0 

at 

se puede demostrar facilmente que se cumple: 


V V= p 
X-B = 0 


Ademas, en el sistema de ecuaciones planteado, hemos de considerar el tipo 
de medio material (e, p,, a): 

J = a£ 

V = e£ 

B = pH 

conocidas como relaciones constitutivas del medio material 

Las ecuaciones de Maxwell permiten calcular los campos £, V, B y H a partir 
de fuentes arbitrarias. Dado que estos campos son importantes debido a su 
action sobre las cargas, los fundamentos de la teoria electromagnetica se 
completan prescribiendo como es dicha accion. Para ello se utiliza la Ley de 
la densidad de fuerza de Lorentz : clada una distribution de carga arbitraria 
caracterizada por p(f) que se mueve a una velocidad v{r) con respecto a un 
sistema de referencia inertial en el cual hay campos £ y B, la densidad de 
fuerza (fuerza/unidad volumen) que actua sobre la distribution / es: 


f = p(£ + vxB) 







0.8. Electrodinamica 27 


Observamos como el primer sumando es una extension del campo elec- 
trostatico al caso electrodinamico. El segundo sumando es la esencia del 
postulado. Generaliza los resultados obtenidos en magnetostatica sobre la 
fuerza entre dos lazos con corrientes estacionarias. Por ultimo, la energia al- 
macenada (electrica y magnetica) por el campo electromagnetico vendra da- 
da por: 

W E = \ [ £-Vdv 
2 Jvt 

W M = x / B ■ 7 idv 

z Jvt 

Como veremos a lo largo del curso, los conceptos clave del planteamiento, 
ampliamente validado por la experiencia, son: 

■ la introduccion cle V y U. 

m la simetrizacion del sistema (J + dT>/dt no tiene manantiales) 

■ la inseparabilidad £ <-> B: observamos como estan acoplados en las 
ecuaciones y que, aunque se verifique J = 0 existiran fuentes de cam¬ 
po magnetico si dV/dt ^ 0, por lo que si existen campos electricos 
variables, estos se convierten en fuentes de campo magnetico variable 
y lo mismo ocurre al contrario: aunque no existan distribuciones de 
cargas, si existen campos magneticos variables dB/dt ^ 0, estos se 
convierten en fuentes de campo electrico. 
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0.9 EJERCICIOS 


0.2 


Sean los vectores 


a = 3x + 2y — z y b = 


—y + 35. Calcule: 


A . a ■ b 

B. a x b 

C. Angulo que forman a y b 

D. Vector unitario en la perpendicular del piano que forman a y b 

Sol : A) -5; B) 5x - 9y - 35; C) 115°; D) -9 y- 35) 


0.3 


Sean los vectores 


a = 3 jx + 2y - jz y b = 


—y + (1 + 2j)z. Calcule: 


A. \a\ y |6| 

B. |o • 6*| 

C. \ Re[a x 6*] 

Sol : A) \a\ = Vli ; |6| = \/6; B) \a -b*\ = 3; C) l/2i?e[a x b*} = x + 3 y 


0. 4 


Indicar las coordenadas esfericas que definen: 


A. La direccion 5 

B. El piano bisector de los ejes xy 

C. La direccion r= — x + y + 5 

Sol: A) 6 = 0; B) & = tt/4; C) 0 = 135° 6 = 54.73° 


0. 5 | Expresar la curva p = 


Sol : x 2 + y 2 = a^/x 2 + y 2 + 


a(l + cos </>) en coordenadas cartesianas. 
ax 


0 . 6 


Expresar en coordenadas cilindricas el vector r = x — y + 25 en la 


direccion definida por </> = 3ir/2 
Sol : r = p + (j>+2z 


0. 7 


Calcule el gradiente de los campos escalares: 


A. xj) = \/x 2 + y 2 + z 2 

B . tp = r 

C- V>=y 

D . ip = 3x 2 + yz 


Para el cuarto caso, calcule un vector unitario que apunte a la direccion de 
maxima variacion de ip en el punto P Q (1, — 1,0). 


Sol: A) 

(6x 


yj x 2 +y 2 +z 2 

y)x + xy; v = ~^{bx 


7 v „ y+ 7 £, 

y/x 2 +y 2 +z 2 y/x 2 +y 2 +z 2 

V) 


5 = r; B) f; C) -4f; D) 







0.9. Ejercicios 


0. 8 Sea A = 3x 2 yx + (x 3 + y 3 )y. Verificar que A es irrotacional. 
Sol: Vxl= 0 


0. 9 Calcular A(£) para r ^ 0. 

Sol : A>J = 1/r 2 Jy[r 2 §^] con ip = 1 /r => A(£) = 0 


0. 10 I Sea F = k^r un campo vectorial central, donde k y Q son cons- 


tantes. Calcular f v V • Fdv siendo V el volumen definido por una esfera de 
radio R centrada en el origen de coordenadas. 

Sol : kQ4n 


0. 11 Sea U{r, 6.<p) = sin 0. Calcular cuanto vale (fi s Uds. 


0. 12 Demostrar que cualquier solucion de la ecuacion V x V xA—k 2 A = 0 


satisface automaticamente la ecuacion vectorial de Helmholtz A A+k 2 A = 0 
y la condicion solenoidal V • A = 0. 

Sol : Aplicando V- a la ecuacion y simplificando => V • A = 0; Desarrollando 
V x V y aplicando lo anterior => A A + k 2 A = 0 
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CAPITULO 1 


El modelo electromagnetico 


1.1 EL CAMPO 
ELECTROMAGNETICO 

1 1 1 _ SOBRE EL CUERPO REAL 

Ecuaciones locales 


Postulamos que las magnitudes que intervienen en el electromagnetismo 
estan ligadas por las ecuaciones locales: 


dV 

y x7f = J + - 

(1.1) 


(1.2) 

= 0 

(1.3) 

V-P = p 

(1.4) 

o 

II 

t ^q 

o 

(1.5) 

J = a£ 

(1.6) 

V = e£ 

(1.7) 

B = pH 

(1.8) 


Las dos primeras son igualdades vectoriales y se conocen habitualmente con 
la designation de primera y segunda ecuacion de Maxwell , respectivamente. 
La tercera, escalar, traduce el postulado de la conservation de la carga 
electrica y, unida a la cuarta y la quinta, tambien escalares, forma el bloque 
de las ecuaciones de la divergencia. Las tres ultimas, vectoriales, introducen 
las constantes materiales cr, e y /i que caracterizan electromagneticamente 
el medio soporte del campo (relaciones constitutivas del medio material). 


Ha de indicarse que, el sistema de ecuaciones asi establecido es compatible 
a pesar de no ser independiente. Vemos que efectivamente no lo es ya que 
(1.41 es consecuencia de (1.11 y (1.31 y (1.5) lo es de (1.21. 


Si tomamos la divergencia de los dos mimebros de (1.11 y tenemos en cuenta 
que el operador V • (Vx) es identicamente nulo, llegamos a: 
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sustituyendo en esta igualdad V • J, de acuerdo con (1.31 y conmutando 
el orden de los operadores V- y d/dt en el segundo termino del primer 
miembro, obtenemos: 

dp d 


o lo que es lo mismo: 


dt dt 


V-£> = 0 




que nos dice que V • T> — p no varia localmente en el transcurso del tiempo. 
Su valor actual en un punto del espacio debe ser igual al que tenia en el 
mismo punto en cualquier instante anterior, p. ej. hace un millon de anos. 
Es natural admitir que este valor era cero y, en consecuencia: 


V V = p 


Repitiendo el mismo razonamiento con (1.2) se llega a (1.5 I. 


Por ultimo, hay que destacar las siguientes observaciones: 


A. Los campos vectoriales y escalares que aparecen en las ecuaciones 
anteriores son reales, por su propia definition. 

B. En lo sucesivo, salvo indication expresa, consideraremos las tres cons- 
tantes materiales <r, e y p como magnitudes escalares independientes 
del punto del espacio y del tiempo, consecuencia de admitir que el 
medio es homogeneo e isotropo. 

C. Todas las leyes integrates del electromagnetismo experimental se de- 
ducen con todo rigor de las ecuaciones postuladas. 


1 . 1.2 


Condiciones de contorno 


Para resolver el sistema de ecuaciones en derivadas parciales (1.1)-(1.8) hay 
que ahadir: 


A. Condiciones de contorno 

B. Condiciones limite 

C. Condiciones iniciales 


que exige cualquiera de estos sistemas para su resolution. Vamos a limitarnos 
a las primeras, pues las segundas, ademas de ser las habituates, no nos seran 
necesarias, y las ultimas, como veremos mas adelante, seran sustituidas por 
la busqueda de soluciones estacionarias (Fourier). 

Las condiciones de contorno relacionan los valores que toman las magni¬ 
tudes del campo electromagnetico en dos puntos infinitamente proximos, 
situados a uno y otro lado de una superficie de discontinuidad del espacio 
soporte del campo, superficie que sera la frontera de dos medios materiales 
electricamente distintos, que no estan caracterizados por las mismas cons- 
tantes a, ey p (habitualmente conocidos como medios 1 y 2). 

Postulamos que las magnitudes que intervienen en cualquier problema de 
campo electromagnetico satisfacen las condiciones de contorno, sobre una 









1.1. El campo electromagnetico sobre el cuerpo real 


superficie de separation S de dos medios distintos (anexo 0: 

n-{V 2 - T>i) = Ps 

ft ■(&-&) =0 

h x (H 2 - Hi) = J s 
h x (£ 2 — E\ ) =0 

en que h es un vector unitario normal a S y dirigido del medio 1 hacia el 
medio 2. Veamos ahora que ocurre en dos casos particulares. 


1.1.3 _ 

Medio material de conductividad infinita. Condiciones de contorno en la 
frontera de otro medio 


Supongamos un medio conductor perfecto, en que a —> oo. La densidad de 
corriente J en un punto cualquiera de tal medio seguira valiendo: 

J = a£ 


y sera infinita si £ es distinto de cero. Como es inadmisible ffsicamente que 
J valga infinito en ningun punto, £ debe ser nulo en todos los puntos del 
medio. Pero entonces, (1.21 nos dice que B no varia en el transcurso del 
tiempo, y retrocediendo entonces en este como lo hicimos antes, deducimos 
que B debe ser tambien nulo. 


Como V = e£ y B = pH, y tanto e como p no pueden ser infinitos, no 
son posibles la creation ni la propagation del campo electromagnetico en 
medios de conductividad infinita. 


El razonamiento parece impecable, pero no lo es si H y T> son nulos en todos 
los puntos del medio, lo que nos dice que J tambien debe ser cero, o sea 
que no pueden existir corrientes de conduction en medios superconductores. 
Si queremos eliminar esta contradiction, debemos aceptar que dV/dt y, en 
consecuencia V, sean distintos de cero, lo que exige que e sea infinita; V 
queda indeterminada, en consecuencia. Como no existen medios materiales 
de valor infinito de /x, si B es cero H sera nulo. 


Consideremos ahora las condiciones de contorno en la frontera S de un 
medio 1, de condictividad infinita, y un medio 2, de conductividad finita. 
Aceptando que J y D no sean nulos, quedan indeterminados, por lo que 
habia que prescindir de las ecuaciones en que intervienen. Por otra parte 
£\ = B\='H\ — 0, B 2 = /.iH 2 y las condiciones (A.l |-(A.5) se reducen a: 


fi ■ B 2 = 0 h-H 2 = 0 (1.9) 

h x H 2 = J s (1.10) 

n x £ 2 = 0 (1.11) 


Observamos que: 


(1.91 afirma que debe ser nula la componente de H 2 segun la normal 
a S 


(1.11) establece que lo mismo debe ocurrir para la componente de £ 2 
segun el piano tangente a S 
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(1.101 da explfcitamente la corriente superficial en funcion de H 2 


Podemos expresar estos resultados de modo sencillo diciendo que en los 
puntos de un medio material de conductividad finita inmediatos a una su- 
perficie frontera con otro de conductividad infinita, el campo electrico debe 
ser normal y el campo magnetico tangencial a dicha superficie. Ademas, la 


corriente sobre este viene dada por (1.101, en que n es un vector unitario 


normal dirigido en el sentido que va del medio conductor perfecto al otro 
medio. 


1.1.4 _ 

Medio material de conductividad nula 


Si un medio material es un aislante perfecto, a es cero y (1.61 nos dice 
que no hay corrientes de conduccion J en ningun punto del medio. Pero 
entonces, (1.31 establece que la densidad cubica de carga no depende del 
tiempo, es decir, que se mantiene igual al valor inicial en todos los puntos 
del medio. Es un resultado que ffsicamente podiamos haber dado a priori, 
pues si el medio no es conductor, no permitira el desplazamiento de cargas 
en su interior y p se mantendra invariable. 


En los problemas que trataremos a continuacion supondremos, salvo en 
casos excepcionales, que es un aislante perfecto el medio para el que plan- 
tearemos y resolveremos las ecuaciones de Maxwell y, en consecuencia, pres- 
cindiremos de J. A pesar de que esta hipotesis solo implica que p sea in- 
dependiente del tiempo, pero no identicamente nula, postularemos tambien 
esta anulacion, apoyandonos en el resultado a que llegaremos en el apartado 
siguiente. Las ecuaciones locales para estos medios aislantes seran pues: 


„ - dV 

Vxn =m 

Y7 X dB 

Vx£= ~m 

V-P = 0 
V-B = 0 


V = eS 
8 = pH 


1.1.5 _ 

Evolucion temporal de la densidad cubica de carga 


Sea un medio material homogeneo e isotropo cualquiera con a finita y no 
nula, para el que son validas las ecuaciones (1.1)-(1.8). De (1.31 y (1.61, 
deducimos que: 


por lo que 


y de (1.41 y (1.71: 


V • J = V • (a£) = cjV • £ = - 


v i=--% 


V • V = V ■ (e£) =eV-£ = p 


dp 

dt 
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se sigue que: 

V-£ = - 
e 

Igualando los segundos miembros de estas ecuaciones, dado que el primer 
miembro es el mismo, llegamos a: 

1 dp p dp a 

a dt e = ’ > dt e ^ 

ecuacion que da la dependencia temporal de p y cuya integration inmediata 
es: 


p(x,y,z,t) = p 0 (x,y,z,t 0 )e "* = p 0 (x,y,z,t 0 )e ‘ 

siendo r = e/cr la constante de relajacion del medio, valor de t que corres- 
ponde a una caida de 1/e (63%). 

En esta solution, p 0 (x,y, z,t) es el valor de la densidad ciibica de carga en 
el punto ( x,y,z) en el instante initial t = 0. Pero este resultado nos dice 
que, salvo el caso extremo de un aislante rigurosamente perfecto, en que a 
sea identicamente nulo, la densidad cubica de carga tiende a anularse en el 
transcurso del tiempo, siendo tanto mas rapida esta tendencia cuanto menor 
sea la cte. r. P. ej. para el Cobre (conductor), r = 1.5 • 10~ 17 s y para el 
cuarzo (aislante) r = 20 dias. 

En consecuencia, salvo si en un instante determinado creamos voluntaria- 
mente una densidad cubica de carga, su evolution temporal sera continua- 
mente decreciente y esta diminution sera un proceso que se habra iniciado 
en el comienzo de los tiempos. Es pues natural suponer que, por pequeno 
que sea a, la densidad cubica de carga es despreciable y admitir que esta 
densidad es cero. 


1 . 2.1 _ 

Soluciones estacionarias 


1.2 EL CAMPO 
ELECTROMAGNETICO 
SOBRE EL CUERPO 
COMPLEJO (DOMINIO DE 
LA FRECUENCIA) 


De la misma forma que se hace en la teorfa de circuitos, cuando se trabaja 
en regimen permanente (soluciones estacionarias) es conveniente imponer 
variaciones de tipo armonico en el tiempo. En ese caso, tal y como ocurre 
en cualquier sistema lineal, excitaciones de entrada con la forma: 


x(t) = cos U)t 


proporcionan salidas del tipo: 

y(t) = Acos(ujt + <f>) 

Dado que el seno o el coseno no son mas que combinaciones de las expo- 
nenciales complejas y e~^ ut , es mas sencillo suponer que la serial de 
entrada es directamente una de estas funciones: 

X(w) = e jwt 

a la que corresponde una salida 

Y(u) = H(u>)X(u>) = Ae j ‘ t> ■ e jult 
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quedando el problema caracterizado por un numero complejo H = Ae^ 
que depende de w. Esta notation se conoce como fasorial. Entonces, para 
trasladar al dominio del tiempo cualquier fasor basta con multiplicar por 
ey tomar la parte real, p. ej. en el caso del fasor Y : 

y(t) = Re[H- A ut ) 


Las exponenciales e Jult son autofuncio- 
nes de los sistemas lineales y los autova- 
lores asociados son los H 


o la parte imaginaria si nuestra referencia es el seno. En cualquier caso, esto 
no es una limitation de cara a trabajar con senales de variation arbitra- 
ria, ya que estas se pueden expresar como suma de senos y cosenos en cut 
(transformada de Fourier). 

Todo lo visto hasta ahora para senales de una dimension, es directamente 
trasladable al calculo del campo EM en regimen permanente, ya que las 
ecuaciones de Maxwell son un sistema lineal, con la particularidad de que 
las magnitudes que manejamos son vectoriales, de forma que los fasores 
son tambien vectoriales. De esta forma podemos trabajar con campos de 
componentes complejas E que se pueden trasladar al dominio del tiempo 
sin mas que tomar: 

£ = Re[E ■ e jMt ] (1.12) 


Una ventaja de este planteamiento es que la variable t y las derivadas en 
t desaparecen, quedando la frecuencia u> como parametro que se conoce 
de antemano (especificado en cada caso). P. ej. si imponemos este tipo de 
solution a la ecuacion del rotational de £ y trabajamos en el cuerpo complejo 
(o dominio de la frecuencia), tendremos: 

V x Re[Ee jut ] = -^(Re[Be jut ]) 


Re[V x Ee jujt ] = i?e[-^(Be 3wt ]) = Re[-ju)tBe jut ] 

que, de forma compacta, se enscribe en el cuerpo complejo como: 

Vx£ = —juB 


El planteamiento presentado es equiva- 
lente a trabajar directamente con las 
transformadas de Fourier de las ecuacio¬ 
nes y los campos 


y tenemos una ecuacion que relaciona directamente los fasores complejos E 
y B a partir de los cuales se pueden calcular los campos en el dominio del 
tiempo £ y B aplicando ecuaciones del tipo (1.12). 


Entonces, el sistema de las ecuaciones de Maxwell junto con las relaciones 
costitutivas del medio y las condiciones de contorno expresado en el cuerpo 
complejo es: 


V x H = J + juiD 

(1.13) 

V x E = —jujjtH 

(1.14) 

V • D = 0 

(1.15) 

V-B = 0 

(1.16) 

J = crE 

(1.17) 

D — eE 

(1.18) 

B = fiH 

(1.19) 


(1.20) 
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junto con las condiciones de contorno: 


ft ■ (D 2 - -Di) = p s 

(1.21) 

ft ■ (B 2 - Bi) = 0 

(1.22) 

n • (J 2 - Ji) = -jurps 

(1.23) 

ft x 02 — H±) = Js 

(1.24) 

ft x (_E 2 — Ei) = 0 

(1.25) 


y resuelto el problema planteado por estas ecuaciones, se obtienen las versio- 
nes de los campos en el dominio del tiempo tomando la parte real o imagaria 
de cada vector multiplicado por e jaJt . 


Si existieran densidades de corrientes debidas a fuerzas no electromagneti- 
cas, y por tanto independientes del campo, la ecuacion (1.131 se escribiria: 


V X H — J ext J jcoD 


1 . 2.2 _ 

Medio material de conductividad infinita. Condiciones de contorno en la 
frontera con otro medio 


Las condiciones de contorno aplicables en este caso son: 


h ■ B 2 = 0 4=> n ■ &2 = 0 

ft x E 2 = 0 
n x H 2 = Js 


(1.26) 

(1.27) 

(1.28) 


Puesto que las partes reales e imaginarias de E 2 y H 2 satisfacen a (1.261 y 
(1.271, podemos traducir en palabras tales ecuaciones diciendo que en los 
puntos infinitamente proximos a un medio que sea conductor perfecto, los 
campos complejos electrico y magnetico deben ser respectivamente normal 
y tangencial a la superficie que limita dicho medio. 


1.2.3 _ 

Medio material de conductividad nula (aislante perfecto) 

Las ecuaciones correspondientes a un medio material de tal caracteristica 
se deducen del sistema anterior sin mas que anular la conductividad cr del 
medio en cuestion. 


1.2.4 _ 

Medio conductor no perfecto 

Es interesante hacer observar que si para un medio conductor, pero no 
perfecto, hacemos en el sistema de ecuaciones la sustitucion 

cr + jure = ju>e c 

las nuevas ecuaciones con esta nueva permitividad son las correspondientes 
a un medio aislante. En consecuencia, si el medio es conductor, podemos 
tratarlo como aislante de permitividad compleja e c (e c € C) mediante: 
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Notar que la parte imaginaria de e c es 
negativa 


Un material con perdidas es aquel en el 
que tan S ^ 0 


Un mismo medio material puede ser ais- 
lante a unas frecuencias y conductor a 
otras 


e c = 

que se suele expresar como: 

e c = e 

De esta forma, el sistema de las ecuaciones de Maxwell en medios ho- 
mogeneos, lineales y con perdidas se puede escribir tambien como: 

Vxff = juie c E 
Vx£ = —jujjiH 




LU 


<7 

- - 3 —, 

ue J 


junto con las condiciones de contorno, que pueden reagruparse como: 


n ■ [e c2 E 2 - e c iEi] = 0 

(1.29) 

n ■ [p 2 H 2 - piHi] = 0 

(1.30) 

h x [E 2 — Ei] = 0 

(1.31) 


en medios lineales y homogeneos. 

Como puede observarse, la parte imaginaria de la permitividad compleja 
esta relacionada con la conductividad y por lo tanto, con las corrientes de 
conduccion y las perdidas por efecto Joule. 

A frecuencias altas, por encima de los GHz, existe otro mecanismo por el 
cual pueden aparecer perdidas o lo que es lo mismo, parte imaginaria en 
e c . Este tiene que ver con calentamiento del material debido al amortigua- 
miento de las vibraciones de los dipolos e iones de que esta constituido 
(microscopico). Esto equivale a clecir que los vectores V y £ no son direc- 
tamente proporcionales (es decir, estan en fase), sino que aparece un cierto 
retardo 6 desfasaje entre ellos que es mas sencillo de expresar en el dominio 
de la frecuencia mediante una parte imaginaria e": 

D = {e'- je")E = [e\f) - je"(f))E = \7e' 2 + e" 2 e~i*E 


Para tener en cuenta este efecto se debe anadir un nuevo sumando a e c : 


e c = e - E 


. <7 



En la practica, las perdidas debidas a er y las debidas a e" son indistinguibles. 
Por ello, una cantidad de interes y ademas inedible, relacionada con esta 
expresion, es la conocida como tangente de perdidas electricas del material 
que se define como: 


tan 6 = 


eoe" + a 
ue' 


de forma que se verifica: 


e c = e(l — j tan J) 

que en funcion de la permitividad relativa e r resulta: 


e c = e r e 0 (l - jtanJ) 


donde e c es la permitividad del vacio. Finalmente, hay que clestacar que 
D = (e' - je")E ± e c E. 
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Por ultimo, hay que senalar que aunque con campos estaticos se pueden 
dividir los materiales en conductores y aislantes dependiendo de la cte. de 
relajacion r, no ocurre lo mismo en la electrodinamica. No obstante, es 
posible establecer un criterio, que dependera de la frecuencia, segiin el cual 
se considera (analizando el termino jue c E en la ecuacion del rotacional de 
H): 


■ Conductor : aquel medio en el cual predomina la corriente de conduc¬ 
tion sobre la de clesplazamiento (predomina la parte imaginaria en 
e c ) ; 

cr + we"(/) » ue'(f) 

■ Aislante (dielectrico): aquel medio en el cual predomina la corriente 
de desplazamiento sobre la de conduction (predomina la parte real en 

e c ): 

we'(Z) >> a + u)e"(f) 


1.3.1 _ 

Introduction 

El estudio electromagnetico se inicia con la electrostatica y sus clefinicio- 
nes de carga electrica y campo electrico, ligadas entre si por el concepto 
mecanico de fuerza. Dada una distribution estatica de cargas sobre conduc¬ 
tores, puede obtenerse energfa mecanica de esta distribucion por alteration 
de la misma o por desplazamiento de los conductores. En consecuencia, a 
la configuration estatica en cuestion le corresponde una energfa mecanica 
potential, de origen electrico, denominada energia electrostatica. 

Para evaluarla se extiende el principio de conservation de la energfa, que 
inclufa ya las energfa mecanica y calorffica, a este nuevo tipo y se calcula la 
correspondiente a una distribucion estatica de cargas, partiendo de la base 
de que sera igual a la energfa mecanica que hay que gastar para llegar a 
esa configuration clesde otra energfa nula, siempre que no hayan aparecido 
energfas de los tipos ya conocidos, mecanica y calorffica. 

Para hacer este calculo se elige arbitrariamente como configuration de energfa 
nula la que presenta ausencia total de cargas electricas a distancia finita de 
los conductores, estando tales cargas a distancia infinita. Esta afirmation 
no es objetable, dado que solo estaremos interesados en conocer diferencias 
de energfa y no sus valores absolutos. 

Pero si es discutible la hipotesis posterior, en la que se afirma que la energfa 
almacenada en el estado final, conseguido a base de traer las cargas clesde 
el infinito hasta los conductores, no depende del modo como se haga este 
aporte de cargas, siempre que se tomen las precauciones necesarias para que 
no aparezcan energfas indeseadas. 

Como no hay posibiliclad matematica de demostrar esta hipotesis, hay que 
introducirla en forma de postulado, que afirma en esencia que la energfa 
electrostatica del sistema es una energfa potencial. 

Para evitar la aparicion de energfas de otro tipo hay que traer las cargas 
en cantidades infinitamente pequenas, desplazandolas ademas de modo que 
pasen por una sucesion indefinida de estados de equilibrio, lo que requiere 


1.3 ENERGIA Y POTENCIA 
EN LOS CAMPOS 
ELECTROMAGNETICOS 
REALES Y COMPLEJOS. 
VECTOR DE POYNTING 
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un tiempo infinitamente grande y Se demuestra que esta energi'a es igual 
a una integral de volumen extendida a todos los puntos del espacio en que 
existan los vectores £ y T> 

I Udv 
■I v 

en la que aparece una densidad cubica de energia U que vale 


U = 


£-V 

2 


(1.32) 


si el medio es homogeneo e isotropo. 

La idea general de Maxwell consistio en postular que este resultado ma- 
tematico era la traduction de una realidad fisica, la del almacenamiento de 
la energia en el medio asilante exterior a los conductores, en los puntos del 
espacio donde existan £ y 2?, en lugar de estar concentrada en donde era 
logico suponer, en las cargas llevadas a los conductores. 

Esta idea, de la que se clesprendieron la definition de campo, la sustitucion 
de la accion a distancia entre cargas por la de accion del campo sobre la 
carga y la del retardo inherente al tiempo necesario para la propagation de 
este campo, demostro ser acertada al comprobarse dicho retardo. 

Pero nada autoriza matematicamente a afirmar que si la energia total es la 
integral del volumen anterior, el elemento de volumen dv contribuye con la 
cantidad UdV, o, lo que es igual, que en el elemento dv esta almacenada 
la energia UdV. Es la hipotesis mas sencilla, pero no la unica posible, ya 
que conocido el valor de una integral de volumen, hay infinitas cantidades 
subintegrales que conducen al mismo valor. Basta, por ejemplo, agregar a 
la integral anterior otra, extendida al mismo volumen y que valga cero, sin 
que sea identicamente nula la cantidad subintegral, para que tuviesemos 
distinta distribution local de la densidad de energia. Cumple esta condition 
la divergencia de cualquier campo vectorial de flujo total nulo a traves de 
las superficies de los conductores y de la superficie que delimite el campo a 
distancia infinita. 

Las mismas objecciones se presentan a proposito del razonamiento, analogo 
al anterior, que lleva a afirmar que la energia magnetostatica esta almace¬ 
nada en el medio con una densidad: 


W= H ^ 5 (1.33) 

si el medio es homogeneo e isotropo. 

Dificultades parecidas se oponen al rigor de la deduction de la formula 
que establece la densidad de energia por unidad de tiempo que se disipa en 
forma calorifica en un medio material conductor recorrido por una corriente 
estacionaria caracterizada por una densidad de corriente J 

q = £ • J (1.34) 


1.3.2 _ 

Energia y potencia en el campo electromagnetico real. Vector de Poyn- 
ting 

A pesar de que como ya hemos dicho antes, los razonamientos seguidos para 
deducir U, W y q se han hecho para fenomenos estacionarios y nada justifica 
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extenderlos a los variables, vamos a postular que tambien en este caso son 
validos los resultados obtenidos. Es decir, supondremos que, cualquiera que 
sean el campo existente y su variation temporal, la energfa electromagnetica 
esta almacenada en el medio, parte en forma de energfa electrostatica, con 


densidad (1.321 y el resto en magnetostatica, con densidad (1.33), teniendo 


ademas lugar la disipacion calorffica local (1.341 si el medio es conductor. 


Con estos postulados previos, pasaremos a presentar la interpretacion energeti- 
ca de una ecuacion local y de su correspondiente forma integral, ambas 
deducidas por Poynting a partir de las ecuaciones de Maxwell. 


Si de la igualdad deducida al multiplicar escalarmente (1.21 por H restamos 


miembro a miembro el resultado de multiplicar (1.11 por £, obtenemos: 


- - - QB - - - dD 

W(-V x£)~ £(-V xH) = -H~-£-J-£--^ 


Agrupando convenientemente los terminos y teniendo en cuenta la identidad 
vectorial: 


llegamos a 


\7-(AxB)=B-\7xA~A-\7xB 




(1.35) 


Multiplicando los dos miembros de esta igualdad por el elemento de volumen 
de dv e integrando el resultado para un volumen cualquiera V , llegamos a : 

J V -(£xH)dv + j £ ■ Jdv = — J £~dv- J H ■ ^dv 

Si el volumen V y su superficie exterior E cumplen las condiciones que per- 
miten aplicar el teorema de Gauss, podemos transformar la primera integral 
de volumen de la igualdad anterior en la integral sobre E del flhjo del vector 
(£ x H). Con esta transformation: 


V • {£ x H)ds+ / £ ■ Jdv = - 

£ JV 


V 


(L36) 


siendo ds el vector normal a E correspondiente al area elemental ds y diri- 
gido hacia el exterior de V. 


Debemos resaltar al llegar a este punto que (1.351 es rigurosamente valida 


en todos los puntos del espacio, si para ellos hemos postulado la validez 


de las ecuaciones de Maxwell, y que otro tanto le ocurre a (1.361, con la 


resticcion de que se satisfagan las condiciones necesarias para la aplicacion 
del teorema de Gauss al volumen V. 

Este rigor matematico va a desaparecer en la interpretacion energetica de 


(1.361 que vamos a hacer a continuation. La energfa electromagnetica W e - 


contenida en el interior de V en un instante cualquiera valdra, de acuerdo 
con la postulation hecha: 


W em = 


£ -V 


dv 


H-B 


dv 


iv 


y, asumiendo un medio homogeneo e isotropo, su incremento en un intervalo 
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de tiempo infinitamente pequeiio At sera: 


dW em A + d 

4tr - = TT A, = S 


£ ■ V 


dv 


H-B 


dv 


IV 


iv 


At = 





At 


si suponemos que V no varfa en el transcurso del tiempo. 
Puesto que: 



a as , a , , av 

e£ ' ei = £ Si (c£) = £ ' St 


y, analogamente, 

d - as 

Ft ^ 2 ) ~ n ' ~dt 

resultando finalmente: 


A W Km = 


[/■ >■ 


[*> 


At 


En consecuencia, el segundo miembro de (1.361 multiplicado por At, es 
igual a la variation en dicho intervalo de tiempo, cambiada de signo, de 
la energfa electromagnetica almacenada en el interior de V. Teniendo en 
cuenta el signo, sera igual a la energfa que ha desaparecido del interior de 
V en el mencionado intervalo de tiempo. 


Por otra parte, el segundo termino del primer miembro, tambien multiplica¬ 
do por At, da la energfa disipada en forma de calor en el volumen V durante 
el intervalo At. 


Ahora bien, si postulamos que en el caso de medios materiales en reposo 
las unicas formas posibles de energfa son la electromagnetica y la calorffica, 
la diferencia entre la desaparecida en la segunda sera una energfa que ha 
debido salir de V hacia otras regiones del espacio, dado que aceptamos 
la validez del principio de conservation de la misma, extendido al campo 
electromagnetico: 


AW Pm = - 


At (h (£ x 7?) • ds + At, £ ■ Jdv 


Entonces, lo mas natural e intuitivo es suponer que tal energfa haya salido a 


traves de la superficie E que limita V y como (1.361 nos dice que la diferencia 
en cuestion vale 


At y (£ x H) 


ds 


es decir, viene dada por una integral de superficie, es logico suponer que 
esta integral nos da el flujo de la energfa electromagnetica a traves de E por 
unidad de tiempo (notar que esta expresion es tanto mas aproximada a la 
verdad cuanto menor sea la unidad de tiempo elegida). 
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Puesto que la energfa que fluye a traves de E por unidad de tiempo es una 
integral de superficie, parece natural que el elemento ds de ella contribuya 
al flujo total con el valor de la cantidad superficial 

{£ x H) ■ ds 


que es el flujo a taves de ds del vector 

S = £xH 


(1.37) 


llamado Vector de Poynting. Afirmar que el flujo de S a traves de ds es 
igual a la energfa electromagnetica que pasa a traves de dicha superficie por 
unidad de tiempo (o, mas correctamente, a la potencia electromagnetica 
que pasa a traves de ds) implica que el flujo a traves de E sea la primera 
integral que aparece en ( |1.36 >. Pero a identico valor de esta intensidad se 
llega si en vez de tomar como vector de Poynting (1.371 tomamos: 


S' = {£ x H) + V x A 


(1.38) 


Hemos comprobado como la energfa 
electromagnetica se propaga. Mas ade- 
lante estudiaremos de que forma y a 
que velocidad ocurre este fenomeno 


siendo A un campo vectorial cualquiera. En efecto 


S' ■ ds = (j) (£ x Tt) ■ ds + (f) V x A ■ ds 
s Tt Tt. 


= (h {£ x H) ■ ds+ (b V • (V x A)dv = (b{£xH)-ds 
Tt Tv Tt 


Localmente, S' es distinto de S y dara lugar a diferente flujo de potencia 
para una superficie no cerrada. 

A pesar de todo, postularemos que el flujo de S = £ xH a traves de cualquier 
superficie, abierta o cerrada, es igual al flujo de potencia electromagnetica 
a traves de la misma, apoyandonos en dos circunstancias: 


■ es el vector mas sencillo 

■ hasta ahora, no ha llevado a ninguna contradiction con la experiencia 


1.3.3 _ 

Energfa y potencia en el campo electromagnetico complejo. Vector de 
Poynting complejo 

Aceptado el postulado anterior, si empleamos la formulation compleja, el 
flujo instantaneo de potencia a traves de una superficie E valdra: 

<j) [Re{Ee jut ) x Re{He jujt )} ■ ds 

cantidad que varfa periodicamente con el tiempo. 

En el estudio de las corrientes alternas se demuestra que no tiene interes 
practico, en general, el conocimiento del valor instantaneo de la potencia, 
siendo suficiente determinar su valor medio . 

De acuerdo con lo desarrollado en el anexo [B] el valor medio del vector 
de Poynting real, correspondiente a una solution compleja E c , H c de las 
ecuaciones de Maxwell, valdra: 

< S >= Re{Ee jujt ) x Re{He ju]t ) = ^Re(Ee jut x H*e~ jut ) = ^Re(Ex H*) 

(1.39) 
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El valor medio (1.391 es la parte real del llamado vector de Poynting com- 
plejo S: 


S=-ExH* 

2 


(1.40) 


Teniendo en cuenta la identidad vectorial compleja 

V • (A x B*) = B* ■ V x A - A* ■ V x B 
y recordando las dos primeras ecuaciones complejas de Maxwell 
V x H = aE + jueE => V x H* = aE* - jujeE* 


obtenemos 


V • S = • (E x H*) = H * ■ V x E ~ E ■ X7 x H*] = 

— aE ■ E* H —jueE ■ E* jupH ■ H* 

2 2 2 


(1.41) 


Analicemos los terminos del segundo miembro. Recordemos que en el campo 
real teniamos: 

■ q = £ ■ J = a£ ■ £= Densidad de disipacion de potencia en el medio 

■ U = l/2£ - V = l/2e£ • £= Densidad de energia electrica almacenada 
en el medio 

* W = 1/2H ■ B = l/2e7Y • H= Densidad de energia magnetica almace¬ 
nada en el medio 


Si hemos hallado la solucion E, H compleja de las ecuaciones de Maxwell, 
los valores que tienen en la realidad fisica q , U, W seran 

q =aRe(Ee-’ ut ) ■ Re(Ee jut ) 

U =^eRe(Ee jujt ) ■ Re(Ee jut ) 

W =^pRe{He juJt ) ■ Re(He jut ) 


Y teniendo en cuenta (B.l ) que a, e y p son constantes en el tiempo y que 
E ■ E* y H ■ H* son niimeros reales, obtenemos los valores medios: 


< q >= < aRe{Ee jult ) ■ Re(Ee jut ) >= ^aE ■ E* 
<U>=< ^ eRe(Ee jujt ) • Re{Ee jut ) >=\eE ■ E* 

< W >= < ^pRe(He jujt ) ■ Re{He juit ) >= ^/iH ■ H* 


Sustituyendo estos resultados en (1.411 


V • S = — < q > +2 jio[< U > 


< W >} 


Como < q >, < U >y<VE> son numeros reales, 

Re(V ■ S) =V ■ Re(S) = ~<q> 

Jw(V • S) =V • Im(S) = 2jiv[< U > - <W >] 


(1.42) 


(1.43) 

(1.44) 
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Si multiplicamos las ecuaciones (1.431 por el elemento de volumen dv, inte- 


gramos los result ados para un volumen V limit ado por una superficie E y 
aplicamos a las integrales de divergencia el teorema de Gauss, deducimos 


(1.45) 

<U > dv — / <W > dv (1-46) 

IV Jv 


Re 

j) S-ds 

-L 

Im 

j) S-ds 

Mi 


Las igualdades (1.451 dan a S categorfa de vector de uso practico, si bien 


la interpretacion es un poco mas complicada que en la version equivalente 
en el dominio del tiempo. Estas ecuaciones indican que: 


A. La parte real, cambiada de signo, del flujo de S’ a traves de una super¬ 
ficie cerrada E cualquiera y hacia su exterior es igual al valor medio de 
la potencia disipada en efecto Joule dentro de E. Si no existen fuen- 
tes de campo, la potencia neta que entra en V se disipa en el medio 
material. 

B. La parte imaginaria del flujo de S' a traves de una superficie cerrada 
E cualquiera y hacia su exterior es igual al doble del producto de la 
pulsacion del campo por la diferencia entre los valores medios de las 
energias electrica y magnetica almacenadas en el interior de E 


Resuelto asi en el campo complejo las cuestiones relativas a la energia y 
potencia en volumenes cerrados, queda todavia pendiente la del flujo de 
potencia a traves de una superficie, abierta o cerrada. Puesto que hemos 
postulado que la potencia que pasa a traves de una superficie abierta es igual 
al flujo a traves de ella del vector de Poynting real, si queremos calcular la 
citada potencia, cuando se ha planteado y resuelto un problema en el campo 
complejo, tendremos que utilizar: 


< P >= Re 

[ S-ds 

= 

/ Ex H* -ds 


[ Jn J 

2 

Un J 


Podemos resumir este resultado diciendo que, en el campo complejo, el valor 
medio de la potencia electromagnetica que pasa a traves de una superficie 
E, abierta o cerrada, es igual a la parte real del flujo del vector de Poynting 
complejo, S, a traves de dicha superficie. 
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1.4 EJERCICIOS 1. 14 El piano z = 0 de un sistema cartesiano de referenda constitu- 


ye la frontera entre dos medios dielectricos de permitividades y £2 res- 
pectivamente. Se mide el campo electrico sobre dicho piano en el lado 1 
obteniendose: 

Ei = Ei(x + y + z) 


A. ^Cuanto vale E 2 en el piano al otro lado de la discontinuidad? 

B. iQue densidad superficial de carga habria que depositar en la separa¬ 
tion de los dielectricos para que el campo tuviese el mismo valor en 
ambos medios? 


Sol: A) E 2 = Ei(x + y+ ^z); B) p s = (e 2 - ei)£i 


1. 15 | Expresar fasorialmente el campo £ = 2 cos (u)t + ir/A)x. 
Sol: A = 2e^/ 4 x 


1. 16 Calcular el vector D en el seno de un medio material con e' r = 10 y 


e" = 0.1 sabiendo que E = x. 

Sol : V = 10e o cos(wt — 9.9 • 10 _3 )a; 


1. 17 | Calcula el vector densidad de corriente lineal en la superficie de un 
conductor perfecto sobre la que incide de forma normal un campo electro- 
magnetico tal que H = H 0 {x + y). 

Sol : J s = H 0 (-x + y ) 


1. 18 


Calcular el vector de Poynting real y complejo del campo: 
E = E a cos (u>t — kz)x 


E a cos (uit — kz)y 


H = 


y la potencia que atravesara una superficie cuadrada de lm de lado situada 
en el piano z = 0. 

Sol: A) 5= ^E 2 0 [\ + icos2(o;t- kz)]z; B ) < S >= z ; C) 5 = 

I - I?2 

±E°. $• D1 P = - 1? 2 - 

2 z, u) -r . „ 


1. 19 


Calcular el vector de Poynting correspondiente al campo: 


E = E n 


a —jkr 


H = 


E n e-i kr 


asf como la potencia que fluye a traves de una superficie de radio R centrada 
en el origen de coordenadas. 


Sol: A) S = 




r- B) P = 


2tt 


E n 









CAPITULO 2 


Propagacion en medio indefinido 


Buscamos soluciones de las ecuaciones de Maxwell para campo complejo en 2.1 ONDAS PLANAS 
medios sin fuentes, homogeneos, isotropos, sin perdidas y con variaciones 
del tipo: 

E,H=f(r)e^ t 

Tomando rotacionales sobre las ecuaciones de Maxwell y simplificando, se 
llega a: 

A E + k 2 E = 0; AH + k 2 H = 0 

donde k = io^fep, se conoce como numero de ondas 6 constante de propaga¬ 
cion. 

Simplificamos el problema para el caso unidimensional (suponiendo que los 
campos solo varian con z y no con x e y). Tomamos arbitrariamente el eje 
x paralelo a E: 

E = E x x 

aplicando la ecuacion Vx£ = —jw J se llega a: 

H = H y y 

esto es, los campos estan en cuadratura espacial. Tambien se verifica que si 
E z = 0, entonces H z = 0 y se cumple que: 

E x H x z = 0 


Busquemos ahora soluciones de esta ecuacion de ondas, ecuacion que se 
convierte en 

dz z 

cuya solucion es de la forma: 

E = [E ie - jkz + E r e +jkz \x 
H = [H ie - jkz + H r e +Jkz \y 


siendo E l , E r , Hi y H r constantes complejas por determinar clependientes 
de las condiciones iniciales y de contorno. 

Si expresamos la solucion anterior en el dominio del tiempo resulta: 

£ = [Ei cos (uit — kz) + E r cos {ait + kz)]x 
H = [Hi cos(u>t — kz) + H r cos (u>t + kz)\y 
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donde hemos considerado que E i: , E r , Hi y H r son constantes reales. 

El termino Ei cos (u>t — kz) corresponde a una onda viajera ( travelling wave) 
que viaja hacia z + . El termino E r cos (cot+kz) es una onda viajera hacia z~. 
La combination de las dos ondas viajeras se denomina onda estacionaria. 
En el dominio de la frecuencia podemos identificar: 

■ e~l kz = Onda viajera positiva 

■ e^ kz = Onda viajera negativa 

Recordamos que la solution obtenida incluye todos los posibles terminos. 
Hasta que no fijemos las condiciones de contorno e iniciales, no quedara 
completamente determinada. Esto es equivalente a fijar el valor del campo 
en t 0 y z Q (o lo que es lo mismo, el argumento absoluto de los cosenos). 

2 . 1.1 _ 

Velocidad de fase 

Tomemos un solo sentido de propagation. Un piano de fase constante es un 
piano movil perpendicular a la direction de propagation: 

tot — kz = cte. — >_L z 

La velocidad a la cual se mueve un punto de fase fija se denomina velocidad 
de fase, Vf y viene clada por: 

dz d /cot — cte\ lo 1 

dt dt \ k ) k ^/Jie 

En el espacio libre: 

1 

Vf = Co = 

y/^oj^o 

donde c 0 es la velocidad de la luz en el vacfo y en cualquier otro medio: 

Co 

v f = , 

V trfJ'r 

La velocidad de fase en un medio mate- Se define la longitud de onda A como la distancia entre dos maximos o 

rial es siempre inferior a la velocidad de mmimos consecutivos. Para calcularla imponemos la condition: 
la luz en el vacfo 

(i uit — kz) — [cot — k(z + A)] = 27t 
de donde se sigue que: 

2t t _v f 

X -T~7 


2 . 1.2 _ 

Impedancia intrmseca 


Buscamos ahora una relacion entre las soluciones obtenidas para E y H. 
Aplicando la ecuacion del rotational obtenemos: 


- jkE x 


0E X 


-juonHy 


dz 




2.2. Ondas pianas en medios con perdidas 
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de donde se sigue que: 

H=-[Eie~ jkz - E r e jkz ]y (2.1) 

siendo 



que se conoce como impedancia intrinseca del medio (tiene unidades de f T) 
y en el vacio es: 

w 120 tt (fl) 

y en cualquier otro medio homogeneo y sin perdidas se escribe: 



V = Vo 






Notar que el sumando que corresponde a &* kz 
con respecto a la solution de E. 


en ( 2.11 lleva el signo cambiado 


Cuando trabajemos en medios con perdidas, debemos emplear la permitivi- 2.2 ONDAS PLANAS EN 
dad y la permeabilidad compleja con lo que la ecuacion de ondas resulta: MEDIOS CON PERDIDAS 

AE-^E; AH-'y 2 H = 0 


donde 

7 = jvy/n c e c 

se denomina constante de propagacion compleja que normalmente se escri- En la rafz compleja que define 7 nos que- 
|^ 0 . damos con el valor principal 

7 = a + j/3 

En este caso, la solution del campo tiene la forma: 

E = \E i e~ z/z + E r e + ~< z }x 
H = [H ie ~ lz + H r e +JZ ]y 

La onda viajera positiva tiene un factor de propagacion: 

e -7* = e -« . e -j0x 

que en el dominio del tiempo tiene la forma: 

e~ az cos (cot — (3z ) 

y observamos como, debido a las perdidas, aparece un factor de amorti- 
guamiento exponential que afecta al modulo de la solution. Por todo ello, 
denominamos: 


■ a= constante de atenuacion (to x ) 

■ / 3= constante de fase (to -1 ) 

Por ultimo, el campo H expresado en terminos de E sera: 

H = —[Eie~ JZ - E r e lz }y 
V 


Las unidades mas frecuentes de a son 
Nep/m 6 db/m y las de /3, rad/m. Un 
neperio es equivalente a 8.68 dB 
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siendo: 

/Tti ju^c 
77 V £c 7 

la impedancia intrinseca compleja del medio con perdidas que podemos es- 
cribir como: 

7 = ju)^p c e c = jaj^/prVoe r eo\/l - j tan 5 ( 2 . 2 ) 

Estudiamos ahora dos casos particulares: 

■ Medio real con bajas perdidas (buen dielectrico) 

■ Medio real buen conductor 


2 . 2.1 _ 

Ondas planas en medio buen dielectrico 


Supongamos un medio buen dielectrico (aislante) con bajas perdidas. Si 
hacemos un desarrollo en serie de Taylor de ( |2.2| ) para el caso tan S « 
llegamos a: 


7 ~ juy/p,e r e 0 



1 

2 


j tan S + ... 


de forma que: 


V 

a 


L0^/pe r e o 
A: tan <5 
2 


resultado que podemos resumir diciendo que: 


■ La constante de propagation es aproximadamente la de un medio sin 
perdidas con la misma permitividad y permeabilidad relativa. 

■ La constante de atenuacion depende directamente de las perdidas (de 
tan S ). 


2 . 2.2 _ 

Ondas planas en medio buen conductor 

Un medio buen conductor es aquel en el que a » toe esto es aquel en el 
que las corrientes de conduction son mucho mayores que las corrientes de 
desplazamiento. 

En este caso, despreciando el termino de la corriente de desplazamiento, 
podemos escribir: 


8 S es la distancia que hay que recorrer 
para que la amplitud del campo caiga 

1/e = 36.8% 


de forma que se cumple: 


a = /3 = 


ujficr 


Definimos 5 S = profundidad de penetracion o espesor de piel como: 


S.= - = 

a 


ujfia 
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La impedancia intrfnseca compleja en el seno de un buen conductor sera: 


V = 





(1 + 3) 


1 

aS s 


Notar que la fase de rj es 45°, lo cual es caracteristico de los buenos conduc- 
tores. 


Resulta sencillo comprobar que el valor medio del vector de Poynting para 
una onda viajera positiva es: 


< S>= \^°L. e ~ 2az z 

2 if 


siendo \E a \ el modulo del fasor del campo electrico. Se observa que: 


■ < S > apunta segun la direction de propagation, como era de esperar 

■ < S > se atenua segun un factor exponencial e~ 2az 


En el caso de un medio sin perdidas (con y £ R) esta expresion se convierte 
en: 


< S>= 


lM 

2 V 


z 


vector que tiene solo parte real. 


Generalizamos los resultados obtenidos para una direction de propagation 
cualquiera h. La ecuacion de ondas se escribe: 


A E + k 2 E = 0 => 


d 2 E 

dx 2 


d 2 E 
dy 2 


d 2 E 

Ih 2 


+ k 2 E = 0 


cuya solution es: 

E = E 0 e- j%p 

siendo: 


k = kh = k x x + k v y + k z z 
f = xx + yy + zz 


el conocido como vector de propagacion cumpliendose ademas que: 


k ■ E a = 0 

E x H x k = 0 

- 1^-1 
H = —k x E = —n x E 
V 1 

8 = E 0 cos(c ot -k-r) 
<S>= &- e - 2ak - p k 

2y 


2.3 VECTOR DE POYNTING 


2.4 ONDAS PLANAS. 

PLANTEAMIENTO 

GENERAL 


k = h es la direccion de propagacion 


etc. 
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2.4.1 _ 

Observaciones sobre las ondas planas 

Notar que: 

■ En cualquier punto de fase constante las soluciones de E y H son las 
mismas. 

■ Es frecuente pintar rayos que siguen la direction de propagation, tal 
y como se hace en optica geometrica. 

■ Matematicamente, el conjunto e^ kz es un conjunto completo de las 
soluciones de las ecuaciones de Maxwell. En la tabla 12.11 se resumen 
las soluciones de la ecuacion de ondas para un sistema de coordenadas 
rectangular. 


Tipo 

Funcion de onda 

Ceros 

Infinitos 

Viajera 

e~iP' z para +z 

0z- 

-j 00 

0 z- 

-*■ joo 


e J @ z para 

0z- 

j 00 

P z- 

-*■ -JOO 

Estacionaria 

cos (/3z) para ±z 

(3z = 

±(n + 1/2)7 t 

0z- 

->■ ±jo0 


sin(flz) para ±z 

Pz = 

±mr 

pz- 

->■ ±i 00 

Evanescente 

e~ az para +z 

az — 

OO 

az - 

-> —OO 


e az para — z 

OLZ — 

— OO 

az - 

-*■ +00 


cosh(az) para ±z 

az = 

± j(n + l/2)7r 

az - 

-*■ ±00 


sinh(a:z) para ±z 

az = 

dzjnn 

az - 

-> ±00 

Viajera 

e ~3~t z = e ~az e -j0z p ara 

7 z — 

OO 

7 z - 

-» —OO 

atenuada 

e 3iz = ^z^Hz para 

7 z — 

—OO 

7 z - 

-> OO 

Estacionaria 

atenuada 

cosh (72) = 
cos(ctz) cosh(/3z) 

-jsin(az) sinh(/3z) para ±z 

7 z = 

±j(n + 1/2)7 t 

7 z - 

■+ ±j 00 


sinh(7z) = 
sin(az) cosh(/3z) 

-jcos(az) sinh(/3z) para ±z 

7 z = 

±jnir 

7 z - 

■+ ±j OO 


Cuadro 2.1: Funciones de onda, ceros e infinitos para funciones de onda planas en 
coordenadas rectangulares 


■ Existen otros conjuntos completes de soluciones que se pueden em- 
plear dependiendo del sistema de coordenadas apropiado a las condi- 
ciones de simetria del problema. En los apendices ?? y ?? se describen 
los correspondientes a las coordenadas cilmdricas y esfericas. 


2.5.1 


2.5 DISPERSION 
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Velocidad de grupo 

Tal y como hemos visto, la velocidad de fase viene dada por: 

dz lo 
Vf= dt = k 


Supongamos una serial muy simple = rp! + 4' 2 donde: 

d'l = cos (kz — Lot) 

^2 = cos[(/c + A k)z — (lo + A Lo)t] 


Entonces: 

>3/ = 2 cos 


zAk — tAco 
-2- C ° S 





t 


Observamos que la amplitud de la serial A es variable y que tenemos una 
const ante de propagation efectiva de valor k+Ak/2 y una pulsation efectiva 
de valor lo + Auo/2. 

Entonces, los puntos en los que tenemos una amplitud constante cumplen 
zAk — tAuo = cte. Podemos defrnir entonces una velocidad de grupo como: 

A LO 


que en el lrmite se expresa: 



Distinguimos dos tipos de medio: 

A. Medio no dispersivo: cuando Vf = v g , o lo que es lo mismo, (3 = k = 
ujy/JIi, esto es, la dependiencia es lineal. En este caso, las componentes 
de distintas frecuencias viajaran a la misma velocidad por lo que se 
conservara la forma de onda y por lo tanto no existira distorsion. 

B. Medio dispersivo: cuando vf ^ v g , y por lo tanto, exista distorsion 

Normalmente se verifrca v g < Vf caso en el que diremos existe dispersion 
normal y v g sera la velocidad a la que viaja la energra. Cuando v g > Vf 
diremos que estamos ante dispersion anomala. Este tipo de dispersion ocurre 
cerca de las frecuencias de resonancia de las distribuciones dipolares que 
forman los medios materiales. 


2.5.2 _ 

Diagrama de Brillouin 


Para representar la dispersion de un medio se utiliza el Diagrama de Bri¬ 
llouin, que no es mas que una representation de u vs. (3 (frgura 3.2). 


En un medio dispersivo general tendremos que Vf ^ v g , siendo: 


Vf = tan 4>i 
v g = tan 4> 2 
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CO 



> 

P 


Figura 2.1: Diagrama de Brioullin 


Si e, p ^ f(io) la serial se propagara sin distorsion. Si el medio tiene perdidas 
tendremos: 

7 = - jtan<5 


y si <5 <<: 


P 



6 

2 


y como se verifica: 

5 = arctan — = f(u>) 

Loe 

se sigue que todos los medios reales son dispersivos. Si ademas <5 >> la 
serial experimentara distorsion y atenuacion, con lo cual sera irreconocible 
a una cierta distancia. No obstante, si el medio es de bajas perdidas, apenas 
tendremos dispersion. 


2.6 POLARIZACION La polarizacion de una onda plana se refiere a la orientation del campo 

electrico, que puede seguir una direccion fija o cambiar con el tiempo. 

Si suponemos que la direccion de propagation es z y fijamos un instante t 
tendremos: 

£ = £ x x + £ v y 


siendo 


£ x = £\ cos ut 
£ y = £ 2 cos (lot + S ) 


Si fijamos un piano z= cte. y observamos lo que ocurre a lo largo del tiempo, 
en el caso mas general, el extremo del vector £ describira una elipse (a;, y) 
como la de la figura |2.2| 

Buscamos la ecuacion de esta elipse. Si en el sistema de ecuaciones anterior 
eliminamos tot se llega a: 

£fy 2 + £ 2 x 2 — 2£i£ 2 cos Sxy — £ 2 £ 2 sin 2 <5 = 0 

ecuacion de la elipse de polarizacion. 

Veamos varios casos particulares: 

A. Si 5 = 0 (componentes en fase): 

£2 

V = 
y £1 

la polarizacion es lineal (rectilinea). Si 5 = n entonces y = ~^x. 
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B. Si d = ±7r/2 (componentes en cuadratura) y £\ = £2 = £'■ 

2,2 c2 

y + x = t 

la polarizacion es circular. 

C. En el resto de los casos, la polarizacion es eliptica. 

Notar que las polarizaciones eli'pticas y circulares pueden girar a derechas o 
a izquierdas , dependiendo del sentido del giro del extremo del vector £. Si 
el vector de propagation apunta segun la perpendicular del piano del papel 
hacia arriba, hablaremos de: 

■ Polarizacion a derechas o positiva: cuando el sentido de giro es anti- 
horario 

■ Polarizacion a izquierdas o negativa: cuando el sentido de giro es el 
horario 

En lo visto hasta ahora, hemos expresado cualquier polarizacion como una 
combination de dos polarizaciones lineales segun los ejes x e y. Dos polari¬ 
zaciones lineales ortogonales forman una base para expresar cualquier otra 
polarizacion. De la misma forma, dos polarizaciones circulares de signos 
contrarios tambien son base (ver figura [2tijl . de forma que a: 



Figura 2.3: Base de polarizaciones circulares 
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£ x = E\ cos cot 
£ y = £2 cos (uit + 5) 


corresponde: 


siendo S' el desfasaje existente entre las dos polarizaciones circulares. 

2 . 6.1 _ 

Representation de estados de polarizacion 

Supongamos que elegimos un sistema de coordenadas de forma que el campo 
E se expresa como: 


—= [£i cos cot + £ r cos (cot + £')] 

V 2 

—= [—£1 sin cot + £ r sin(o;f + <5')] 

V 2 


£ x = £\ cos (cut) 

£ y = £\ cos (uit + S ) 

campo al que corresponde la elipse de polarizacion es la de la figura |2.4| 



Definimos: 

R = Razon axial, siendo \R\ = OA/OB, tomando R > 0 para la polariza¬ 
cion a derechas y R < 0 para la polarizacion a izquerdas 

e = Angulo de elipticidad (— 7 r /4 < e < 7t/4) 

r = Angulo de inclinacion (0 < r < it) 

7 = Angulo de aspecto (0 < 7 < tt/2) 

La razon axial cumple R\ > 1 

Las representaciones mas habituales de los estados de polarizacion y sus 
principales ventajas son: 
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■ Elipse cle polarizacion: parametros (e, r) 6 ( 7 , 5). La elipse es facil- 
mente construible. 

■ Vector complejo de polarizacion: Facil cle calcular. Preserva la fase en 
la interaccion onda-antena. 

■ Coeficiente {ratio) de polarizacion: pg (complejo). Util para represen- 
tar medios depolarizantes. 

■ Esfera de Poincare. Todos los estados de polarizacion se representan 
en una unica esfera. 

■ Parametros de Stokes: Si, s 2 > S 3 (reales). Facilita la evaluation de la 
interaccion onda-antena. 

Nos centraremos en los cuatro primeros. 

Parametros de la elipse de polarizacion 

En las expresiones que relacionan (e, r) y 
( 7 , 5) aparecen funciones trigonometri- 
cas inversas en las cuales hay que respe- 
tar el signo del cuadrante de la solucion 


Para las representaciones basadas en la elipse de polarizacion, las relaciones 
entre los distintos parametros son: 


e = — arcsin(sin 2j sin 5) 


7r 7r 

-< e < - 

4 — — 4 


junto con: 


r = - arctan(tan 2 q cos S) 


0 < T < 7T 


7 = - arc cos(cos 2 e cos 2t) 


5 = arctan 


tan 2 e 
sin 2 r 


°< 7 <| 


Vector de polarizacion 

Supongamos un campo de la forma: 

£ = £ 1 cos u>tx + £2 cos (ujt + 6)y 
que en el dominio de la frecuencia tiene como fasor complejo: 

E = £\x + £2 e^y 

Definimos el vector de polarizacion como un vector unitario e que tiene la 
direccion de E. Como ademas se cumple: 

. £2 

tan 7 = -- 
£1 

resultara: 

e = cos 7 x + sin 7 e? s y 


Coeficiente (ratio) de polarizacion p 7, 


Es un numero complejo que se define como: 


PL = 


Ey 

E„ 


—e jS 

£1 ' 


Notar que e es un vector complejo uni¬ 
tario y que incluye toda la informacion 
de polarizacion 
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2.7 ONDAS PLANAS EN 
CAMBIOS DE MEDIO 


Suponemos dos medios homogeneos e isotropos caracterizados por (ei, /r 1; ) 

y (e 2 ) IJ- 2 , 02 ) que tienen como superficie de separacion el piano S (dioptrio). 
Podemos plantear dos casos: 

A. La superficie de separacion es infinita (o muy grande en comparacion 
con A). Postulamos una solucion que consiste en tres ondas planas: 
incidente, reflejada y transmitida 

B. La superficie de separacion es del orden de A: ocurre un fenomeno de 
difraccion. Este fenomeno se estudiara posteriormente. 

En este apartado estudiaremos el primero de los casos analizando dos situa- 
ciones: 

A. El campo E incide normalmente a S (direccion de propagacion k per¬ 
pendicular a S) 

B. El campo E incide de forma oblicua a S (direccion de propagacion k 
oblfcua a S)) 

El procedimiento en ambos casos consiste en postular unas soluciones que 
relacionamos mediante las condiciones de contorno en la superficie de dis- 
continuidad. 


2.7.1 _ 

Incidencia normal 

Suponemos, sin perdida de generalidad, que el campo electrico incide sobre 
S orientado segun x. Postulamos la existencia de: 

■ Una onda reflejada en el medio 1 (z < 0) 

■ Una onda transmitida (o refractada) en el medio 2 (z > 0) 



Figura 2.5: Incidencia normal 
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La onda que excitamos en el medio 1 sera: 

Ei = E a e~ llZ x 
Hi = — e~^ z y 

m 

y definimos las ondas rcflejada: 

E r = E 0 Te llZ x 

H r = ——Te 1lZ y 
Vi 

y transmitida: 

E t = E 0 Te~ l2Z x 

H t = —Te~ l2Z y 
Vi 

donde falta por encontrar los niimeros complejos T y T conocidos como: 

■ r= coeficiente de reflexion del campo electrico en la discontinuidad 

■ T= coeficiente de transmision del campo electrico en la discontinuidad 

r y T relational! los fasores de las ondas rcflejadas y transmitidas con el de 
la onda incidente en la discontinuidad. 

Si ahora aplicamos las condiciones de contorno en dicha discontinuidad: 
n x (E 2 — Ei) = 0 y fix (H 2 — Hi) = 0 obtendremos: 


de donde se sigue que: 


T + 1=T 


i -r t 


CN 

1 


r _ V2 - m 
m + m 

(2.3) 

T _ 2? ?2 
r )2 + m 

(2.4) 


Los campos en el medio 1 seran: 

E 1 = E i + E r = E 0 e~ llZ (l + re 27lZ )f 

= Hi + H r = — e“ 7lZ (l - Te 2 ' 1lZ )y 
Vi 


Podemos introducir un coeficiente de reflexion generalizado para cualquier 
punto z, clefinido por: 

r(«) = r e 27iz = r(o)e 27iz 

de forma que resulta: 

Ei = E 0 e~^ z [l + T{z)}x 
Hi = — e -^ z [l-T{z)]y 

m 

expresiones correspondientes a una onda estacionaria. Pasamos a estudiar 
dos casos particulares. 


Notar que podi'amos haber definido unos 
coeficientes de reflexion y transmision 
analogos para el campo magnetico 


Se pueden definir unos coeficientes de re¬ 
flexion y transmision en potencia equi- 
valentes a los definidos para el campo 
electrico 
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Notar que la distancia entre maximos (o 
mfnimos) consecutivos es de Ai/2 y en¬ 
tre un maximo y un mfnimo consecutivo, 
Ai/4 


Medios 1 y 2 sin perdidas 

Suponemos que 01=02 = 0. Entonces tendremos: 

71 = 3 Pi = jkoy/Pr l£rl 

72 = jP. 2 = jkoy/Pr 2ti2 

siendo ademas 771 y 772 € 7?., por lo que T(0) sera tambien un numero real 
(con signo). El modulo del campo electrico sera: 

\e x \ = |£ o ||i + r(0)e 2 ^ i2 | = \E 0 \yJi + r(0) 2 + 2r(o) cos 2 k lZ ( 2 . 5 ) 

La figura [2~6| representa esta ecuacion, en lo que se conoce como diagrama 
de onda estacionaria 


El 


.. r 

X 

* - 

/2 

-►: 

A 


1 1 

■^max 

E min 


◄— ► 

< -► 

N 

II 

O 

1 


A -M- 

A.,/4 \!2 


Figura 2.6: Diagrama de onda estacionaria 


Es interesante introducir un nuevo parametro s, conocido como relation de 
onda estacionaria 6 ROE, dado por el cociente entre el maximo y el mmimo 
valor del campo electrico: 


ROE =s= 


E„ 


E„ 


1+|E| 

i — |r| 


que se puede obtener facilmente a partir de (2.51 y que toma valores en el 
intervalo 1 < s < oo. 


La ROE es un parametro facilmente medible (basta un medidor de campo) 
que proporciona directamente el modulo del coeficiente de reflexion en la 
discontinuidad: 


l r (o)l 


s — 1 
s + 1 


En el caso particular en que iji = r/ 2 , obtendremos T(0) = 0, s = 1 y 
|S(^)| = E 0 = cte. como corresponde a una onda viajera positiva en el medio 
1 (no existe onda viajera negativa ni, por tanto, onda estacionaria). 

Calculamos ahora el vector de Poynting en los dos medios: 


Si = \Ei x H* = ^^(e~ jklZ + r(0)e- jklZ )(e~ jklZ + T(0)e- jklZ )*z 

2 2771 


2771 


( 1 - | r ( 0)| 2 + 2jT(0) sm2k!z)z 


So = -Eo x H* 


\Eo\ 2 

2772 


\T\ 2 z 


( 2 . 6 ) 
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Observamos como Si tiene parte real (energia que se propaga) y parte ima- 
ginaria (elergia reactiva almacenada) mientras que S 2 tiene solo parte real. 
Se puede demostrar facilmente que las partes reales de los dos vectores cal- 
culados coinciden (lo cual es logico, ya que la energia que se propaga en la 
direction z desde el medio 1 pasa al medio 2). 


Medios 1 sin perdidas, medio 2 con perdidas 

En este caso, el diagrama de onda estacionaria vendra dado por: 

\E X \ = |£ o ||l + r(0)e 2jfel2 | = \EoWl + |r(0)| 2 + 2|r(0)|cos(2& 1 « + <t>) 
donde se ha tornado un T(0) complejo, dado por: 

r(o) = |r(o)|e^ 

Se observa como aparece un desfasaje adicional <fi en el coseno dado por la 
fase del coeficiente de reflexion en la discontinuidad. Senalar que el concepto 
de ROE sigue siendo perfectamente aplicable. 


Medio 1 sin perdidas, medio 2 conductor perfecto 

En este caso, no existe onda transmitida por lo que tendremos: 


T = 0 

r = -1 


El modulo del campo tiene la forma: 

\E \| = \E 0 \y/l + (—l) 2 — 2cos2fci0 = 21 E 0 \\ sinfcizj 


que corresponde a una onda estacionaria pura. En este caso los valores 
minimos del campo electrico son nulos, por lo que tendremos s —> 00 . Notar 
que el campo se debe anular en la discontinuidad (ver figura 2.7). En este 



Figura 2.7: Reflexion sobre conductor perfecto 

caso, la parte real del vector de Poynting Si es cero, como corresponde a 
una situacion en la cual no hay propagation de energia segiin z. 
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Caso general: Medios 1 y 2 con perdidas 

En este caso r(^) es un niimero complejo cuyo modulo depende tambien de 
z: 

r(z) = r(0)e 27lZ = r(0)e 2ai V i/3lZ 

y tanto la onda incidente como la reflejada se ven atenuadas. En la figura 
|2.8| se representa el ’’diagrama de onda estacionaria”que corresponde a esta 
situacion. Se observa como ya no se puede hablar de una ROE ya que los 
valores maximos y mi'nimo de campo electrico no se mantienen constantes 
al variar z. 



Figura 2.8: Reflexion (caso general) 


2.7.2 _ 

Incidencia normal con cambio de medio 

Suponemos una situacion en la cual tenemos presentes multiples medios, 
cuyas superficies de separation son pianos perpendiculares a la direction de 
propagation z. 



Figura 2.9: Incidencia normal con multiples cambios de medio 


Si tomamos la superficie que corresponde al cambio entre los medios n y 
n + 1 (donde ahora existira onda transmitida y reflejada procedente del 
interfaz con el medio n + 2), podemos calcular el coeficiente de reflexion 
T n {z n ) aplicando las condiciones de contorno, de forma que se cumple: 


S:+ r n+ iW 
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Como se puede ver, el planteamiento del problema es complicado cuando se 
usan coeficientes de reflexion y que en cada discontinuidad hay que recal- 
cular dicho coeficiente segun la formula anterior. 


Para tratar este problema, es conveniente introducir una impedancia defi- 
nida como: 


Z n {z) 


z x E n (z ) 
H n (z) 


conocida como impedancia en un piano de z constante que tiene dimensiones 
de ft, siendo E n {z) y H n (z) los campos totales (incidente+reflejado). En 
la ecuacion anterior tomamos un signo positivo para la impedancia (que 
relacionamos con la direccion i + ). Un signo negativo tambien hubiera sido 
valido. 


La definition de Z n (z) no es un cociente 
de vectores, sino de sus medidas (com- 
plejas) ya que numerador y denominador 
son dos vectores para Ielos 


La ventaja de esta nueva impedancia estriba en que relaciona los Et y Ht , 
campos que son continuos (incluido en las discontinuidades), por lo que 
cualquier relation entre ellos, tambien lo sera. 

Operando en la definition de la impedancia se sigue que: 


Zn(Zn) Tj n 


i + r n (z) 

i-r n { z ) 


_ Z n (z ) - Tj n 
Z n (z) + r] n 


Observamos que si r n = 0 (ausencia de onda reflejada), entonces Z n (z) = 
y n . Sustituyendo el valor de los coeficientes de reflexion, 


Z„(z) = r] n 


i + r n (z')e 2 ^( z - z ') 

1 - r n (z')e 2 ^(- z ~ z 'l 


que puede expresarse como: 


Z n (z')ch'y n (z - z') - p n shg n {z - z') 
rjnCh'Yntz - z') - Z n (z')shry n (z - z') 


y en el caso de un medio sin perdidas: 


Z n {z) = Vi 


Z n {z') cos @ n (z - z') - jy n sm.l3 n (z - z') 
Vn cos f3 n (z - z') - jZ n (z') sin j3 n {z - z') 


Por ejemplo, para una incidencia normal con dos cambios de medio, tal 
y como se indica en la figura 2.10 y aplicando la continuidad de Z(z), 
tendriamos: 


Pi(0) = 


^l( 0 )- 7 ?l 

-Zi(o) + Vi 

_ Z 2 (d) cos p 2 d + ji 72 sin fad 

1 ^ 2 V 2 cos (3 2 d + jZ 2 (d) sin [3 2 d 

Z 2 (d) = Z 3 (d) = V3 


2.7.3 _ 

Incidencia oblicua 

Suponemos ahora un caso en el cual la direccion de propagation de la onda 
plana incidente no coincide con la perpendicular a la superficie S. Situamos 
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d 

Figura 2.10: Incidencia normal con dos cambios de medio 



el origen de coordenadas en cualquier punto de dicha superficie, cuya normal 
viene dada por el vector h. A1 igual que hicimos en el caso de la incidencia 
normal, postulamos la existencia de tres ondas: 

■ Una onda inidente: Ei segiin la direccion hi 

■ Una onda reflcjada: E r segiin la direccion h r 

■ Una onda transmitida: E t segiin la direccion ht 

cuyas expresiones son: 

Ei 
Hi 

y definimos las ondas reflejadas: 

E r 
H r 


= E 0 e~ llfliJf! 

= —hi x Ei 
Vi 

= E 2 e- llflr ' ? 

= —hr X E r 

Vi 
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y transmitida: 

E t = 

1 

H t = — n t x E t 

m 


Leyes de Snell 

Imponemos que los exponentes de las ondas sean los mismos en la superficie 
S, definida por la ecuacion n ■ f = 0. Dado que podemos expresar r como 
(ver Stratton) 

f = (n • r)h — h x {hx r) 
sobre la superficie S sera: 


r = —h x (n x r) 

con lo que podremos escribir: 

7 ifb • h x (n x r) = 71 n r ■ hx {fix r) 

'yihi ■ h x {n xr)= 7 2 n r ■ hx {hx r) 

y dado que: 

hi ■ h x {h x r) = {hi x h) • {h x r) 

resulta: 

{hi x h — h r x h) • h x r = 0 
(71 hi x h — 72 ht x h) ■ h x r = 0 

de donde podemos extraer las siguientes conclusiones: 

■ Los vectores n, hi, h r y ht. estan en el mismo piano, que denominare- 
mos piano de incidencia, que es perpendicular a S. 

■ Primera Ley de Snell: 

sin 9 r = sin 0i =>■ 9 r = 9i 


■ Segunda Ley de Snell: 


71 sin 9i = 72 sin 9 t 


que en un caso sin perdidas resulta: 

y/fhJi sin 9i = y/n 2 e 2 sin 9 t 


Ecuaciones de Fresnel 

Aplicamos ahora las condiciones de contorno para determinar la relation 
entre los campos en cualquier punto de la superficie S: 

h x {Ei + E r ) = h x E t ", h x {Hi + H r ) = h x H t 


En los visto hasta ahora, la orientacion de Ei era arbitraria. Distinguimos 
dos casos: 





66 PROPAGACION EN MEDIO INDEFINIDO 


A. Ei es perpendicular al piano de incidencia, E % = E± 

B. E t es paralelo al piano de incidencia, = E\\ 

de forma, que cualquier campo puede descomponerse en esta nueva base 
como: 

E = E\\+ E± 


Caso 1: E, perpendicular al piano de incidencia (Ei = E±) En este 
caso se verifica h ■ E± = 0. Tendremos ademas que n ■ E t = h ■ E r = 0. Si 
aplicamos las condiciones de contorno llegamos a: 

? 7 2 cos 6i — r]i cos 9 t 
ip cos 9i + r) i cos 9 t 

2r]2 cos 9., 

r /2 cos 9i + r]i cos 9 t 

Si aplicamos la segunda Ley de Snell, obtenemos que para el angulo 9 t se 
cumple: 



que en un caso sin perdidas se convierte en: 




y vemos que este valor puede llegar a ser complejo, incluso aunque tratemos 
con medios dielectricos sin perdidas. No obstante, si sustituimos este valor 
en las ecuaciones anteriores utilizando la segunda Ley de Snell, solventamos 
de momento el problema y obtenemos una version mas practica en funcion 
unicamente del angulo de incidencia 9,: 


V 2 cos 9i - 771 \j 1 - (^) sin 2 9 r 
r/2 cos 9i + 771 sin 2 9 t 

2r]2 cos 9i 

?? 2 cos 9i + 771 sin 2 9i 

y podemos decir que los coeficientes complejos que multiplican a E]_ en cada 
caso indican que E' ± y E no estan en fase con E]_. 

Caso 2: E t paralelo al piano de incidencia (Ei = E\\) En este caso 
H sera normal al piano de incidencia, esto es n ■ Eli — n • H r = h ■ H t = 0 y 
procediendo analogamente llegamos a: 

r _ -Ej _ -771 cos 9 j + 772 cos 9 t 
^ Ey r/i cos 9 i + r/2 cos 9 t 


t_l = — 
E]_ 


T± = 

E\ 
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T = E ± = 2^2 COS Oj 

" E jj 771 cos 6i + 772 cos 0 t 

con el mismo problema de antes, de forma que es conveniente escribir: 


E \\ _ 

-771 cos 9i + 772 \j 1 - sin 2 9i 

Ei 

1 t \ 2 

II 

771 cos 9, + 772 Y 1 - ^ J sin 2 9i 


2 r /2 cos 9i 

Ei 

1 / \ 2 

II 

771 cos 9i + i 72 Y 1 - ( ^ J sin 9i 


y observamos que El y Et no estan en fase con El y que (T_l, T±) ^ (r||, Ty) 


Comentarios sobre los resultados obtenidos: A la vista de los re- 

sultados anteriores, podemos afirmar que una polarization lineal cualquiera 
(incidente) dara lugar a ondas rcflejadas y transmitidas con polarization 
eliptica salvo en el caso de incidencia normal (Oi = 0) en el cual no se 
distingue E±_ y E\\ y resultan las expresiones conocidas: 


ri = t_l = r(o) 


m - m 
V2 + m 


Ti = T ± = T{ 0) 


2 r/2 

m + vi 


Medios dielectricos 

Particularizamos los resultados obtenidos para el caso <r-\ = 02 = 0. Supo- 
nemos, ademas, que /j,-\ = 7*2 = /i 0 - Entonces, la segunda Ley de Snell se 
escribe: 


sin^t sin6*t n\ 



donde hemos introducido el valor 

n = \f^r 

conocido como indice de refraccion del medio junto con: 

m 

n 2 1 = — 

_ n2 _ 

indice de refraccion relativo de 1 con respecto a 2 
Observamos que: 

■ Si e 2 > e\ entonces 6 t < &i, esto es, el rayo se acerca a la normal a S 

■ Si ei > £2 hay dos posibilidades: 

• V0.j < arcsin(n2/7ri) sera 0 t > 6i 

• \/6i > arcsin(?i2/ni), 9t sera complejo. Excluimos, de momento 
esta posibilidad 
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Para valores de 9 t reales las formulas de Fresnel proporcionan coeficientes 
de reflexion y transmision tambien reales, por lo que las ondas transmitidas 
y reflcjadas estaran en fase (0°) 6 contrafase (180°) con respecto a la onda 
incidente. 


Angulo de Brewster: Veamos ahora si existe algun angulo para el cual 
no existe onda reflejada. Es facil comprobar que para la polarizacion perpen¬ 
dicular al piano de incidencia no es posible anular el fasor E\_. Sin embargo 
la ecuacion de Ty se anula cuando se verifique: 

—r] i cos 6i + rj 2 cos 9 t = 0 

que combinado con la primera la Ley de Snell conduce a un angulo de 
incidencia: 


9i = 9b = arcsin 


£2 


ei + e 2 


/ ^2 

= arctan . / — 


valor que se denomina angulo polarizante 6 de Brewster y que existe solo 
para la polarizacion paralela al piano de incidencia. 


Reflexion total: Volvemos ahora al caso en el cual obtenemos angulos 9t 
complejos. Esto ocurre a partir de un cierto angulo conocido como angulo 
critico dado por: 


a • 712 • l e 2 

9 C = arcsin — = arcsin ^ / — 

ni 


para el cual 9t = 7r/2, esto es, la onda transmitida sigue una direction 
paralela a la superficie S (y los pianos de fase constante ya no progresan 
hacia la direction z + ). 

El problema para angulos mayores al critico estriba en que los angulos 
obtenidos en el medio 2 pierden su significado fisico, si bien la solution ma- 
tematica sigue siendo valida. El problema esta en la suposicion de onda 
plana que postulamos al principio del apartado. Veamos esto con deteni- 
miento y estudiemos que tipo de soluciones obtenemos. 

Cuando 9t sea complejo tendremos: 

. Q 7i • a Oix+jPx 

sin 9 f = — sm 9; = -sm 9a 

72 a 2 +j/3 2 

y en el caso de medios sin perdidas, a± = a 2 = 0 por lo que sera: 


cos 9 t = — jn 2 \\J sin 2 9 t — nf 2 


(2.7) 


ecuacion en la cual la raiz lleva el signo que corresponde a la condition de 
que el campo nunca puede hacerse oo. 

Supongamos que n = x de forma que todos los puntos del medio 2 tienen 


valores positivos de x (ver figura 2.12l. Entonces: 


72h • r = jio^Je 2 [i 2 (x cos 9 t + zsm9 t ) = 

= juy/exg^-jxy/ sin 2 - n\ 2 + zsin 9 r ) 
y por tanto, el campo transmitido sera: 

E t = E to e~' l2fl - f = E to e~ a2X - jl3z 


( 2 . 8 ) 
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Figura 2.12: Reflexion total en una superficie coincidente con el piano xy 


donde E to se obtiene de la Formula de Fresnel correspondiente y: 


a 2 \/ sin 2 9 t - n\ 2 = — \ sin 2 0, 

v Vfl v 

U) 

(3 sin 0.j = -sin 0, ; 

v f i 


i ~ n 


12 


y observamos como el campo definido por ( 2.81 incluye dos terminos: 


Un factor e~ a2X que depende de x y que decae exponencialmente segun 
x —> oo lo cual confirma la election del signo negativo de la rafz de 
(2.71, y que define un conjunto de superficies equiamplitud en pianos 
paralelos a la discontinuidad (ver figura 2.13) 

Un factor e~^ z que depende de z y tiene aspecto de onda, a pesar de 
que su velocidad de fase: 


Vfl 

v„ = —-— 
sin 9i 

depende del angulo 9 t y que define un conjunto de superficies equifase 
perpendicular a las superficies equiamplitud y a la discontinuidad. 


S 


0 


Pianos 

equiamplitud 




Pianos 

equifase 


z 



t 


x 


Figura 2.13: Onda de superficie en medios sin perdidas para reflexion total 

Mediante las ecuaciones de Fresnel se puede llegar a las expresiones de los 
campos en los medios 1 y 2 y analizando el vector de Poynting en los dos 
medios se puede llegar a las siguientes conclusiones: 

A. El flujo de energfa de la onda reflejada coincide con el de la onda 
incidente 

B. El vector de campo E t no es nulo 


El origen del problems de los angulos 
complejos esta en que la suposicion de 
ondas planas (excepto para el caso de 
incidencia normal) es incompatible con 
la presencia del dioptrio 
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C. No existe flujo medio de energfa en el medio 2 (el de menor fndice de 
refraction) segun la componente normal a la discontinuidad, si bien 
existe una componente normal instantanea de flujo de energfa no nula, 
dado que los campos no son nulos 

D. Existe un flujo medio de energfa en el medio 2 en la direccion paralela 
a la discontinuidad que cae muy rapidamente cuando nos alejamos de 
ella. 

Por lo tanto, podemos decir que la excitation en el medio 2, en el caso 
de reflexion total, es una onda no transversal (no plana) confinada en la 
vecindad de la superficie de discontinuidad. Las soluciones de este tipo se 
conocen como ondas de superficie. Estas soluciones van a tener componentes 
transversales y longitudinales con respecto a la direccion de propagation (a 
lo largo de la discontinuidad). Por todo ello, podemos decir que: 

■ Los angulos complejos corresponden a ondas no planas, que presentan 
componentes longitudinales de campo 

■ Las componentes longitudinales estan relacionadas con la energfa re- 
activa almacenada en el medio 

■ Existe propagation por onda de superficie lo que da lugar a superficies 
equiamplitud y superficies equifase 

Los angulos complejos estan relaciona- 
dos con ondas no transversales (no pla¬ 
nas). 

Medios con perdidas 

En este caso, las leyes de Snell y Fresnel siguen siendo formalmente validas, 
pero nuevamente, el hecho de que 9 t sea compleja, conduce a una interpre¬ 
tation distinta. Supongamos que el medio 2 es conductor. De acuerdo con 
la Ley de Snell, tendremos: 


sin 9 t = (a — jb ) sin 9i 


de donde se sigue que: 


cos 9t = \J 1 — (a 2 — b 2 — 2jab) sin 2 0,; = pe 21 
por lo que el campo sera: 


Et — E to e 


-px—j(—qx-\-Oiiz sin Gi ) 


siendo: 


p = p{(3 2 cos 7 + a 2 sin 7) 
q = p(a 2 cos 7 — /? 2 sin 7) 


y ahora: 

■ Las superficies de amplitud constante son los pianos px = cte 

■ Las superficies de fase constante son los pianos — qx + aizs\n9i = cte. 

y observamos como las dos familias de pianos no coinciden y nuevamente la 
solution no es una onda plana (ver figura 2.141. 
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Pianos equiamplitud ... ' ' 

y equifase 



© 

s 

0 

V? 

vp : 

Z 

Pianos — 

equiamplitud. 

X 

\ 

\ 

- \ 


© 


Figura 2.14: Refraccion en un medio con perdidas 


Podemos definir una Ley de Snell modificada, que ahora depende del angulo 
de incidencia: 


n(9i) 


sin 9i 
sin d* 


1 

(X1 



sin 2 6i + q 2 


con lo cual la velocidad de fase sera: 


V2 (9i) 


U) 



sin 2 Oi + q 2 


Vf 1 
n(6i) 


Como caso particular, si el medio 2 es un buen conductor [a 2 >>), tendre- 
mos: 


Ot2 = P2 


UV2H2 


Se puede demostrar que en este caso, 
existe un angulo analogo al de Brewster 
para el cual el coeficiente de reflexion pa¬ 
ra la polarizacion de la onda paralela al 
piano de incidencia alcanza un mfnimo 
(no un nulo) 


que implica: 


0.2 




0 n{ 9 i) —> 00 Vf2 —► 0 


por lo que al ailmentar CT2 (6 disminuir /), los pianos de fase constante se 
alinean con los de amplitud constante y la propagacion en el seno del conduc¬ 
tor ocurre practicamente segun la direction a la normal a la discontinuidad 


(figura 2.15). 


Dielectrico 

© 

O 


. '!’ z 

Conductor 



© 


r 


Figura 2.15: Refraccion en medio buen conductor 


2. 21 Una onda plana con una frecuencia de /=3 GHz se propaga en un 
medio material infinito con e r = 7 y /x r = 3. Calcular A, Vf y 77. 

Sol. A = 0.0218m; v p = 6.55 • 10 7 m/s; 77 = 246.8H 


2.8 EJERCICIOS 
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2 . 22 Calcular la profundidad de penetration, S s , en el Aluminio (a = 
3.86 • 10 7 S/m) a / = 10 GHz. 

Sol : 6 s (Al) = 8.14 • 10~ 7 m 


2. 23 


Calcular la constante de atenuacion en tierra seca (e 7 = 4, a = 10 


-3 


S/m) y en mar (e 1 = 80, a = 4 S/m). 


/ 

a (dB/m) Tierra 

a (dB/m) Mar 

15 KHz 

0.06 

4.22 

150 KHz 

0.20 

13.37 

1.5 MHz 

0.56 

42.24 

15 MHz 

0.80 

132.6 

150 MHz 

16.83 

389.1 !!! 


2. 24 Determinar la polarization de una onda plana en los siguientes casos: 


A. E = E a (x + y)e~ jk ° z 

B. E = E a (x + jy)e- jk ° z 

Sol : A) Lineal; B) Circular a izquierdas 


2. 25 


Describir el estado de polarization de las siguientes ondas: 


A. E = E 0 cos (kz — ut)x — E 0 sin(fc2 — u>t)y 

B. E = E 0 cos (kz — ut — n/2)x — E 0 sin(kz — u>t)y 

C. E = E 0 cos (kz — ut)x — E 0 cos (kz — u>t)y 

y determinar sus componentes magneticas y la energia media transportada 
en un periodo por unidad de superficie. 

Sol: 

A. Circular a derechas. H = j!^{cos(u}t — kz)y — sen{u>t — kz)x\ ; < 

S>=T2W|£o| 2 * 

B. Lineal con direction: u= = j^§^[x + y]cos(cot —kz+ |); 

<S>=T^\Eo\ 2 Z 

C. Lineal con direction: u = x — ~^y- H = j§§^[x + y]cos(u>t — kz); 
<S>=j^\E 0 \ 2 z 


2 . 26 


La amplitud compleja del campo electrico de una onda plana mo- 


nocromatica que se propaga en el vacio esta dada por: 


E = [(V 3 + j)x + 2 jy\e 12 °^ 


Determinar la frecuencia de la onda, el vector unitario en la direction de 
propagation, la polarization, la expresion del campo magnetico, la razon 
axial y la densidad de potencia media. 
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Sol : /=1 GHz; k = z; Elfptica a izquierdas; H = 2jx + (\/3 + 

j)y)e - - 720 ”'/ 3z A/m; ra = \/3; ^z W/m 2 


2. 27 


Dado el campo electrico: 


—; r— A 207T . 207 r 

A = V3a;cos(a;f-— 2 je 3 

O 


(x + 2y) sin(w< 


207 T , 20tt .. 


que se propaga en un medio conductor, clescriba su polarizacion y calcule 
la constante de atenuacion en dB/m. 

Sol : Elfptica a izquierdas; a = 57,86ndB/m 


2. 28 


El campo magnetico en una cierta region del espacio responde a la 
siguiente expresion: 

H = H 0 (x + jy)e~i kz 


Calcular el campo electrico instantaneo, la curva descrita por el vector cam¬ 
po electrico en un punto fijo del espacio y la densidad de potencia media e 
instantanea. 

Sol : E = —1207ri7o [cos(ujt — kz)y + sen(u>t — kz)x]\ Polarizacion circular a 
izquierdas; < S >= QQ-rrHlz 


2. 29 


Una onda plana monocromatica se propaga en el vacfo segun z 
y esta constituida por la superposicion de dos componentes linealmente 
polarizadas cuyos valores instantaneos en z = 0 estan dados por: 


Ei(0,t) =3 cos cotx 
E2(0,t) =2cos(wt + 7 r/ 2 )y 


Determinar la relacion entre los valores maximo y mfnimo del campo total, 
la polarizacion y el vector de Poynting medio asociado a la onda y a cada 
una de sus componentes por separado. 

Sol : RA = 3/2: Elfptica a izquierdas; P av = W/m 2 Pi = 

W/m 2 P 2 = 2 iM 


2. 30 Determinar la amplitud compleja del campo electrico de una onda 
plana monocromatica con polarizacion elfptica a izquierdas sabiendo que 
se propaga en el vacfo segun una direction que se considera como eje z\ la 
relacion entre el valor maximo y el valor mfnimo del campo electrico es N; 
la direction segun la cual se mide el maximo valor del campo electrico se 
toma como eje y\ el valor medio de la densidad de potencia que transmite 
la onda es P(W/m 2 ) y el campo electrico en el piano z = 6 tiene direction 
x en el instante t = 0. 




2. 31 Encontrar el vector de polarizacion de una onda plana si e = 30° y 


r = 135 °. 
Sol: e = 


2. 32 


Calcular el vector de polarizacion de un campo con polarizacion 
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eli'ptica a derechas cuya razon axial es 3 dB y su semieje mayor corresponde 
al eje x. 

Sol: e = 0.816x — j0.577y 


2. 33 | Calcula el coeficiente de reflexion de una onda plana sobre un piano 
conductor infinito de Cu (u = 5.813 • 10 7 S/m) para /= 1 GHz. 

Sol: T « 1 Z179.9°; T « 6.181 • 1CT 5 Z45° 


2. 34 


Una onda plana que se propaga en un medio de e r = 4, /i r = 1 y cu- 


yo fasor de campo electrico es E l = 2 • 10 ~ 3 e~ : ’ kz yV/m incide normalmente 
sobre un medio igual al vaefo. Determinar el campo magnetico incidente, 
el coeficiente de reflexion y la ROE, los campos reflcjados y transmitidos 
(electricos y magneticos) y las densidades de potencia media incidente, re- 
flejada y transmitida. 


Sol: 


A. Hi = -MM e -i k *x{A/m) = (A/m) 

B. Ti(^) = §e? 4fc ° z ROE=2 

C. E r = |0.002 e j2koZ y{V/m)] # = §0.002 e-^ k ° z y{V/m); 

H r = ei 2k ° z x(A/m); H l = ^§e~3 k ° z x(A/m) 


D. 


PI 


10~ 6 

30tt 


z(W/m 2 ) Pl v 


io ~ 6 

270tt 


z(W/m 2 ) P* v 


8 - 10~ 6 

270tt 


z(W/m 2 ) 


2. 35 | Una onda plana que se propaga en el espacio libre incide normalmen¬ 
te sobre un material dielectrico sin perdidas. La reflexion producida origina 
en el espacio libre una onda estacionaria, y por medio de mediciones se de- 
termina que: la relacion de ondas estacionaria (ROE) es 2; la distancia entre 
dos maximos consecutivos de campo electrico es 6 cm y el primer maximo 
del campo electrico esta situado a 3 cm del dielectrico. Calcular: 


A. La frecuencia 

B. La constante dielectrica relativa de la lamina 

C. El porcentaje de energfa incidente que atraviesa la lamina 

Sol: f= 2.5 GHz; e r = 4; 88.9% 


2. 36 | Una onda plana de frecuencia 3 GHz se propaga en el vaefo e incide 
normalmente sobre una superficie dielectrica de permitividad 16eo y espesor 
5/8cm, a la que sigue otra de permitividad eo y 5cm de espesor, tras la que 
se halla una nueva lamina de 4eo y 5/4cm de grosor a continuacion de la 
cual se tiene de nuevo el vaefo. 


A. ^Cual es el porcentaje de la potencia incidente en el primer medio que 
se transmite al ultimo? 

B. zQue medio o medios podrfan eliminarse de la estructura dada man- 
teniendose el mismo porcentaje de potencia transmitida? 
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Sol : 64%; 2 y 3 


2. 37 Una onda plana incide oblicuamente desde el vaci'o sobre un conduc¬ 


tor perfecto como se indica en la figura. La amplitud compleja del campo 
incidente es: 

E = 5xe-^^ v+z) 


Y 


(1) 

(2) 

9- 

A 

tx 



A. Calculese la expresion de los campos totales en funcion del tiempo. 

B. ^Existe algun punto en el cual se anule el campo electrico para cual- 
quier valor del tiempo? 

Sol : E(t ) = 10sen(\/27r^)sen(a;t — s/2-k y)x; Los pianos z = 
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CAPITULO 3 


Ondas guiadas 


Planteamos el problema de la transmision de energfa electromagnetica entre 3.1 INTRODUCCION 
dos puntos cumpliendo dos exigencias: 

■ Debe llegar la totalidad o la mayor parte de la energfa enviada desde 
el primer medio 

■ El sistema debe permitir el paso de las potencias que sea necesario 
transmitir en la practica 

En el caso de las bajas frecuencias, el problema se resuelve empleando dos 
conductores paralelos separados entre si una distancia pequena comparada 
con A. Por este motivo, no se producen perdidas por radiacion. 

Si se quiere mantener esta situation al elevar la frecuencia, es necesario 
disminuir dicha distancia, lo que limita la potencia maxima transmisible por 
la formation de arcos electricos. Todo ello, limita la transmision por lfneas 
bifilares hasta frecuencias del orden de los MHz, lo que obliga a recurrir a 
otros artificios para suprimir la radiacion de los conductores paralelos. 

Una primera solution es el cable coaxial en el cual uno de los conductores 
engloba completamente al otro. Para soportar el conductor central se puede 
utilizar un medio dielectrico que rellene el espacio entre los conductores o 
simplemente rodajas de dicho material aislante. 

Generalizando lo anterior, introducimos unos elementos nuevos, las gufas 
de ondas 6 guiaondas que, en general, consisten en un tubo cilfndrico de 
paredes metalicas y lleno de aislante, aire o el vacfo en las cuales el soporte 
de la energfa que se transmite es el campo electromagnetico que, para una 
determinada frecuencia, se propagara a lo largo del tubo si se han elegido 
adecuadamente las dimensiones de la seccion recta de este ultimo. 

La propagation estara tanto mas confinada dentro de la gufa cuanto mas alta 
sea la frecuencia de trabajo porque se reducira la profundidad de penetration 
del campo e.m. en el conductor metalico que lo constituye. 

Las gufas se conocen segun la forma de la pared que constituye su seccion. Si 
en el interior de la gufa no hay mas que un aislante homogeneo e isotropo, 
se dice que la primera es homogenea. En caso contrario, se especifica la 
geometrfa con un calificativo especial como p. ej. guia parcialmente rellena 
de dielectrico. 
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3.2 GUIAS CILINDRICAS 
HOMOGENEAS. 
PLANTEAMIENTO TEORICO 
DE LA PROPAGACION 


En coordenadas cartesianas £1 = x, 
£ 2 = y; en coordenadas cilfndricas £1 = 

p, & = <f>; 


Existen otros mecanismos para conseguir la propagacion guiada basados 
en la propagacion de ondas superficiales sobre un conductor metalico. Un 
ejemplo son las lineas microstrip, que consisten en un piano conductor infe¬ 
rior, una capa de sustrato dielectrico y una tira de conductor situada sobre 
este, estructura que permite la propagacion guiada si las dimensiones son 
las adecuadas. 

En primera aproximacion hacemos la abstraction matematica de considerar 
que el conductor que forma la guia es perfecto. Como se ha mendionado 
anteriormente, esto es tanto mas cierto cuanto mas alta sea la frecuencia ya 
que la profundidad de penetration sera menor. En consecuencia, no existe 
campo e.m. en el metal que forma la guia y el problema se reduce a hallar 
la propagacion en un medio homogeneo isotropo aislante (o conductor im- 
perfecto) y que llena la region del espacio interior a una superficie cilindrica 
donde esta en contacto con otro medio, conductor perfecto. 

Para simplificar el calculo, suponemos que el dielectrico no tiene perdidas, 
por lo que se trata de un aislante perfecto. 

Dada la simetria de traslacion, caracterizamos cualquier punto P del espacio 
mediante tres coordenadas (^1,^2, z) de forma que: 

■ Los puntos situados en una recta paralela a las generatrices del cilindro 
en cuestion tienen las mismas coordenadas (^1,^2)- 

■ Todos los puntos situados en un piano perpendicular a dichas genera¬ 
trices tienen la misma z. 

■ La diferencia entre las coordenadas z de dos puntos es igual a la dis- 
tancia geometrica que separa los pianos que pasan por tales puntos y 
que son perpendiculares a las generatrices del cilindro. 

Se trata entonces de resolver las ecuaciones de ondas: 

A E - 7 2 0 E = 0 
AH - 7 2 0 H = 0 

siendo 7^ = — w 2 /ie = — k 2 anadiendo las condiciones de contorno que impo- 
ne la guia a los campos E y H. Junto con el campo electromagnetico en el 
interior de la estructura, podriamos calcular la densidad lineal de corriente 
segun: 

J s = n x H 

Podemos escribir el campo E como (para el campo H el desarrollo es analo- 
go): 

E — E t + E z 

donde E t es la componente tangencial y E z la componente longitudinal, con 
respecto al eje de la guia, por lo que resultara: 

A E t + A E z - 7 2 Q E t - 7 2 0 E z = 0 

En sistemas que presenten simetria de traslacion, A E t es un vector trans¬ 
versal a z mientras que A E z tiene la direccion de z, por lo que: 

A E t - 7 2 Q E t = 0 

A E z - 7 2 q E z = 0 
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Aclemas se verificara: 


A E z — (A Ez)z — ( A t E z + 


d 2 E z 

dz 2 


siendo A t la parte de A que no contiene derivadas con respecto a z, luego: 


A E z — A Et + 


d 2 E z 

dz 2 


- 7 %E t 


ecuacion a la que se le puede aplicar la tecnica de separation de variables 
para escribir: 

E z =F e {^,^)A e {z) 

obteniendose: 

A t F e til' 


F e 


- 7 2 = 0 

A P 7o 


lo que supone que: 


K 

A P 


= 7 


donde q 2 es una constante de separacion y la solution para E z es de la 
forma: 

E z = F e (^,^)(Ae-^ + Be^) 

Siendo F e (£ 1,^2) la solucion de la ecuacion: 

A t F e - JcF e = 0 


con: 


222 
7 C = lo ~ 7 


Si se elige una de las dos exponenciales, tendremos: 

E z = F e {^^ 2 )e-^z 
y operando de la misma manera para H z : 

H z = Fhitu&e-^z 


donde, como hemos dicho anteriormente, F e y Fh son dos soluciones inde- 
pendientes de la ecuacion: 

A t F - 7 2 C F = 0 

cuyas constantes de separacion pueden ser distitintas y clependeran de las 
condiciones de contorno que apliquemos a cada campo. 

Atendiendo a las expresiones de los campos, podremos clasificar las solucio¬ 
nes en: 

■ Modos Transversales electromagneticos 6 TEM: E z = 0 y H z = 0 

■ Modos Transversales Magneticos 6 TM: H z = 0 

■ Modos Transversales Electricos 6 TE: E z = 0 

quedando justificados los calificativos de transversal magnetico, electro- 
magnetico y electrico dados a los modos en cuestion puesto que se anulan 
las componentes longitudinales del campo o los campos correspondientes y 
estos son vectores transversales. 

Pasamos a estudiar cada uno de los casos introduciendo una impedancia 
que relacione los campos transversales E y H de cada modo de la misma 
forma que hicimos con las ondas planas. 
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Modos TEM 


Este caso corresponde a la solucion trivial de las ecuaciones 
Los campos unicamente tienen componentes transversales, 
ecuacion a resolver es: 




cuya solucion se puede expresar como: 


para F e y F h . 
por lo que la 


E = F e (£ i,6)e- 


siendo la funcion F e la solucion de: 

V t x F e = 0 


junto con las condiciones de contorno correspondientes. Observamos como 
F e nos proporciona la variation transversal del campo en cualquier section 
recta de la gufa. 

Se demuestra, ademas, que se verifica: 

H = —z x E 

Vo 


Ztem 


Z x E t 
H t 


= Vo 


que coincide con la impedancia intrinseca del medio material para una onda 
plana. 


Modos TM 

En este caso, F[ esta contenido en pianos perpendiculares a z. Los campos 
son: 

E z = F e (H 1 ,&)e^ z z 

clonde F e cumple la ecuacion: 


A t F e - Y c F e = 0; 


2 2 2 
7 C = lo ~ 7 


H t = - J ^t.E z x 2 

E t = ^X7 t E z 
7c 


Ztm = 


z x E t _ 
H t 


Modos TE 

En este caso, E esta contenido en pianos perpendiculares a z. Los campos 
son: 

H z = Fbitu&e-^z 
donde Fh cumple la ecuacion: 

A t F h - 7 c F h = °i 7c = 7o - 7 2 
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E t = - 3 -^V t H z x £ 


H t = -EV t H z 
7c 


Zte 


Z x E t 
H t 


7 


Notar que como E z y iL, son independientes, las soluciones con E z ^ 0 
y iL z 0 tambien lo seran y podran obtenerse como superposition de 
las anteriores, esto es, no se trata de una familia separada de soluciones 
distitintas de las que ya se han descrito. 


Comentarios sobre las soluciones 


La solution general del problema es combination lineal de las soluciones 
TEM, TE y TM, siendo las condiciones de contorno las que determinan las 
constantes de separation 7 y los coeficientes de la combination correspon- 
dientes a cada caso. De esta forma: 

■ Las condiciones de contorno laterales definen la forma especffica de 
variation de los modos con las coordenadas £ 1 , £1 y z, esto es, las F e , 
Fh y 7, cuya determination es un problema bidimensional en (£ 1 , 62 ) 


Notar que la tecnica de separacion de la 
variable z no funciona si la frontera de 
la gufa no es un cilindro recto o bien las 
condiciones de contorno dependen de z, 
ya que en ese caso, F(£ 1,^2) dependerfa 
tambien de z 


■ Las condiciones en los pianos z = cte. correspondientes al generador 
y la carga determinan cuantos y cuales son los modo que se propagan 

Para cada geometrfa de guia, existe un sistema de coordenadas que facilita 
su solution, y que viene determinado por las simetrias existentes en las 
condiciones de contorno. 

Asi, para el caso de una guia con section rectangular, es apropiado emplear 
las coordenadas cartesianas. Para el caso de la guia cilindrica, es natural em¬ 
plear las coordenadas cilindricas, que pasamos a definir, suponiendo ademas 
la condition de que el sistema obtenido sea ortogonal, por ser los unicos que 
permiten hallar soluciones sencillas. 


Las constantes 7 C son tambien incogni¬ 
tas del problema. Desde un pun to de vis¬ 
ta matematico, las ecuaciones de los F e 
y Fh plantean el problema de la deter- 
minacion de las funciones propias F y 
de los valores propios 7 C del operador de 
Laplace bidimensional. 


3.2.1 _ 

Analisis de la variation con z 

La constante de propagacion que hemos obtenido viene dada por: 

7 = \hl-ii 

dependiendo 7 C de las condiciones de contorno que gobiernan F e y Ff, . De la 
misma manera que hicimos con las ondas planas, denominamos a 7 constante 
de propagacion, que dependera tambien de la frecuencia y que normalmente 
expresaremos como: 

7 (w) = a(u) + j (3 (u) 

siendo a la constante de atenuacion y /3 la constante de fuse. De la misma 
forma, podemos definir: 

■ A(w)= ongitud de onda 

■ Vf( cu)= velocidad de fase 
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■ Vg(uj) = d/3/du= velocidad de grupo 

Lo novedoso con respecto al planteamiento de las ondas planas es que solo 
el conocimiento de las condiciones de contorno que corresponden al modo 
cuya solution calculamos fija el valor de la constante de propagation 7 . 

3.2.2 _ 

Condiciones de contorno laterales 

Las condiciones de contorno en z=cte. determinan los modos que es nece- 
sario considerar para obtener la solution completa del problema. 

Las condiciones de contorno laterales en (£ 1 , £ 2 ) determinan las caracteristi- 
cas de cada modo (a traves de 7 C ) y la forma de los campos. 

Las ecuaciones de F e y F ^ solo tienen solucion sencilla cuando se pueda 
aplicar la tecnica de separation de variables y esto se puede hacer cuando 
las condiciones de contornop se den sobre pianos de = cte. 

Los tipos de condition de contorno mas habituates son: 

A. Conductor perfecto (cr —► 00 ) 

B. Conductor no perfecto (buen conductor) 

C. Problemas abiertos (Ej. la Knea bifilar). Se supone un conductor en el 
infinite que cierra el problema 

D. Cualquiera de los anteriores cuando existan discontinuidades de medio 
en superficies distintas de z=cte. como p. ej. la Knea microstrip. 

y unicamente el primero admite una solucion anali'tica sencilla. 

3.2.3 _ 

Condiciones de contorno de conductor perfecto 

Supongamos una situation como la de la figura|3TT| En este caso, el conduc¬ 


es —> CO 



Figura 3.1: Condicion de contorno de conductor perfecto 


tor perfecto impone: 


n x E t = 0 
E z = 0 
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Modos TE 


Aplicando las condiciones anteriores se llega a: 


dF h 

dn 


c 


= 0 


Modos TM 


En este caso directamente: 

F e \ c = 0 


Modos TEM 

En este caso tendremos: 

h x E t = 0 =7 n x FEt = 0 
Como ademas se verifica: 

Vt X pEt — 0 

podremos decir que: 

pEt = — V t 0 

siendo </> = £ 2 ) una funcion escalar que verifica la ecuacion de Laplace 

bidimensional: 

A t </> = 0 

por lo que el problema coincide con el planteamiento de un problema elec- 
trostatico bidimensional, esto es, una vez encontrado <j), el campo electrico 
vendra dado por: 

E t = -V t 0 • 

Por ello, y dado que Pfa es normal a C, esta curva debe ser una linea de <f> 
constante, asi que en los conductores podremos imponer <j >| c = cte. 

Una primera conclusion importante es la siguiente. Si la frontera del con¬ 
ductor es simplemente conexa la unica solucion posible es <j> = cte. en toda 
la region y por tanto FEt = 0. Por tanto, en una region limitada por un 
conductor perfecto simplemente conexo no puede existir modo TEM. 


3.2.4 _ 

Caracterfsticas de los modos 


Se puede demostrar que cuando las condiciones de contorno son las de un 
conductor perfecto, se cumple que 

7c <0 


Por tanto, la expresion q 2 = 7 ^ — q 2 sera positiva o negativa dependiendo 
de la frecuencia, siempre y cuando el medio material no presente perdidas. 

El valor de la pulsation para el cual la expresion anterior cambia de signo 
se conoce como pidsacion (frecuencia) de corte y viene dado por: 


w c 


1 

VJF 


\[-ii 


y las dos posibilidades son: 
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Si f < f c las impedancias de modo son 
imaginarias puras por lo que el flujo me¬ 
dio de potencia es nulo 


■ Si u> < u c (/ < fc) => 7 2 > 0 y por tanto 7 = a por lo que el 
modo no se propaga, esto es esta al corte. 

■ Si to > uj c (/ > f c ) => 7 2 < 0 y por tanto 7 = j/3 por lo que el 
modo se propagara. 


En funcion de la frecuencia de corte, podemos escribir: 


7 = ±7o 




Zte = ± 



Z T m = ±7 



ft 

P 


correspondiendo el signo ’+’ al caso / > f c y el a / < f c 


Podemos representar esta situacion en un diagrama de Brillouin (figura 3.21 
donde las curvas de a y j3 vienen definidas por las ecuaciones: 



Observamos que: 



( v ^) 2 



( v /= 7!) 2 


■ El modo TEM no posee frecuencia de corte 


El diagrama representado corresponde a un unico modo. El diagrama 
completo debe incluir todos los posibles modos que se pueden excitar 

(fig.pl 


A partir de este diagrama podemos calcular la velocidad de fase y de grupo 
que corresponde a cada modo en cualquier frecuencia. Notar que la velocidad 
de fase es siempre superior a la velocidad de grupo (salvo en el caso TEM) 
motivo por el que se conoce a estas ondas como ondas rapidas. 

Algunos comentarios sobre las soluciones encontradas: 
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Figura 3.3: Diagrama de Brillouin para todos los modos 


■ Si la frecuencia de trabajo es inferior a la frecuencia de corte del modo 
con menor — 7 ^ ningiin modo se propagara. 

■ El modo con menor f c se denomina modo dominante y el resto modos 
superiores. 

m En los sistemas capaces de soportar modos TEM, no existe frecuencia 
de corte absoluta. 


■ Se puede demostrar que el modo de frecuencia de corte menor es 
siempre un TE. 

■ Se puede demostrar que en un sistema capaz de soportar modos TEM, 
el primer modo superior es siempre un TE. 


El ancho de banda en el cual tenemos un linico modo propagandose 
se denomina ancho de banda monomodo B (ver figura 3.41 



P 


Figura 3.4: Ancho de banda monomodo B 


3.2.5 _ 

Ortogonalidad de los modos 

Dadas dos soluciones de las ecuaciones de Maxwell para una gufa de ondas, 
en que los campos electrico y magnetico sean E, H y E', H' respectiva- 
mente, definimos como producto escalar de ambas soluciones los numeros 
complejos: 
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/ H ■ H' da 

Jt. 

siendo E una section recta cualquiera de la gufa. 

Se puede demostrar que estos productos escalares se anulan para dos cua- 
lesquiera de las soluciones particulares estudiadas, esto es, los modos TM, 
TEM y TE. Para ello se debe realizar un estudio de las autofunciones del 
operador A. Empleando el lenguaje habitual de calculo vectorial, diremos 
en consecuencia que los modos distintos son ortogonales. 

Ademas, el conjunto de los modos obtenidos es completo. Esto significa que 
cualquier solution de campos definida en un piano 2 = cte de la guia se 
puede escribir como una combination lineal de soluciones TE, TM y TEM: 

E = E TM,TEMfniEi + Y,TEnjEj 


3.2.6 _ 

Flujo medio de potencia a traves de una section recta de la guia 


El flujo medio de potencia sera igual a la parte real del flujo del vector de 
Poynting complejo: 

S=^ExH* 

a traves de tal section recta. Designandola como E y tomando en ella como 
sentido positivo de su normal el del vector unitario z, el flujo de S valclra: 


S ■ da = / (S ■ z)da = - 
Jy. 2 , 


[E x H*]da 


desarrollando los campos en sus componentes tangenciales y longitudinales, 
se llega a: 


I S-dff= -J [Et xH* t ]da 


Ademas hay que destacar que: 


■ En el vector de Poynting solo influyen las componentes transversales 
de los campos 

■ Solo existe propagation de energfa cuando el modo no esta al corte. 
En caso contrario, existira energfa reactiva (energfa em. almacenada 
en la gufa). 

■ Cuando varios modos se propagan simultaneamente, no existe trasvase 
de energfa entre ellos. Esto se debe a las propiedades de ortogonalidad. 


3.2.7 _ 

Estudio de la gufa rectangular 

El convenio habitual es llamar a al lado Se trata cle un sistema invariante con z cuya seccion es: 
largo, que se coloca en el eje x . 

Dado que las condiciones de contorno se dan en pianos x = cte y = cte es 
apropiado utilizar las coordenadas cartesianas. Entonces, en la resolucion 
de las ecuaciones de F e y Fh- 

d 2 F d 2 F 


dx 2 dy 2 


llF = 0 
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Figura 3.5: Gui'a rectangular 


podremos aplicarle la tecnica de separacion de variables: 

F(x.y) = X(x)Y(y) 
que sustituida en la ecuacion da: 

X"Y + XY" - 7 2 c XY = 0 

que descomponemos en: 

x - 


Y" 

— - -k ' 2 
Y ~ v 


junto con: 


,V * = -h Z _ k 2 
Ic 

siendo k x y k y constantes de separacion. Las soluciones de esta ecuacion 
son: 


X(x) = A sin k x x + B cos k x x 
Y(y) = C sin k v y + D cos k v y 


Ahora aplicamos las condiciones de contorno obtenidas en el apartado |3.2.3| 
para cada modo. 

Modos TE en una guia rectangular 

La condicion es: 


dF h 

dn 


C 


= 0 


que en este caso se traduce en: 


dXY 

dx 

dXY 

dy 


y=0,b 


x=0,a 


= o 

= o 














88 ONDAS GUIADAS 


De la primera de ellas se extrae que: 

C = 0 

T17T 

sin kyb = 0 => k y — — 

y de la segunda: 

A = 0 

sin k x a = 0 =4> k x = - 

a 

por lo que la solution para el campo magnetico longitudinal sera: 
ft irnr mr 

Hz(mn) COS X • COS y • 6 Z 

K ’ a b 

siendo la constante de propagacion: 

7 

Notar que: 

■ Para cada pareja de valores m y n existe un modo TE mn solution. 

■ La solution con n = in = 0 no esta permitida 




donde: 



Notar que la relation entre las componentes transversales de E y H es la 
impedancia de modo. 

Modos TM en una guia rectangular 

La condition es: 

F e \ C = 0 

que en este caso se traduce en: 

XYx=0,cr,y=0,b = 0 
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De donde podemos extraer la solution 

f! „ . ITITT . U7T 

E~(mn) = Q sin-x • sm —ye lmn z 

y a b 

siendo la constante de propagation: 


Notar que: 



■ Para cada pareja de valores my n existe un modo TM mn solution. 

■ Los indices my n deben ser ambos distintos de 0. 


La solution completa es entonces: 


E z =P mn sin(“W^V— 


a 


) V b ) 

( nrvKX \ . 
cos (-) sm 

(- ° 




p _-7mnTO7T 

2 -^Enn 

„ -Jmnmr . / mwx\ (mry\ 

Ey =-- p rnn Sill ( - ) COS ( ) e Tm " 2 

y bkl 


a / V b ) 




donde: 



Notar que la relation entre las componentes transversales de E y H es la 
impedancia de modo. 

La tabla siguiente representa los modos TE y TM con frecuencias de corte 
mas bajas para el caso a = 2b. Notar que existen modos degenerados, esto 
es, con la misma frecuencia de corte. 



TE 10 

TEqi, TE 20 

TEn, TMn 

TE 21 , TM 2 1 

f c /fcT Eiq 

1 

2 

2.24 

2.82 


Cuadro 3.1: Frecuancias de corte para a = 2b 

En general, si a > 6, tal y como venimos haciendo, el modo fundamental es 
el TEio para el cual 

f I = 1 — 

Jc ' TE10 y/JIe2a 

o, en funcion de las longitudes de onda de corte: 

/ > fc => A < A c 


siendo A c : 


C 


2 


A, 


i 


b 2 
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que para el modo fundamental resulta: 

^ c \te 10 = ^ a 

Esta ecuacion podemos interpretarla en terminos de semilongitudes de onda 
que caben en a. Asi, para saber si el modo fundamental se propaga, hay que 
comprobar si una semilongitud de onda es inferior a a. Esto es aplicable a 
otros modos que tengan un mdice 0. 

La expresion completa de los campos para el modo fundamental TEio sera: 

A partir de la expresion de H pode¬ 
mos calcular las densidades lineales de 
corriente J s aplicando la condicion de 
contorno que los liga 

H z 


H x 


siendo: 

7 = 



=P cos e~ lz 

=^ a Psin( 7 rX )e-- z 
it \ a / 

=0 

7 « . / 7TX \ 

=—Psm — e ' 

7r V a / 

=0 


3.2.8 _ 

Estudio de la gui'a circular 

En este caso es conveniente aplicar el sistema cilfndrico de coordenadas. La 


|.i,e 



Figura 3.6: Gufa circular 


ecuacion a resolver es 

1 d f dF\ 1 d 2 F 2 _ n 

r dr \ dr) r 2 dcj) 2 c 

y podremos aplicarle la tecnica de separation de variables segun: 

F(r,(j>) = R(r)H(</>) 

que sustituida conduce a 
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i?" + -R' + 
r 




R = 0 


siendo k$ una constante de separation. Las soluciones de estas ecuaciones 
son: 


H((j>) = A sin k^(f> + B cos k^(j> 

R(r) = CJfc^v/^Tfr) + DN k<j> (t^r) 


Las funciones de Bessel desempenan un 
papel equivalente al de los senos y cose- 
nos del caso rectangular 

Modos TM en una guia circular 

Debe ser: 

= 0 => 

siendo p nm la raiz n-esima de J n = 0. El campo longitudinal resultante es: 

E z (nm) = PJn (^^ r ) C0S ’ e ~ ^ 

y la frecuencia de corte esta dada por: 

/.i _ 1 Pnm 

J C\ r p fi] c\ /- 

nrn 27 T^/fle a 


y dado que la solucion es una funcion periodica en <f>, k$ debe ser entero. J y 
N las funciones de Bessel de orden n de primera y segunda especie. Por otra 
parte, N n presenta una impropiedad en el origen, asi que debe ser D — 0. 


Ahora aplicamos las condiciones de contorno obtenidas en el apartado 3.2.3 
para cada tipo de modos. 


Modos TE en una guia circular 

Procediendo de forma analoga: 

H z (nm) = QJn C0S n( t> ’ e ^ ^ ^ 

siendo p' nm la rafz n-esima de J ' n =0 (el valor to = 0 esta prohibido). y la 
frecuencia de corte esta dada por: 

1 p’ 

r I _ Jrnm 

J C I r p ]\/[ / 

nm 27 Ty/JIe a 

Un modo particularmente es el TEoi, 
primer modo circular, para el cual los 
campos no dependen de (f). Esto tiene 
aplicacion en la construction de juntas 
rotatorias de antenas de exploracion. 


siendo el modo dominante el TEn (p' n = 1.84). 

3.2.9 _ 

Soluciones para medios dielectricos con perdidas 


Cuando trabajamos con un dielectrico con perdidas, la constante de propa¬ 
gacion va a ser compleja a cualquier frecuencia. En efeto: 

7 2 = 7o - 7c 2 = - Je") 7c) = (« + i/3) 2 


Esto significa que ya no existen modos completamente al corte. Ademas 
las impedancias de modo van a tener parte real e imaginaria a todas las 
frecuencias. En este caso se define la frecuencia de corte como aquella para 
la cual a = /3 (ver figura 3.7) Se demuestra que la frecuencia de corte 
obtenida de esta manera coincide con la que resultarfa de analizar la guia 
rellena de un dielectrico con el mismo e! y sin perdidas. 
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Figura 3.7: Gui'a rellena con dielectrico con perdidas 


Soluciones para medios dielectricos con bajas perdidas 

La constante de propagation puede escribirse como: 

7 = shl ~ 7c = v / -w 2 /i 0 e 0 e r ( 1 - j tan 5) - 7 2 (3.1) 

Si tan (5 <<, podemos emplear el desarrollo en serie de Taylor: 

\/a 2 + x 1 ~ a + — 

2 a 


de forma que (O resulta: 

7 = \/-fc 2 (l - jtam5) - 7 2 = \/-fc 2 — 7 2 + jfc 2 tan <5 
_ - 2 , jfcotan,5 tan S 

siendo (3 la constante de fase correspondiente al caso sin perdidas. Por lo 
tanto podremos escribir: 


Old 


kg tan S 

2/3 




donde = uj^/pe r e 0 


k 0 tan 8 

2 


(T£/, TM ) 
(TEM) 


3.2.10 _ 

Condiciones de contorno para conductor con perdidas 

Cuando la guia este delimitada por un medio conductor real, existiran cam- 
pos en el seno de ese conductor (o al menos en una distancia tanto menor 
cuanto mejor conductor sea). En cualquier caso, los campos estaran confi- 
nados en una region pequena que dependera de 8 (profundidad de penetra- 
cion). Esta situation se indica en la siguiente figura para una guia de section 
rectangular. 

En estas circunstancias: 

■ No es posible encontrar soluciones que cumplan las condiciones de 
contorno salvo para modos hibridos (con E z y H z simultaneamente) 


















3.3. Los modos TEM 
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Figura 3.8: Gufa rectangular con perdidas en el conductor 


■ Las condiciones de ortogonalidad dejan de ser validas. 

■ Existe propagation a todas las frecuencias. Ademas, las frecuencias de 
corte ya no son las mismas. 

Dado que la resolution de este problema es muy compleja, se suele hacer 
una aproximacion de medio buen conductor para que sea tratable, de forma 
que tomamos: 


(X — (Xg + a c 
P = Pg 

siendo a g y P g las que corresponden a una gufa con conductores perfectos y 
a c una constante de atenuacion relacionada con las perdidas por efecto Joule 
en los conductores que se calcula en cada caso en funcion de la geometrfa. Se 
puede demostrar, que para un buen conductor, una aproximacion aceptable 
viene dada por: 


a W c _ \§ c ^g(O) • H(0)*dl 

2W T Re ff s E t ( 0) x H t (0)* ■ ds 

que depende del modo que estemos analizando. Asf, para el modo TEio de 
una guia rectangular, esta expresion se convierte en: 


Oic 


R s 

aPbPgkrj 


(2bTT 2 + a 3 k 2 ) 


y para el modo TEn de una gufa circular: 


R, 


aPgki'j 


n' A — 1 

P n l 


siendo 

Rs = 

la resistencia superficial del conductor. 

Los modos TEM tienen dos propiedades interesantes: 3.3 LOS MODOS TEM 

A. Se propagan a cualquier frecuencia 

B. La velocidad de fase es constante y coincide con la velocidad de grupo 
(si no hay perdidas) 
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Para este tipo de modos tendremos: 

E t = -V^e - 7 “ 2 


donde: 

A0 t = O 


junto con: 


4>ci = </>i; 4>C2 = <h\ ■■■ 



Figura 3.9: Analisis de los modos TEM 


Si integramos E t a lo largo de una li'nea cualquiera contenida en un piano 
2 = cte que una los dos conductores, obtendremos: 



-V t 0e-^.df= 


y oe -^ 2 


donde: 


V 0 — 4>2 — 4>i — 


que no clependera del camino elegido. 
Por lo tanto, la expresion: 



V t4> • dl 


V(z) = V 0 e~^ z 


define de forma univoca la diferencia de potencial entre conductores en 
cualquier seccion recta de la gufa (notar que depende de z). 

Por otra parte, la corriente total que recorre el conductor 2 en una seccion 
recta vale: 



h x H t ■ dl = 


C 2 


n ■ E t 
C2 V 


dl = zI 0 e~ loZ 


donde 


In = 




Jc 2 V 

Ademas, si aplicamos la Ley de Ampere al un contorno que envuelva al 
conductor exterior: 


H ■ dl = 


J c idl 


Cl 


dC‘2 d:l 


C2 



dado que el ultimo sumando es nulo, se verificara que la corriente en los dos 
conductores es la misma pero con sentidos opuestos. Entonces 


I(z) = / 0 e -7 ° z 
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define de manera univoca la corriente en cualquier seccion recta de la gufa. 

La relacion entre la tension y la corriente calculadas sobre una misma seccion 
recta vale: 


7 V 

Zo = 7 = v- 




dl 


—V(d> • hdl 


= 17 ■ cte. 


C2 


donde cte., factor adimensional, depende de la geometria del problema. Z a 
se denomina impedancia caracteristica. 


Los parametros 7 0 y Z Q se demominan parametros primarios de una linea 
de transmision ya que a una guia capaz de soportar un modo TEM se 
le denomina linea de transmision. Como hemos visto, podemos utilizar el 
formalismo de tensiones y corrienes que se usa en teorfa de circuitos. La 
representation grafica mas frecuente, independientemente del tipo de linea 
de transmision es: 


Observamos que las formas V = 
V( 0 )e —72 y I = l( 0 )e~ 72 correspon¬ 
dent a on das de tension y corriente que 
se propagan a lo largo de la estructura 


▲ 


o 


I 

-► 


V 


Zo 


o 


◄- 


I 


■c 


Figura 3.10: Linea de transmision 


La potencia media transmitida por el modo TEM viene dada por: 

W T = l -Re U E t x H* ■ ds = 

expresion equivalente a: 


1 VV* 

W T = -Re(VI*) = —Re 


= U ^Re{Z 0 ) 


La teorfa clasica de circuitos a partir de 
elementos concentrados deja de ser vali- 
da cuando el tamano de los elementos 
empieza a ser comparable a la longitud 
de onda 


que usaremos cuando utilicemos el modelo circuital y que depende de la 
impedancia caracteristica. 


3.3.1 _ 

Estudio de la guia coaxial 

El modo TEM se calcula mediante la solucion electrostatica de Ad) = 0 
(coordenadas cilindricas): 

13 / <9d>\ 13 2 $ 

p dp \ dp ) p 2 86 2 


junto con la condicion de contorno: 

d>(a, 4>) = V 0 ; 


®(b, </>) = 0 
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Figura 3.11: Lfnea coaxial 


cuya solution es: 


De aquf se sigue que: 


que junto con 


*(p,9) = V 0 


In b/p 
In b/a 


Z 0 = — In 
2tt 



/3 = uiy/pe r e 0 


caracterizan completamente a la lfnea sin perdidas. 

Se demuestra que el primer modo superior es el TEu cuya frecuencia de 
corte viene dada, de forma aproximada por: 


Sc 


c 2 
2-7T.y/fy a + b 


3.3.2 _ 

La lfnea como elemento circuital 


La lfnea de transmision describe rigurosamente el comportamiento del mo¬ 
do TEM independientemente de las dimensiones del sistema y la frecuencia. 
Pero la lfnea es tambien un elemento circuital. En este apartado caracte- 
rizaremos la lfnea de transmision como cuadripolo simetrico con perdidas 
(ver figura 3.121 


V, 



V, 


Figura 3.12: La lfnea de transmision como cuadripolo simetrico 
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Segun teoria de circuitos, la matriz ABCD tiene la forma: 

a b 


c a 


siendo a 2 — be = 1. Imponiendo 


V 2 

f Z o = 

1 2 


y teniendo en cuenta que: 


Vi 


—^ = v^a 2 + Vbc = e 7oi 
V 2 


se sigue que: 


6 = Z 0 sh"f 0 l 
sh"/ 0 l 

C= ^7 

a = ch^ 0 l 

cuyo circuito equivalente en T es: siendo: 

Z Z 


T 

I 


Figura 3.13: Equivalente en T de una Ifnea de transmision 




ch'yj, — 1 
sh"f Q l 


Y = 


sh'yol 


Sin embargo, dado que 

7o = 

tiene infinitas soluciones, no se puede establecer un circuito equivalente 
con elementos concentrados. No obstante, si particularizamos para el caso 
l « A, podemos desarrollar en serie las funciones hiperbolicas para llegar 
1 circuito de la figura [Xl4| y dado que: 


podemos escribir: 


Z 0 = K 





Y = lA = 
p z a 


= ju[iK 


1 

2 




Notar que la constante K introducida no 
es k = LOy/Jie 


por lo que para lineas de pequeno tamano si es factible utilizar un modelo 
de elementos concentrados como el representado en la figura |3~.15| 

Los parametros: 
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ZJU2 Zjt!2 , z.e 




s ' 




L 

V 

T 

X 


_ 





Figura 3.14: Circuito equivalents de una h'nea de transmision (l << A) 


U/2 U/2 



Figura 3.15: Circuito equivalents de una h'nea de transmision (l « A) 


■ L= inductancia por unidad de longitud 

■ G= conductancia por unidad de longitud 

■ C= capacidad por unidad de longitud 


dados por: 


T W 

L = jiK = p — 

V 

n e " // V 

g = “k = u1 ‘ Z 




se denominan parametros secundarios de la Knea de transmision y la rela¬ 
tion con los parametros primarios de la linea viene clada por: 


juL 


G + juC 



7o = VJ uj L(G + juiC) = Z S Y P 


siendo Z s e Y p la impedancia serie y la admitancia paralelo por unidad de 
longitud, respectivamente. En el caso de bajas perdidas estas expresiones 
resultan: 



7o = G^- + juspLC 

Por otra parte, los parametros secundarios pueden relacionarse entre si me¬ 
diant e: 


^=4 

G = u^-C 
e 
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donde C 


C = 


(fi c 2 —• hdl Q t 


SfVtt-dl 



es la capacidad estatica por unididad de longitud, parametro inedible o cal¬ 
culable resolviendo la integral anterior. Por lo tanto, vemos como la carac- 
terizacion de una linea de transmision se reduce a calcular la capacidad del 
problema electrostatico asociado. 

Como casos particulares interesantes, notar que: 

■ si Z a » predomina L 

■ si Z a « predomina C 

con lo que los circuitos equivalentes son: circunstancia que puede usarse al 


Li 


ce 


Zo » 


Zo « 


Figura 3.16: Circuitos equivalentes de una linea de transmision (l « A) (casos 
particulares) 

fabricar circuitos de microondas. 


3.3.3 


Linea de transmision con conductor no perfecto 

Las perdidas dan lugar a una componente E z en la superficie del conductor 
dada, aproximadamente por (efecto pelicular): 


Ez\c=^rfixH(0)\ c 

ad 


por lo que la solucion se parece mas a un modo TM (por tanto, en rigor ya 
no se trata de modos TEM). En tal caso, se habla de modo cuasi- TEM. 

En este caso introducimos en el modelo circuital una resistencia serie por 
unidad de longitud R, que esta relacionada con las perdidas en los conduc- 
tores. 


Ri/2 Li/2 


Li/2 Ri/2 



Figura 3.17: Linea de transmision con perdidas en dielectrico y conductores 
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Trabajando con las componentes transversales de campo, podemos definir 
los parametros primarios como: 

7 o = \/(-R + jtoL) (G + juiC) = \J Z s Y p 


R + jujL 
G + juC 



y en el caso de bajas perdidas en dielectrico con buen conductor, se puede 
realizar la aproximacion: 


7 =7 0 + a c 
Z 0 =r]K (1 — j 


\ ujy/uje 

esto es, para R « ojL y G « u>L, resulta: 

7 = jujyJLC ~ ju>VLC 
o lo que es lo mismo: 



R 

loL 



a 

P 




3.3.4 _ 

Cierre de la li'nea con una impedancia terminal. Impedancia transformada 


De momento no hemos especificado los valores de las tensiones V+, V a y de 
las corrientes /+, I~ que corresponden a las ondas de tension y corriente. 

Para determinarlas habra que anadir dos condiciones de contorno complejas, 
p. ej. los datos referentes a los dos extremos de la linea. 

Desde el punto de vista del empleo practico de una li'nea de transmision, en 
uno de sus extremos estara colocado el generador, con la mision de enviar 
potencia electrica (modulada por information o no) hacia el otro extremo, 
donde se captara en una impedancia terminal. Los valores de esta ultima 
y de la fuerza electromotriz y resistencia interna del primero, determinaran 
las const antes anteriores. 


Como se desprende de las deducciones posteriores, no estaremos, en general, 
interesados en obtener tanta precision y nos bastaran las consecuencias que 
implican la consideration exclusiva de la condition de contorno relativa a la 
impedancia terminal. Vamos pues a limitarnos a ella. 


Supongamos que en OO' (figura 3.181 termina la li'nea y que hemos conec- 
tado entre sus extremos OO' una impedancia terminal Zl. Si tomamos el 
origen de la coordenada z en OO ', las de los puntos de la li'nea seran negati¬ 
ves. Para evitar errores, vamos a introducir en lugar de la coordenada 2 otra 
que llamaremos /, igual a la 2 cambiada de signo, que sera siempre positiva 
y nos dara la distancia de un punto de la li'nea hasta su termination. 


Podemos decir que V+ da lugar a una onda incidente sobre dicha termina- 
cion y que V~ produce onda reflejada en la misma, con lo que se justifica 
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y, Zo 

z L 

◄- 

-► 


1 0 z 

Figura 3.18: Li'nea terminada 


introduzcamos una nueva notacion Ai y A r en vez de la anterior. Fijar la 
impedancia de carga es equivalente a imponer el cociente V(0)//(0) en ; = 0 

En el caso sin perdidas el voltaje total en la li'nea sera: 

V(z) = V o + e~ j0z + V~e +j0z 


y la corriente total 


I(z) 


y en z = 0 tendremos: 



YL 

z n 


e j0z 


, E(0) _ E Q + + V- 

L 1(0) V 0 + - Vo~ 


Introducimos un coeficiente de reflexion en tension T dado por: 


r _ Vq~ Zl - 

Vcfl Zl + Z a 

a partir del cual podemos escribir: 

V(z) = V+[e~ j0z +Te +j0z ) 

V{z) = Y^[e~i 0z -Te +j0z \ 

Notar que si Zi, = Z D sera T = 0 por lo que no existira onda reflejada. En 
este caso se dice que hay adaptacion de impedancias. 

La potencia media transmitida a lo largo de la li'nea sera: 

w T = \iv* = i(i - |r| 2 ) 

potencia que se entregara totalmente a la carga puesto que no hay perdidas. 
Observamos que: 


■ Si r = 0 => hay adaptacion, no existe potencia reflejada y toda la 
potencia incidente se transfiere a la carga. 

■ Si |r| = 1 => no tendremos potencia transmitida a la carga (tendre¬ 
mos emergi'a electromagnetica almacenada en la li'nea pero no habra 
energia viajera). Esto ocurre en dos casos: 

• r= -1 => cortocircuito 
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• r = 1 => circuito abierto 

Es frecuente trabajar con el valor conocido como perdidas de retorno dado 
por: 

RL = —20 log |T| 

En cuanto a la variation del modulo de la tension a lo largo de la linea 
tendremos: 


\V{z)\ = \V+\\l + Te 2 ^ z \ = \V+\\l + \T\e^ e ~ 20l) \ 


de forma que podemos definir una relation de onda estacionaria (ROE) de 
la misma forma que lo hicimos en las ondas planas: 


ROE = 


*min 


i + |r| 

i - l r l 


y tendremos dos casos extremos: 

■ ROE = 1: carga adaptada 

■ ROE = oo: carga completamente desadaptada 


El coeficiente de reflexion en cualquier punto de la linea sera: 

T(l) = r(0)e" 2 ^ = r(0)e 2j/3z = 


y la impedancia de entrada en cualquier punto de la linea: 


„ V(-l) V+e-M’ + V-e+V’ v 

in I(-l) V+e-iP* - Voe+V* '’ 

que operando, resulta: 

v , n _ v Zh cos 01 + jZ r) sin 01 

^in\ l) &O ry Ol , ■ ry ni 

Z Q cos pi + jZl sm pi 


La impedancia de entrada Zi n equivale 
al tramo de Imea, cerrado con Zl que 
hemos suprimido. Por lo tanto, podemos 
decir que es la impedancia equivalente a 
la parte eliminada o que es la impedancia 
transformada de Zl por la longitud de 
Imea l. 



(3.2) 


Figura 3.19: Impedancia de entrada Zi„ en cualquier punto de la linea 


Es de especial interes detallar los valores que toma Z, n en los casos parti- 
culares de Zl correspondientes a cortocircuito, circuito abierto y Z Q : 

■ Z L = 0 => Z in = Zo tan ((51) 

. Z L -»• oo ==> Z in -> Z 0 cot((3l) 
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Z r. — Z n 


7 _ 7 


Este ultimo resultado nos dice que si una linea de transmision se remata 
con su impedancia caracteristica, la impedancia transformada por cualquier 
longitud de linea es igual a dicha impedancia caracteristica. Observamos 
ademas, que las situaciones de impedancia 0 6 co se alternan con una pe- 
riodicidad de A/4 para los casos en que |T| = 1 


|V| 




|V| 

2|VJ 


AB1ERTO 


Figura 3.20: Lfnea terminada en corto o abierto 


Para completar este apartado, en el caso en que existan perdidas, la expre- 
sion (3.2) se convierte en: 


Z^chyl + Z 0 shjl 
Z a chyl + Z^shyl 


3.3.5 _ 

Voltaje, intensidad e impedancia relativas 


El cierre de una linea con Zl deja arbitraria P 0 + y tanto V(l ) como 1(1) 
resultan proporcionales a esta constante. Es natural que refiramos entonces 
V(l ) e 1(1) a V+, o mejor aun, a las ondas correspondientes a V/ + . De- 
fimimos los llamados voltaje e intensidad relativos en el lugar de la linea 
de coordenada l, V'(l) e I'(l) por las correspondientes ondas incidentes, es 
decir: 


no = 


v(i) 

Vfei 1 


= —e -27 ' + 1 
V+ 


IV) 


m 

Vo+e^/Zo 



e" 27i + 1 


De las definiciones se deduce inmediatamente: 

V’(l) 1 V(l) Z in (l) 

I’(l) Z a 1(1) Z 0 

resultando logico que introduzcamos la impedancia relativa en el punto l 7 
Z' in (l) cociente de Zi n (l) y Z 0: que valdra: 

7 , ,]S = Z L cosh (yl) + sinh(7 1) 
in{ 1 Z L sinh( 7 0 + Z 0 cosh( 7 l) { J 


Si tomamos la impedancia relativa correspondiente a Z l : 
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(3.61 se escribira: 


= Z' L cosh( 7 Z) + sinh( 7 l) 
in[> ^sinh( 7 0+cosh( 7 0 


Los casos particulares tratados al final del apartado anterior, se expresan 
entonces como: 


■ z 'l = 0 => Z L( 1 ) = tanh( 7 Z) 

■ Z' L —> oo => Z' in (l) —■> Z 0 coth( 7 Z) 

■ -^Z, = Z o —* Z in{l) = -Zo 


3.3.6 _ 

Lfnea equivalente a un modo 

Para un modo arbitrario la solution de los campos en sus componentes 
transversales viene dada por: 


E t = e{x,y){A+e- Jl3z + A~e^ z ) 


H t = h(x,y)(A + e~ jl3z - A~e jf3z ) 

donde e y h expresan la variation transversal del modo en cuestion. Ademas, 
se cumple que: 

r, s z x ?(x, y) 
h{x,y) = --- 

dado que E t y H t estan relacionadas por la impedancia de modo Z m . 

Si nos fijamos unicamente en la onda progresiva tendremos: 

E t = e{x,y)A + e~ J P z 
H t = h(x,y)A + e~^ z 

Buscamos ahora introducir una lfnea de transmision equivalente a un modo 
de manera que en una section dada por la coordenada z, se cumpla: 


P+= f s. d a= l -V{z)I{zY= l -VI * 


y para ello, expresamos los campos transversales como: 

e{x,y)_ v+e - j0z 


(3.4) 


E t = 


H t = 


Ci 

h(x,y) r~t ~ r il3z 

c 2 


Sustituyendo en la expresion (3.4 (: 


P+ = 


V+I+* 


1 


CiC: 


e x h* ■ zds 


1C 2 Js 


Dado que pretendemos calcular la potencia segun (3.4), tendra que ser: 

C-i C* = f exh* -zds 


(3.5) 
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Dado que tenemos dos constantes y una sola ligadura, debemos especificar 
una condicion mas para que queden completamente determinadas. Veamos 
que ocurre con la impedancia caracterfstica de la lfnea equivalente: 

_ V+ _ V~ _ A+C i _ Ci 
° ~ 7+ ~ 7^ ~ A-C 2 ~~ C 2 


Por lo tanto, si damos valor a esta ultima expresion, el problema estara ce- 
rrado. Las dos opciones mas comunes consisten en tomar: 


6 bien: 


7 _ Cl 

^ 0 — ^ — ^rr 

^2 


02 


Notar que la impedancia no queda determinada de forma absoluta, pues si 
multiplicamos las ondas de tension y corriente por a y 1 /a respectivamente 
(a £ 1Z por simplificar), esto no afectaria a la condicion (3.4( impuesta 
inicialmente, mientras que Z 0 se veria alterada segiin: 


Zo = 


aV+ 

l/al+ 


= a Z„ 


Esto no hace sino reafirmar que la Z Q no queda determinada y podemos 
fijarla segun nos convenga. Una vez que fijemos Z 0 , dado que la constante 
de propagation 7 es la que corresponde al modo, podemos calcular el circuito 
equivalente de la guia segiin: 


Z s = 7 Z 0 


Y = — 
P Z n 


impedancia serie y admitancia paralelo, respectivamente (ambas por unidad 
de longitud). 

Ejemplo: Modo TE\q 

Para este modo la cte. de propagation vale: 


7 = - ( fife ) 2 


Supongamos que elegimos: 

Z„ = Z m (TE„) = V1 _ ” {f/fc)2 

Se demuestra que las constantes Cl y C 2 valen: 



Por lo que se refiere al circuito equivalente, tendremos: 

Z s = jup 

Y p = jiue-je — 

LO 

y el circuito equivalente que resulta es (por unidad de longitud): 
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Propagante 



A1 corte 


Figura 3.21: Linea equivalente al modo TE\q 


3.3.7 _ 

Discontinuidades en gui'a 

De manera muy resumida podemos comentar: 

■ Siempre que tengamos una discontinuidad debida a un cambio de geo- 
metri'a se excitaran los infinites modos que la estructura es capaz de 
soportar. 

■ Cuando tengamos un cambio de dielectrico no aparecen modos dis- 
tintos del incidente. El problema de la onda reflejada puede tratarse 
mediante la linea de transmision equivalente. Esto es aplicable tam- 
bien a una terminacion en cortocircuito. 

■ Si sobre un obstaculo inciden dos modos que se propagan puede produ- 
cirse un intercambio energetico entre ellos (a pesar de las propiedades 
de ortogonalidad) debido a que cada modo es capaz de generar el otro. 
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3.3.8 _ 

Voltaje, intensidad e impedancia relativas 


El cierre de una li'nea con Zl deja arbitraria V+ y tanto V(l) como 1(1) 
resultan proporcionales a esta constante. Es natural que refiramos entonces 
V(l) e 1(1) a V a + , o mejor aun, a las ondas correspondientes a V+. De- 
fimimos los llamados voltaje e intensidad relativos en el lugar de la b'nea 
de coordenada l, V'(l) e I'(l) por las correspondientes ondas incidentes, es 
decir: 


V'(l) = 


vi!) 

Vo+e ^ 1 


Kr 

v 0 + 


o~ 2 1 l 


+ 1 


I'(l) 


m 

V+e^/Zo 


vr 

Vo+ 


e~ 2jl +1 


De las definiciones se deduce inmediatamente: 

V'(l) = = Z m (l) 

I'(l) Z a 1(1) Z a 


resultando logico que introduzcamos la impedancia relativa en el punto l, 
Z' in (l) cociente de Zi n (l) y Z 0l que valdra: 


Zl cosh(7^) + Z Q sinh(7Z) 
Zl sinh( 7 /) + Z a cosh( 7 /) 


Si tomamos la impedancia relativa correspondiente a Zl'- 

Z' T = ^ 


(3.61 se escribira: 


, = Z' L cosh( 7 1) + sinh( 7 Q 

in{ ) Z’ L sinh( 7 Z) + cosh( 7 0 


Los casos particulares tratados al final del apartado anterior, se expresan 
entonces como: 


■ Z'l = 0 =► Z' in (l) = tanh( 7 Z) 

. Z' L ->oo^ Z' in (l) -> coth( 7 /) 


Zl — Z 0 —» Z in (l) — Z, 
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3.4 EJERCICIOS 


3. 39 Sea una gufa rectangular rellena de aire en banda X cuyas dimen- 


siones son o=2.286 cm y b =1.016 cm: 


A. Calcular las frecuencias de corte de los cuatro primeros modos propa- 
gantes. 

B. ^Cual es la atenuacion en dB que corresponde a lm de gufa cuando 
trabajamos a una frecuencia /=10 GHz? (suponga que las paredes son 
de cobre) 

Sol : 


Modo 

m 

n 

fc (GHz) 

TE 

i 

0 

6.562 

TE 

2 

0 

13.123 

TE 

0 

1 

14.764 

TE, TM 

1 

1 

16.156 


A. a c = 0.11 dB/m; «d= 1.2 dB/m 


siendo a = a c + a,d 


3. 40 


Sea una gufa de seccion cuadrada de 


lado a: 


A. Calcular las frecuencias de corte de los modos 

B. Comprobar que existen dos modos dominantes 


Sol:: 

A r y/m 2 -\-n 2 

Jc ~ 2 a^/JCe 

B. TEi o y TE 0 i 

C. Wteio = ^2^° P 2 \ P ara H z = Pcosn/ax 


3. 41 | Se pretende transmitir un tono puro de frecuencia 11875 MHz por 
una gufa rectangular de conductor perfecto y dimensiones 10x6mm rellena 
de un determinado dielectrico. Calcular la permitividad relativa de este 
dielectrico para que nuestra frecuencia de trabajo este un 5 % por debajo 
de la frecuencia de corte del primer modo superior. 

Sol: 


A. 


= 4 


3. 42 


Sea una gufa rectangular de dimensiones a x b: 


A. Calcular los valores de x para los que las amplitudes de H x y H z del 
modo dominante son iguales. 

B. Calcular para que frecuencias dicho valor de x es a/4 


Sol: 
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A. x = ^ arctan ^ y cualquier valor cle y siendo (3 = — (n/a) 2 

B. / = V2fc 


3. 43 Considere un cable coaxial de tipo RG-142 con a = 0.035” y 6=0.116” 
y dielectrico con e r = 2.2. partir de que frecuencia aparece el primer mo- 
do superior TEn ? 

Sol: 


A. ,/c = 17 GHz 


3. 44 


Sea un cable coaxial de radios a y 6 (a < 6). Calcular: 


A. La impedancia caracteristica 

B. La potencia transmitida 

Sol: 

A. Z 0 = ^~ In 

y2 

para una onda incidente 


3. 45 


Calcular: 


A. Las frecuencias de corte de los clos primeros modos propagantes en 
una gufa circular con a=0.5 cm y e r =2.25 

B. Si la gufa es de plata y el dielectrico tiene perdidas (tan 6 = 0.001) 
calcular la atenuacion en dB correspondiente a un tramo de linea de 
50cm operando a /=13.0 GHz. 


Nota: Para el modo dominante: 


Oi c 


Re 

akrj/3 



(to 


SW: 

A. f c (TE n ) = 11.72 GHz-, f c (TM 0 1 ) = 15.31GHz; 

B. ad = 0.47 Nep/m; a c = 0.066 Nep/m; a = 2.38 dB 


3. 46 


Sea una gufa rectangular de plata (a = 6.1710 'S/m) de dimensio¬ 
ns a=0.86 mm y 6=0.43 mm. Calcular la potencia maxima que se puede 
transmitir trabajando clentro del ancho de banda monomodo si el maximo 
campo admisible es 30 KV/m 

Sol: 


A. P rnax = 2.08 KW 
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3. 47 | Una lfnea de transmision de Z 0 
dancia de valor Zl = 50 — 10 j O: 


500 se termina con una impe- 


A. Calcular el coeficiente de reflexion en el piano de la carga 

B. Calcular la ROE 

C. Calcular el coeficiente de reflexion en cualquier punto de la lfnea 

D. Calcular la impedancia vista en cualquier punto de la lfnea. 

E. Calcular la potencia que fluye por la lfnea. 

F. Calcular cuanto hay que desplazar la impedancia de carga hacia el 
generador para conseguir una situation de adaptation de impedancias 

Sol: 

A. |r(0) |=0.0995; Zr(0) = -84,2°; RL= 20 dB 

B. ROE=1.22 

c. r(/) = r(o)e-J' 80,ri 

ta 7 (1\ _ 7 Zl+jZq tan (31 

1J ’ ~ Z 0 +jZ L tan(3l 

e. p=gi(i_|r(o)| 2 ) 

F. No es posible 


3. 48 


Una impedancia de carga de valor Z^ = 130 + j’90 fi termina una 


lfnea de transmision de 50U de 0.3A. Calcular: 


A. El coeficiente de reflexion en la carga 

B. El coeficiente de reflexion a la entrada de la lfnea 

C. La impedancia de entrada 

D. La ROE de la lfnea 

E. Las perdidas de retorno 

Sol: 

A. r(0) = 0.598121.8° 

B. T(l) = 0.598|_165.8° 

C. Z in (l) = 12.7 + j5.820 

D. ROE=3.9 

E. RL= 4.46 dB 


3. 49 Una lfnea de transmision sin perdidas esta teminada por una carga 
de 1000. Si la ROE de la lfnea es 1.5, determinar los dos posibles valores 
de la impedancia caracterfstica de la lfnea. 

Sol: 
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A. Z 0 = 66.70 6 Z Q = 1500 


3. 50 


Una linea de transmision de impedancia caracteristica Z c = 500 
esta terminada en una carga de impedancia Zl = 25/0. Calcular a que dis- 
tancias l\ y I 2 el coeficiente de reflexion vale: 


A. 

r (h) = 1 


B. 

r (k) = -1 


Sol: 



A. 

h = 0.125A + 


B. 

l 2 = 0.375A + 

”5 


3. 51 


Una impedancia de carga de valor Zl = 80 
una linea con Z a = 1000 usando un tramo de linea 
l e impedancia caracteristica Z\. Calcular Z\ (real) 


+ j20 debe adaptarse a 
sin perdidas de tamano 

y 1. 


Sol : 


A. Zi = \/60O; l = ^anCeo A 
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CAPITULO 4 


Radiacion electromagnetica 


El problema basico consiste en determinar el campo electromagnetico pro- 4.1 INTRODUCCION 
ducido (radiado) por una estructura (antena) en cualquier punto del espacio 
partiendo de los parametros geometricos de la antena y de las condiciones 
radioelectricas de emision y reception (potencia y frecuencia). 

Basicamente, se utilizan dos procedimientos: 

■ Se calculan los potenciales electrodinamicos a partir de la distribucion 
de corrientes en la antena emisora (conocida o postulada). 

■ Se calcula el campo a partir de la distribucion de campo em conocida 
de antemano sobre una superficie cerrada S que encierra en su interior 
a todo el sistema radiante (ver apendice). 


Supongamos una distribucion conocida de cargas p y corrientes J en un 


volumen V' tal y como se indica en la figura 4.1 



4.2 CAMPO RADIADO EN 
FUNCION DE LAS 
CORRIENTES 


Figura 4.1: Distribucion de cargas y corrientes 

Deseamos calcular los campos E = E(r) y H = H(r). Para ello es conve- 
niente introducir unas funciones vectoriales nuevas que nos serviran para 
simplificar estos calculos. 

De la misma forma que hicimos en magnetostatica, dado que V • B = 0 
podemos introducir un nuevo vector A conocido como potential vector ( po- 
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tencial vector generalizado, potencial vector retardado o potencial de Lorenz) 
de forma que se cumpla: 

B = V x A 

Dado que A no queda determinado de forma unica, poseemos un grado de 
libertad, de manera que el A' definido como: 


A' = A + Vip 


tambien verificara la ecuacion. Este grado de libertad lo emplearemos cuan- 
do nos convenga. 

Sustituyendo en la ecuacion del rotacional de E se obtiene: 

Vx£ = —juB = —ju;(V x A) => V x (e + jtvAj = 0 
Por lo tanto, podremos escribir: 

E + jtoA = — V$ 


A1 valor $ se le conoce como potencial electrodinamico (6 potencial generali¬ 
zado 6 retardado). Notar que las cantidades A y $ son funciones del tiempo. 
Entonces, E y H pueden obtenerse facilmente mediante una diferenciacion 
por lo que los campos buscados se calcularfan como: 


E = — V$ — jcoA 
B = V x A 


(4.1) 


El siguiente paso es relacionar los potenciales electrodinamicos con las fuen- 
tes de los campos p y J. Operando sobre las ecuaciones de Maxwell, se llega 
a: 

/\A H - uj^pcA — — pjJ -\- V H - V • A^j (4-2) 

A4> + V • juA= (4.3) 

ecuaciones diferenciales locales que relacionan los potenciales electrodinami¬ 
cos con las fuentes de los campos, pero son enormemente complicadas. 
Ademas, en ellas aparecen los dos potenciales entremezclados. Para sim- 
plificarlas, podemos imponer una condicion de separabilidad aprovechando 
que A no esta univocamente definido: 


La ventaja de los potenciales electro¬ 
dinamicos es que sirven para obtener los 
campos mediante una simple derivacion 


jojpe® + V • A = 0 


Esta condicion se conoce como condicion de Lorentz. Aplicandola se llega a: 


Notar que en el caso estatico, las ecua¬ 
ciones se reducen a AA — pj, ecua¬ 
cion del potencial vector magnetostatico 
yA$ = —-, ecuacion de Poisson de la 
electrostatica 


A A + co 2 p,eA = —pj 


A 4> + w 2 /xe4> = - - 
e 


(4.4) 

(4.5) 


de manera que basta con solucionar la ecuacion correspondiente a A ya que 
se puede despejar de la condicion de Lorentz, calculandose entonces el 
campo electrico como: 


E = - juA - j 


. V(V - A) 


cope 


(4.6) 
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Por lo tanto, hemos llegado a un conjunto de ecuaciones que relacionan de 
forma implicita, dy$ con J y p (ecuaciones diferenciales). Buscamos ahora 
una relation explfcita entre estas cantidades, esto es, clos expresiones que 
nos proporcionen Ay $ en funcion de J y p (ecuaciones integrales). 


La ecuacion vectorial (4.41 se puede resolver formando tres ecuaciones es- 
calares. Para ello, descomponemos A en componentes rectangulares: 


A A = A A x x + AA v y + AA z z 


(4.7) 


Las coordenadas rectangulares se em- 
plean porque los vectores unitarios car- 
tesianos no son funcion de las coordena¬ 
das, por lo que el laplaciano puede ser 
factorizado 


donde k — co^fjJe. 

Las ecuaciones anteriores son identicas en su forma. Si resolvemos una de 
ellas las otras se pueden resolver de forma muy sencilla. En primer lugar 
buscamos la solution correspondiente al impulso unidad para, posteriormen- 
te, buscar la solucion correspondiente a una fuente arbitraria formada por 
una coleccion de infinitas fuentes puntuales. 

La ecuacion diferencial para una fuente puntual se convierte en: 

A'L + k 2 ^ = —6(x)S(y)S(z) (4-9) 

donde \f r es la respuesta a una excitation puntual en el origen de coordena¬ 
das. Notar que, a pesar de que la corriente tiene un caracter puntual, tiene 
caracter vectorial, por lo que lleva asociados una cierta direction y sentido. 
Si suponemos que la corriente esta dirigida segiin el eje z, esto es: 


Por consiguiente, tendremos que resolver: 

AA X + k 2 A x = —yJx 

A Ay k 2 A v — —(4.8) 
AA Z + k 2 A z = -yj z 


y = A z 


(4.10) 


para los puntos que no sean el origen de coordenadas, la ecuacion diferencial 
resultara: 

A 4' + k 2 ^ = 0 (4.11) 

Esta ecuacion se conoce como Ecuacion Escalar Homogenea de Helmholtz y 
corresponde a la ecuacion escalar de ondas. Dado que nos encontramos ante 
un problema de simetri'a esferica, 4' tendra unicamente dependencia radial 
y las soluciones tendran la forma ^—. Si nos quedamos unicamente con 
la solucion que se aleja del origen de coordenadas (la interpretation en el 
dominio del tiempo se realizara mas adelante) sera: 

p-jkr 

4/ = — (4.12) 

47T?’ 

Si la fuente se encuentra localizada en otro punto distinto al origen de 
coordenadas, la ecuacion anterior se transforma en: 


e —jkR 
47ri? 


(4.13) 


siendo R = |r — r |, | la distancia entre el punto donde se localiza la fuente, 


r', y el punto de observation P, r, tal y como se muestra en la figura 4.1 


Para una corriente arbitraria que tenga direction segun z, el potential vec¬ 
tor estara tambien dirigido segun z. Si consideramos que la fuente como 
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una coleccion de fuentes puntuales en un volumen V' ponderadas por la 
distribucion J 2 , la respuesta A z sera la suma de las respuestas puntuales 
dadas por (4.131: 


Az ~ 47T 


J, 


o-jkR 

R 


dv' 


(4.14) 


que junto con las ecuaciones analogas asociadas a las coordenadas x e y, 
dan la solution total: 



p -jkR 

J(r-)—— dv' 

V' *1 


(4.15) 


En el dominio del tiempo, estas expresiones tienen la forma: 




1 f p(r , ,t±R/c)_ dv , 


R 


A(r,t) = 7 - / J JAl±m dv ' 


47r 


R 


donde: 


c = 


\//^ 


tiene dimensiones de velocidad y corresponde a la velocidad de propagation 
de la perturbation electromagnetica en el medio en cuestion y R = \r — r*j 
corresponde a la distancia del elemento radiante al punto bajo estudio P. 

La solution que corresponde al signo menos, (t — R/c) se denomina solution 
retardada y la del signo mas, solution adelantada. Normalmente se elige 
la solution retardada ya que los potenciales en un instante t dependen del 
comportamiento de las cargas en instantes de tiempo previos a t (lo con- 
trario no tiene sentido fi'sico). Tenemos entonces, dos expresiones para los 
potenciales retardados: 


A(r, t) 


47T6 Jy, R 

p f J(f',t-R/c) , 

Tn Jy, - R - dV 


(4.16) 

(4.17) 


Estas ecuaciones resuelven el problema de encontrar los campos producidos 
por una distribucion arbitraria de carga y corriente. 

Observamos dos conclusiones importantes: 


A. la contribution de un elemento dv' al campo en un punto P depende 
de la excitation en el instante t — R/c. Por lo tanto, la perturbation 
se propaga con una velocidad c que depende del medio. 

B. Para un elemento de corriente o carga centrado en el origen de coorde¬ 
nadas, por simplificar, la contribution a la integral total tiene la forma 
J(t — r/c)/r 6 p(t — r/c)/r. Esto significa que el valor del modulo de 
los campos sera el mismo para todos los puntos situados a la misma 
distancia del origen r, o lo que es lo mismo, que existe simetria esferi- 
ca. Matematicamente, este tipo de expresion corresponde a una onda 
esferica , de tipo vectorial o escalar. 


Indicamos con un submdice A que los 
campos E y H se han obtenido a partir 
del potencial vector A 


lesumiendo, para calcular los campos generados por una distribucion de 
:orriente J se debe calcular A a partir de (4.141 y posteriormente aplicar: 
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H = 


-Vxl 


(4.18) 


y calcular E mediante (4.61, aunque suele ser mas simple emplear: 


E = — (V xH-J) 

jue 


en la region de las fuentes, y: 


E = 


—V x H 
jue 


(4.19) 


(4.20) 


en cualquier otro punto. 

Es habitual introducir un planteamiento paralelo utilizando corrientes magneti- 
cas (sin existencia fisica) como fuentes del campo radiado. Esto es asf ya que 
en algunas situaciones simplificadas pueden encontrarse corrientes magneti- 
cas ficticias que sirven para modelar un problema de una forma mas senci- 
11a. En el anexo A se indica el planteamiento completo del campo electro¬ 
magnetic© radiado por una antena con corrientes electricas y magneticas. 


4.2.1 _ 

Campo radiado por un dipolo infinitesimal 


Supongamos un dipolo infinitesimal de tamaho l muy pequeno comparado 
con A recorrido por una corriente I Q (fasor) y de grosor despreciable (ver 
figura 4.21. En este caso no existen corrientes magneticas por lo que el pro¬ 
blema se simplifica. Calculamos el campo em radiado segun el procedimiento 


z 



Figura 4.2: Dipolo infinitesimal 

expuesto. En primer lugar el potencial vector vendra dado por: 




J 


, p -jkR 

dv'=^~ 


V' 


R 


47T 


0 -jkR 

77/ ^ 

.. dl = 2-— 
R dirr 


.Lllo 


A/2 

1 - 1/2 


dz ' = Z^A e -3kr 
47rr 


Ahora calculamos el campo em a partir de las formulas vistas en el aparta- 
do anterior. Es conveniente realizar previamente un cambio a coordenadas 
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esfericas. Entonces: 

II J lp-jkr 

A r = A z cos 0 = - cos 9 

47t r 

II T lp-jkr 

A g = -A z sin9 = -sin 9 

4n r 

A,/, = 0 


y observamos que no existe dependencia con <f> por lo que tambien se veri- 
ficara que d/d(j) = 0 esto es, tendremos simetrfa de revolution con respecto 
a Se dice entonces que el campo obtenido es omnidireccional u omniaci- 
mutal. Entonces: 

H = 

[ir 


9 dA r 

d~r (TAl) - W 


de donde se sigue que: 


H r = H e = 0 

. kl 0 l sin 9 




4irr 


+ jkr 


0 -i kr 


y aplicando ahora (C.l): 


E r = rj- 


IJ cos 9 


Eg = jrj- 


E<b = 0 


2tt r 2 
kl n l sin 9 


1 + 


jkr 


0 —jkr 


47T r 


1 + 


jkr (hr) 2 


0 -jkr 


Notar que en la solution completa existen componentes longitudinales de 

E. 

Estudiemos ahora dos casos particulares en funcion del valor del producto 
kr. 


■ Si r << A (campo cercano) despreciamos los terminos en r/A: 

_ I 0 l cos 9e~j kr 

Er = ~ m 2^kP> 

. I 0 le~^ kr sin 9 
8 = 4tt kr 3 

Ecf, = 0 

H r = H e = 0 

I a le~ 2kr sin 9 
^ 4irr 2 

Observamos que los campos no incluyen el factor de propagation (la 
solution del campo magnetico coincide con el creado por una corriente 
rectilinea) 

■ Si r » A (campo lejano) predomina 1/r frente al resto: 

E r = 0 

„ . kl 0 l sin 9e~^ kr 

Ee = JV -:- 


E ( j > = 0 
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H r = 0 
H 0 = 0 

TT kl„l sin Oer jkr 
- 7=1 - 


y vemos que: 

A. El campo es transversal a la direccion de propagation 

B. E esta en fase con H 

C. E es perpendicular a H 

D. Localmente la onda esferica es asimilable a una onda plana. Ademas 
se verifica que: 

-* r x E 

H = - 

V 

donde rj es la impedancia intrinseca del medio. 


Por ultimo, la densidad de potencia radiada se calculara a partir del vector 
de Poynting: 


W = i Re{E x H*) = f^-\E\ 2 = f| 


kl 0 l 


47T 



Introducimos una serie de parametros que se emplean para caracterizar las 
antenas como elementos radiantes. 

4.3.1 _ 

Centro de fase 

Es el centro de las superficies equifase (esferas) a que pueden ser asimiladas 
a gran distancia los frentes de ondas de la radiation emitida por la antena. 

4.3.2 _ 

Intensidad de radiation 

Es la potencia radiada por una antena en una determinada direccion por 
unidad de angulo solido: 

U = r 2 W rad 

siendo: 

W = i. Re(E x H*) 


Para campo lejano se cumple: 


U^-[\E e \ 2 + \E^\ 2 } 


La potencia total radiada por la antena sera: 

n /* 27T n 7T 

Prad = j)UdSl = J I 


4.3 PARAMETROS 
CARACTERISTICOS DE UNA 
ANTENA 


La intensiadad de radiacion puede ser 
distinto para cada tipo de polarizacion, 
esto es se puede dar por separado para 
cada una de las componentes que cons- 
tituyen la base de polarizaciones 


0 


0 


U(&, <f>) sin 6d6d(j) 
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y para una fuente puntual isotropica (que radia igual en todas direcciones): 


Prad = (p UodtO = Uo 47T 


de donde se sigue que: 


U n = 


Prad 

4-7T 


Prad = (J) U(Kl = 


U sin Odddcj) 


0 Jo 


Se conoce como diagrama de radiation (en potencia) a la version normali- 
zada de la intensidad de radiacion: 

r(0, <t>) = U 


U„ 


4.3.3 _ 

Directividad 


Cuando la antena de referenda es 
isotropica, los decibelios empleados se 
denominan dBi 


Se define la funcion directividad como el cociente entre la intensidad de 
radiacion en una direccion con respecto a la intensidad de radiacion de una 
antena de referenda: 

D g =V 
9 U 0 

siendo U Q la intensidad de radiacion de una antena de diferencia que suele 
ser una fuente isotropica, si bien en ocasiones es comiin utilizar un dipolo 
6 una antena tipo bocina. Para una antena isotropica resulta: 

° 9 = 47 V P^d 

A1 valor maximo de esta funcion se le conoce como directividad: 


D = 


U„ 


U n 


= 47r- 


U n 


Pr 


ad 


Es frecuente expresar este valor como: 


D = 


47r 
^,4 


siendo: 


o2n 


4 = 


r(0, 4>) sin 9d6dcj) 


en funcion del diagrama de radiacion normalizado. 

Cuando el haz es suficientemente estrecho, el angulo solido Q 4 se puede 
expresar de forma aproximada: 

_ 4 7r _ 471- 25000 

0i( o )0 2 (°) 

donde 9\ y 62 son angulos pianos que corresponden a la anchura del haz 
principal a -3 dB. 


Notar que: 


■ La directividad da informacion de como se distribuye la energfa radia- 
cla por el espacio (en cada direccion) 

■ La directividad no aporta informacion sobre la potencia real que esta en- 
trando a la antena 
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4.3.4 _ 

Ganancia 

La directividad solo tiene en cuenta la distribution espacial de la energia 
radiada pero no aporta information sobre la energia que entra en la antena. 
Definimos la funcion ganancia como: 


siendo la potencia en bornas de la antena (concepto circuital). A1 valor 
maximo de esta funcion se le denomina ganancia. 


4.3.5 _ 

Eficiencia de radiation de una antena 


Podemos relacionar Pi n con P ra d mediante una eficiencia de radiacion de la 
antena: 


Prad 


C't Pm 


Esta eficiencia debe incluir perdidas ohmicas y desadaptaciones: 

ti C c CdC r — Ccd{\ |r| ) 


siendo: 


■ e c = perdidas debidas a los conductores 

■ ed= perdidas debidas a los dielectricos 

■ (1 — |T| 2 )= Perdidas debidas a desadaptaciones 

■ e cd= perdidas debidas a los dielectricos y los conductores (globalmen- 
te) 


donde T es el coeficiente de reflexion que se calcula segun: 

p _ A:n Z Q 

Zin + Z a 

siendo Z ln la impedancia equivalente a la antena (vista circuitalmente): 


Zin — Ra + JXa 


con 

Ra — Rr H~ Rl 

suma de: 

■ Ra : Resistencia de radiacion. 

■ Rl : Resitencia asociada a las perdidas 

La potencia que se disipa (de forma ficticia) en Ra es la potencia que se 
radia. 

Si consideramos que la antena esta recorrida por una corriente I , tendremos: 

P r = I 2 R, 

Pj = I 2 Rl 

P r R r 

Pr + Pj Rl + Rr 


por lo que: 
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Figura 4.3: Circuito equivalente a una antena 


4.3.6 _ 

Polarizacion 

Son aplicables directamente los conceptos vistos en ondas planas. 


4.3.7 _ 

Teorema de reciprocidad 


En lo visto anteriormente considerabamos la antena unicamente como trans- 
misor de energfa. Sin embargo, la antena tambien puede captar energia em. 
La relacion entre las propiedades de una antena receptora con las de una 
emisora vienen dadas por el teorema de reciprocidad. 


Supongamos una camara anecoica (sin eco) en la cual situamos dos antenas 
que pueden transmitir y recibir y que operan a la misma frecuencia (ver 
figura 4.41. Suponemos tambien que cada uno de los sistemas esta adpatado. 


A 



A’ 



Camara 

anecoica 


Figura 4.4: Teorema de reciprocidad 


En estas condiciones: 


■ Si A emite Pi, A’ recibe p\. 

■ Si A’ emite P 2 , A’ recibe p 2 . 


El sistema formado por ambos conjuntos equivale a un cuadripolo (determi- 
nado por los pianos de referenda S y S’). Como el cuadripolo es recfproco, 
el coeficiente de transmision de A a A’ es el mismo que de A’ a A: 


Pi 

Pi 


\A 2 = 


P2 

P2 
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Entonces, si las antenas estan ligadas a sistemas respectivamente adaptados 
en los que no existe mas que un tipo de ondas, la relacion entre la poten- 
cia recibida por el receptor y la emitida por el transmisor es constante e 
independiente de cual sea el sistema emisor o receptor. 

Por tanto, se pueden intercambiar en todo momento los papeles de emi¬ 
sor y receptor manteniendose sus caracterfsticas con independence de la 
operation que este realizando. Asi, p. ejlos diagramas de radiation en 
transmision y recepcion coinciden. 


4.3.8 _ 

Apertura efectiva (6 Area equivalente de absorcion) 


Una antena se utiliza en recepcion para capturar energfa electromagnetica 
de una onda. 

Se puede definir una apertura equivalente a la antena independientemente 
de la forma de esta: 


Potencia entregada a la carga Pt 
Densidad de potencia incidente Wi 


\It\ 2 Rt /2 

Wi 


La apertura efectiva no tiene por que coincidir con la apertura (superficie) 
fisica. 


Cuando la antena es plana, la relacion entre la apertura fisica (Ay) y la aper¬ 
tura efectiva (A e y) se conoce como eficiencia de apertura e ap , verificandose 
que: 

A e y = e ap Af 


Se demuestra que: 

0 < e ap < 1 

dependiendo e ap de como sea la distribucion de campo (en amplitud y fase) 
en la apertura que define la antena. 


4.3.9 _ 

Relacion entre directividad y apertura efectiva maxima 

Supongamos una situation como la de la figura |4.5| 


Tx Rx 



At, Dt " * Ar, Dr 


Figura 4.5: Directividad y apertura efectiva maxima 


Para la antena 1, a una distancia una distancia r tendrfamos: 


W t = W 0 D t = 


PtDt 

47rr 2 


La potencia entregada a la carga en el receptor sera: 


P r = W t A r 


P t D t A r 


El valor e ap = 1 se consigue para una 
distribucion de campo constante en am¬ 
plitud y fase a lo largo de la apertura 


47rr 2 
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entonces: 

D t A r = ^(47rr 2 ) 

It 

Si usamos 1 como transmisor y 2 como receptor y el medio es lineal, pasivo 
e isotropico: 

D r A t = ^(47rr 2 ) 

de donde se sigue que: 

Dt _ D r 
At A r 

y para los valores maximos sera: 

D t _ D r 
Atm A r m 


y esto es valido para cualquier par de antenas que consideremos. Se demues- 
tra que, para un dipolo infinitesimal se verifica: 


A 


em 


3A 2 

8n 


de donde se sigue que: 

\ 2 

A — _n 

n em — , o 

47T 

ecuacion que liga los parametros D 0 (caracterfstico de la transmision) y A em 
(de la recepcion). 


Si hay perdidas asociadas a la antena, la apertura efectiva maxima sera: 


\ 2 \ 2 \ 2 

Aem ~ Cf-j — Do = e c d( 1 — |r| 2 )-— Do = -— Go 

47T 47T 47T 


4.3.10 _ 

Ecuacion de Friis 

Supongamos dos antenas separadas una distancia r: 


■ Antena transmisora con potencia P y ganancia Gt 

■ Antena receptora con ganancia Gt 


A la distancia r, la densidad de potencia sera: 

PtG t 


W = 


47rr 2 


La potencia capturada por la antena receptora sera: 

P'-WA. = w£G'-™>£ 

47t (47t r) z 

conocida como Formula de Friis. Esta ecuacion puede expresarse como: 

2 


Pr 

Pt 


= e cdt ecdr{l- |r t | 2 )(l - |r r | 2 ) 


A 

AirR 


D t D r let e 


*\2 


PLF = 


e t ■ e 


*\2 


siendo: 
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el conocido como factor de perdidas de polarizacion que tiene en cuenta el 
desacoplo debido a la diferencia de polarizaciones entre las antenas trans- 
misora y receptora (se demuestra por reciprocidad) obtenida a partir de los 
vectores de polarizacion de las antenas emisora y receptora. 

En terminos de ganancias, esta ecuacion se escribe: 


Pr 

Pt 



2 

G t G r \e t • e* 


Cuando nos encontremos suficientemente lejos de la antena podemos cal- 4.4 APROXIMACION DE 
cular el campo total debido a una estructura radiante como la suma de CAMPO LEJANO PARA 
las contribuciones de elementos infinitesimales individuales (en virtud de la DIPOLOS FINITOS 
linealidad de las ecuaciones de Maxwell), realizando las siguientes aproxi- 
maciones: 

R « r ( modulo ) 

R « i — z' cos 6 (fase) 

de forma que sustituiremos (en las ecuaciones donde corresponda) 



Figura 4.6: Aproximaciones para campo lejano 


e —jkR 


R 


por: 


g—jk(r—z' cos 0) 


r 

que, de forma mas general puede escribirse como: 


p -jkr 

° „ikrri 

-gj 

r 


siendo r) el vector que apunta al elemento que crea el campo (Ver figura 


4.71: 


El procedimiento es entonces: 


Notar que el termino e^ kz COB 6 introdu¬ 
ce un cambio de fase relacionado con la 
diferencia de caminos. 


A. Dividir el problema 

B. Aplicar a cada elemento individual el factor corrector de fase 

C. Sumar (integrar) todas las contribuciones 
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A 



-» 


Figura 4.7: Aproximaciones para campo lejano (2) 


Si los elementos individuals son dipolos infinitesimales, tendremos (para el 
campo lejano): 



de forma que el campo total sera: 


r 1 / 1 Jrp-jkr rl/l 

Eg = dE e =jr ] —— sm6[ I{z')e jkz cos0 dz'} (4.21) 

J-U 2 47r ^ J-l/2 


y observamos que el campo tiene la forma: E=Campo elemento individual 
x [Factor funcion de geometria y excitacion] 

4.4.1 _ 

Aplicacion 

Aplicando todo lo anterior, se demuestra que para un dipolo corto (A/50 < 
l < A/10) 


i=4 4 IJ ;~ 3kr 

2 47rr 


cuya solucion de campo lejano es: 


Eg « jrj — -sin 9 


E r « Eq ~ H r « Hg ~ 0 



Eg/H ^ « r\ 


y la resistencia de radiacion: 



y para un dipolo finito recorrido por una corriente senoidal: 


I e = zsm[k(l/2 ± z')] 
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se obtiene: 



siendo la intensidad de radiation: 



\I 0 \ 2 fcos(fcZ/2cos#) — cos(fcZ/2) 


87t 2 sin 9 


Todo lo anterior es aplicable a los casos en los que tengamos multiples 4.5 ARRAYS 
fuentes iguales. A estas estructuras se las conoce como arrays. Los arrays 
pueden ser: 

■ Lineales 

■ Pianos 

■ Conformes 


Lineal 



Plano 



Conformado 


Figura 4.8: Tipos de arrays 


En todos los casos se trata de modificar las caracteristicas del campo radiado 
variando la geometria del array y la excitation de los elementos individuales 
Ii (fasores). El campo radiado sera: 


E t = + A 2 I 2 e jkp2f + ... + A N I N e jki>Nf 


siendo A, el campo radiado por cada elemento individual para una excita¬ 
tion unidad. 

Si todos los elementos radiantes son iguales, resultara: 


E t = A^he^ + I 2 e? krtf + ... + Iffe?**”*] 
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donde FA se conoce como factor de array y depende de la geometrfa y 


la excitacion (analogo al corchete de la ecuacion 4.211, de forma que el 
diagrama de radiation completo sera: 


= ri{9, (f>)\fa(9, <f)\ 2 


donde r, es el diagrama de radiacion del elemento y fa el factor de array 
normalizado. 

Atendiendo al tipo de excitation (iluminacion), un array puede ser: 


■ De iluminacion uniforme: los mismos A 

■ De fase uniforme: los mismos <pi 

■ De fase progresiva (rampa de fase): (ft = iAcf 


4.5.1 _ 

Arrays lineales uniforme equiespaciado 

Supongamos un array lineal de elementos isotropos equiespaciados. Fijamos 
la fase en el origen a 0° 



Figura 4.9: Array lineal equiespaciado 


El factor de array vendra dado por: 

N-l 

FA = I 0 + i ie i kdcose + l 2 e? 2kda » e + ... = £ 

n—0 

Supongamos una excitacion uniforme en amplitud y progresiva en fase: I n = 
A n e^ na . En tal caso tendremos: 

N-l 

FA=J2 A n e Mkdcose + a) (4.22) 

n— 0 

Introducimos una nueva variable 4': 

4' = kd cos 9 + a 

Entonces: 

FA=J2 Anj”* 

n—0 

funcion compleja que depende de 9. Observamos que: 

■ FA es una funcion periodica 27r en 4>. 

■ FA es una funcion de 9 pero no de </> (hay simetrfa de revolucion) 

■ Su estructura queda fijada en el intervalo 9 £ [0, 7r] 





4.5. Arrays 


129 


Caso particular: Array de amplitud uniforme y fase progresiva 


Supongamos que A 0 = A\ = A 2 = ... y que existe un desfasaje lineal 
progresivo a entre los elementos. En este caso, la suma de la serie (4.22 ) es: 


FA = Ao^^-^ 


2 


sin(lV4'/2) 
sin /2 


y prescindiendo del factor que afecta a la fase y normalizando se sigue: 


m 


1 sin(A^4'/2) 
N sin4'/2 


que se puede interpretar a partir de la figura siguiente, donde se representa 
un caso correspondiente a N = 4 elementos (con simetrfa de revolution 
en <f>): de forma que variando IV, d y a podemos disenar el diagrama de 



Figura 4.10: Margen de visibilidad 

radiation que queramos. 

A1 intervalo de variation de ip (valores barridos al movernos en la circunfe- 
rencia inferior) se le llama margen de visibilidad. 

Cuando el margen de visibilidad incluya mas de un lobulo del diagrama en 
4' diremos que aparecen grating lobes (multiples lobulos principales), tal y 
como se muestra en la figura |4.11 

Distinguimos dos casos: 

■ Array broadside: el maximo apunta a la direction 6 = 90° 

■ Array endfire: el maximo apunta a la direction 8 = 0° 6 9 = 180° (el 
haz es de tipo pincel) 
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Figura 4.11: Grating lobes 


Exploration del haz principal 

El maximo del factor de array ocurre cuando ’5 = 0. Si queremos que el haz 
apunte a una direction 0 o tendremos que aplicar: 

a = —kd cos 9 0 


4.5.2 _ 

Arrays pianos 

En este caso, el factor de array viene dado por: 

N 

FA = J2 -1 

n—0 


que se puede expresar como: 


FA = FA x FA y 

ya que el array piano es equivalente a un array lineal cuyos elementos indi- 
viduales son arrays lineales. 
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Calcular la directividad, la potencia radiada y la resistencia de ra- 
de un dipolo infinitesimal (l » A) 

Sol: 


| 4 . 53 

diacion 


4.6 EJERCICIOS 


A. U = A a sin 2 9 

B. Prad = A 0 8 tt/3 

C. R r = 80(7rZ/A) 2 ; D 0 = 3/2 


4 . 54 | Calcular la directividad de una antena cuya intensidad de radiation 

U = U n sin 0 


viene dada por: 


Sol : D = 4/tt = 1.27 


4 . 55 Encontrar la ganancia de una antena resonante cuya impedancia 
de entrada es 73f l cuando se alimenta con una ltnea de transmision de 
impedancia caracteristica 50 O suponiendo que el diagrama de radiation de 
la antena es: 

U ( 9 , 4>) = B a sin 3 9 


Sol: G = 2A4dB 


4 . 56 Estimar la directividad de una antena cuyos anchos de banda a -3dB 
en dos pianos perpendiculares son 9\ = 15° y 6*2 = 20°. Calcular la eficiencia 
de apertura si la antena tiene un area fisica de 0.1m 2 si f = 3GHz. 

Sol: D « 19 dB- e ap « 0.6 


4 . 57 


El campo electrico de una onda incidente viene dado por el fasor: 


Ei = xE 0 e- jkz 


y se recibe con una antena cuya polarization es lineal en esa direction segun: 

E a = (x + y)E(r, 9,4>) 

Calcular las perdidas por desacoplo de polarizacion (en unidades naturales 

y dB) 

Sol : l = 0.5 (L = 3 dB) 


4 . 58 Calcular la potencia recibida por un dipolo A/2 cuando se transmiten 
P=10 W a /=150 MHz con otro dipolo A/2 si R= lKm en los siguientes 
casos: 


A. Cuando ambos estan alineados 

B. Cuando el dipolo receptor esta girado 45° 

C. Cuando son perpendiculares 

Sol: 
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A. P r = 6.83 • 10 -7 W (-31.6 dBm) 

B. P r = 3.41 • 10" 7 W (-34.6 dBm) 

C. P r = 0 W 


4 . 59 


Calcular el diagrama de radiation correspondiente a los siguientes 


arrays: 


A. Dos fuentes isotropicas con amplitud identica y misma fase, separadas 
d = A/2 

B. Dos fuentes isotropicas con amplitud identica y fase opuesta, separa¬ 
das d = A/2 

C. Dos fuentes isotropicas con amplitud identica y desfase 90°, separadas 
d = A/4 

D. Dos fuentes isotropicas con amplitud identica y en fase separadas, 
d = A 


Sd: 

A. 

B. 

C. 


fi{6) = cos (| cos 0 s ) 
f 2 (d) = sin (|cos6>) 


f 3 (0) = cos 


^(cos0-l) 


D. 


fi(9) = cos(7rcos 9) 




APENDICE A 


Deduccion de las condiciones de 
contorno del campo EM 


Sean dos medios materiales, que llamaremos mediol y medio 2, de constan- 
tes di, ei y /ii y (i 2 , £2 y 2 respectivamente, separados por una superficie 
S (figura A.ll. Imaginativamente, tracemos las normales a S en todos los 


Medio 2 (ct 2 , Ll,) 



Figura A.l: Discontinuidad de medios materiales 

puntos de una linea cerrada situada sobre S y que delimite un area infini¬ 
tesimal AS. Estas normales definen una superficie reglada que en primera 
aproximacion sera un cilindro de generatrices perpendiculares a S. Conside- 
remos ahora el volumen infinitesimal interior a este cilindro y comprendido 
entre dos secciones rectas infinitamente proximas a S, una situada en el 
medio 1 y la otra en el 2. Supongamos, ademas, que la distancia entre estas 
secciones es tal que el area lateral del cilindro asf formado es un infinitesimo 
de orden superior al de AS. 

Si designamos con V 1 y V 2 los valores que toma V en las dos secciones rectas 
citadas, valores que corresponderan a puntos infinitamente proximos entre 
si y situados en los medios 2 y 1, respectivamente, el flujo total a traves del 
cilindro en cuestion valdra: 

{'D 2 -h-'D 1 - h)AS = {V 2 - Vi) ■ hAS 

si n es un vector unitario normal a Sen un punto cualquiera de AS. Hemos 
prescindido del flujo a traves de la superficie lateral por ser un infinitesimo 
de orden superior al escrito. 

Aplicando el teorema de Gauss (si es que puede aplicarse), este flujo sera igual 
a la integral de la divergencia de "D extendida al volumen del cilindro en 
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cuestion. Pero como (1.51 nos dice que, para cualquier volumen V: 

L 


V • Vdv = / pdv = Q 


esto es, igual a la carga contenida en el interior de V, llegamos a: 

(V 2 - 2>!) ■ fiAS = Q 


Si existe una distribucion volumetrica de carga, con una densidad cubica p 
finita, el segundo miembro de la igualdad anterior es un infinitesimo des- 
preciable frente al primer miembro, y entonces: 

(V 2 - Vi) • fiAS = 0 => (V 2 -V 1 )-h = 0 


Pero si sobre S' existe una distribucion superficial de carga, de densidad p s , 
(2? 2 - V\) ■ fiAS’ = p s AS' 


de donde se sigue que: 

n ■ {V -2 — Vi) = p s 


(A.l) 


alterando para mayor simetrfa del conjunto de formulas el orden de los 
factores del producto escalar. 


Procediendo identicamente con el vector B y haciendo uso de (1.51 se ob- 
tiene: 

(A.2) 


n ■ (i3 2 — B\) = 0 


Repitiendo el razonamiento para J empleando (1.31, se deduce que: 


(J2 


Jx) ■ nAS = 


d(AQ) 

dt 


donde A Q es la carga contenida en el interior del cilindro. 


Si la distribucion de carga es volumetrica, el segundo miembro es despre- 
ciable frente al primero, pero si es superficial, A Q = p s AS , y la igualdad 
anterior se expresa: 




(A.3) 



Figura A.2: Condiciones de contorno 


Para deducir las restantes condiciones de contorno, consideremos un rectangu- 
lo infinitesimal, situado en un piano ortogonal a S' y formado por dos lados 
de longitud A l, paralelos a S uno en el medio 1 y otro en el medio 2, y 
por dos lados normales a S y cuya longitud sea un infinitesimo de orden 
superior a Al. 











Tomemos un sentido de circulation de este rectangulo, que caracterizamos 
con el vector unitario r, dirigido segun uno de los lados de longitud A l, 
segun muestra la figura |Ad2) A este sentido le corresponde un sentido posi- 
tivo de la normal al rectangulo, que definimos con el vector n 0 . Las posibles 
orientaciones del rectangulo estan en correspondence biunivoca con los po¬ 
sibles vectores f, unitarios paralelos a S o con los vectores h 0 . Ademas los 
tres vectores f, fi 0 y fi forman un triedro trirrectangulo, siendo: 

f = h 0 x n 


La circulation de H a lo largo de dicho rectangulo valdra: 


H 2 ■ tAI -H 1 -tAI = t- (H 2 - Hi)Al 


siendo TL\ y Hi los valores de H en los lados del citado rectangulo paralo- 
los a S y situados en el medio 1 y en el medio 2 respectivamente. Hemos 
prescindido de la circulation a lo largo de los lados normales a S por ser un 
infinitesimo de orden superior. 


Si es aplicable a este area abierta el teorema de Stokes, dicha circulation 
sera igual al flujo del rotacional de H a traves del rectangulo, flujo que, de 
acuerdo con (1.1) consta de dos terminos: 


A. flujo de J= Intensidad de la corriente que pasa a traves del rectangulo 

B. flujo de dV/dt= Variation temporal del flujo de I? a traves del mismo 


El segundo siempre sera un infinitesimo de orden superior a A/, por ser 
proportional al area del rectangulo. En cuanto al primero, si las cargas 
electricas en movimiento dan lugar a una densidad superficial de corriente, 
J s , tambien sera un infenitesimo de orden superior a A l, pero si hemos 
admitido la posibilidad de una distribution superficial de carga sobre S', 
tenemos que aceptar la de su movimiento, que clara lugar a una corriente 
a lo largo de S. Para cartacterizar esta corriente sera necesario un vector 
densidad de corriente lineal J s tangente a S y tal que la intensidad que 
pasa a traves de un elemento de lfnea ds en el sentido del vector unitario a, 
normal a ds , valdra J s ■ ads. 

En consecuencia, el flujo de corriente producido por J s a traves del area 
del rectangulo y en el sentido de n 0 valdra J s ■ n 0 A/, que igualado a la 
circulation de H a lo largo del contorno conduce a: 

T ■ ( U-2 - Hi) = Js - ho 


Sustituyendo f por h a x n y alternado cfclicamente el orden de los tres 
vectores que aparecen en el producto mixto obtenido, llegamos a: 

{h 0 x n) ■ (H 2 - Hi) = [h x (H 2 - H{)\ ■ h 0 = J s ■ n 0 

Pero esta igualdad es valida para cualquier orientacion del rectangulo, es 
decir, para cualquier vector n 0 . Como tanto n 0 como n x (H 2 — Hi) y J s 
son vectores tangenciales a S, para que la igualdad anterior sea valida para 
cualquier n D , es condition necesaria y suficiente que: 


n x ( H 2 - Hi) = J s 


(A.4) 
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El mismo razonamiento, aplicado a £ y (1.21 lleva a: 


x =0 


(A.5) 


El sistema de ecuaciones (A.l l-( A.51 da las condiciones de contorno bus- 
cadas. Veamos que, como consecuencia de ellas, S es una superficie de dis- 
continuidad para los vectores del campo electromagnetico £, D, B y TL. 
Tomamos, p. ej. (A.lI y (A.5), particularizando la primera para p s = 0: 


h-(V 2 - Pi) = 0 
n x {£• 2 — £\) = 0 


Notar que: 

■ la primera expresa la continuidad de la componente normal de T> al 
atravesar S, lo que esta justificado puesto que h ■ V 2 es igual a la 
proyeccion de T>o sobre la direction del vector unitario n. 

■ la segunda afirma que la continuidad de la componente tangencial de 
£\ en efecto, si h es un vector unitario normal en un punto a una 
superficie S, (h x A) es la componente de A segun el piano tangente 
a S' en el mismo punto, cualquiera que sea A. En consecuencia, la 
componente tangencial de £\ coincide en valor, direction y sentido 
con la de £ 2 , lo que implica que £\, £2 y la normal a S estan en un 
mismo piano. 

Llamando £ 2 ™ y £ 21 , £1 n y £ 11 , a las componentes normales y tangenciales de 
£2 y £ 1 , respectivamente, y teniendo en cuenta que V = eV , las condiciones 
anteriores pueden escribirse: 


£ 2^2 n — Cl£l n 
£ 2 1 = £it 

Los angulos 02 y 0i que £ 2 y £\ forman con la normal seran tales que (figura 

All): 


, a £21 
tan 02 = -pr- 
£2 n 

, a £ u 

tan 0i = —— 
C'ln 

por lo que se cumplira: 

tan 0 2 _ £2 
tan 0i ei 



Figura A.3: Componentes del campo electrico 













Por lo tanto, salvo si e 2 = ei, sera 61 ^ 0 2 y las Kneas del campo electrico 
se refractan al pasar de un medio a otro. Como, ademas: 


l ^| 2 = |£fj + 1 ^ 2*1 = ~i £ ln + \£lt\ 
e 2 


Es evidente que £2 7 ^ |£i|. La superficie S es una superficie de disconti- 
nuidad para £, discontinuidad que afecta tanto al valor de las componentes 
de £ como al de las derivadas primeras. No ofrece dificultad repetir el ra- 
zonamiento para los restantes vectores , V, H y £>, llegandose al mismo 
resultado. 


Pero si S es una superficie de discontinuidad no se puede aplicar el teorema 
de Gauss al cilindro empleado para obtener (A. 1), (A.2) y (A.31 ni el de Sto¬ 
kes al rectangulo introducido para (A.41 y (A.5). Los resultados cleducidos 
invalidan el metodo de deduccion. 


Desde el punto de vista del rigor, por lo tanto, es preferible postular las 
condiciones de contorno (A.l l-( A.5), como lo era tambien postular el sistema 

([lT|i-([l 8 |. 
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APENDICE B 


Promedio de valores complejos 


Recordemos que si A es una magnitud escalar o vectorial compleja, que B.l VALOR MEDIO DE UN 

varfa periodicamente con el tiempo, su valor medio < A > en un periodo r PRODUCTO 

vale 

< A >= - [ Adt 
T Jo 


y que 


Re(A) 


A + A* 
2 


simbolizando con A* el complejo conjugado de A. 

Sean ahora A\ y A 2 dos campos complejos variables en el tiempo segun la 
ley 

27T 

A 1 (x,y,z,t) = A / 1 (x,y,z)e JU]t -, u>= — 

T 


y pasamos a calcular el valor medio del producto 


Re(Ai) © Re(A 2 ) 


en que simbolizamos con 0 el producto en el campo complejo si Ai y A 2 son 
magnitudes escalares y los productos escalar y vectorial si son vectoriales. 

Puesto que, en virtucl de la propiedad distributiva de tales productos con 
relacion a la suma 


Re{A x ) © Re(A 2 ) 


A 1 + A\ ^ A 2 + AI 
2 ® 2 


1 

4 


(Ai © A 2 + Ai © A* 2 + A{ © A 2 + Al © A* 2 ) 


sera 


< Re{A x ) © Re(A 2 ) >=< ~(A t & A 2 + Ai © A^ + A{ © A 2 + A\ © A* 2 ) >= 


1 

4 


(< Ax © A 2 > + < A x © Al > + < A{ © A 2 > + < A\ © Al >) 


ya que el valor medio de una suma es la suma de los valores medios de los 
sumandos. 
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Pero, 


< A-! © A 2 > =< A[ © A 2 e 2juit >= 


-A'iQA ' 2 


2 n 


to 

2 n , 

2 jut 1 27r /^ 


r2n/cj 


A\ © A'e 2juJt dt = 


= 0 


L 2 J i 

<AlQA* 2 >=< 4* © A 2 *e~ 2jut >= 0 

u r 2v/u 

< A 1 © A 2 > =< Al © A 2 >= — / A[ © A 2 dt = 

27t J o 

—Ai © A 2 * / dt = A\ © ,4 2 * = Ai © ^2 

^ </0 


En consecuencia 

Re(Ai) © i?e(A 2 ) 
en virtucl de ser 


-Re(Ai) © Re(A 2 ) 


1 

4 


(.41 QA* 2 + AlQ A 2 ) 


1 Ai © A% + (^li © A* 2 )* 

2 2 


',^(,1, ©.,1$) (B.l) 

^ (Ai © j4 2 + A* © A 2 ) 



APENDICE C 


Corrientes magneticas 
equivalentes 


Aunque las corrientes magneticas no existen, es conveniente introducirlas 
ya que, matematicamente son utiles para tratar el problema (aparecen al 
hacer uso de ciertos teoremas cle equivalencia, como se explica en el anexo 
siguiente). 

Por ello, introducimos un nuevo vector M= Densidad de corriente magneti- 
ca. De esta manera simetrizamos las ecuaciones de Maxwell: 

Vx£ = — M — jco^H 

Entonces, podemos definir un nuevo potential vector F de manera que: 

Ep = — -V x F 
e 

que por un procedimiento analogo al anterior conduce a: 

Hp = - jvF 

donde <p m es el potential electrodinamico magnetico, verificandose: 

A F + k 2 F = —eM 


Si imponemos: 


resulta finalmente: 



Hp = -juF - — V(V • F) 

ujfie 

por lo que la solution completa del problema viene dada por (indicando el 
subfndice A 6 A si la fuente de campo es una corriente electrica o magnetica): 


E = Ea + Ep 
H = H A +H F 


—juA - 

Vo. 

<1 

< 

A) — -V x F 

(C.l) 


cj/ie 

€ 

-juF - 

- —V(V 
uj/ie 

■ F) — -V x A 

(C.2) 




142 CORRIENTES MAGNETICAS EQUIVALENTES 


que se pueden expresar tambien como: 

E = Ea + E f = —— V x Ha - -V x F 
jue e 

H = Ha + Hp = - V x Ep H— V x A 

jujfi n 


como funcion de Ha y Ep, siendo A y F las soluciones de: 



e —jkR 

R 



e ~jkR 

R 


dv' 

dv' 


(C.3) 

(C.4) 



APENDICE D 


Principio de equivalencia 


Veamos otra forma de calcular el campo radiado por una apertura con una D.l CAMPO A PARTIR DE 

distribution de campo conocida. LA DISTRIBUCION EN UNA 

Segiin el principio de equivalencia, es posible reemplazar una antena de aper- APERTURA 

tura con corrientes equivalentes que producen campos radiados identicos a 

los de la antena original. El principio de equivalencia se deriva del hecho de 

que una solution de las ecuaciones de Maxwell junto con unas determinadas 

condiciones de contorno, es la linica solution posible. 

En el problema original, los campos que satisfacen las ecuaciones de Maxwell 
en la region exterior al volumen V y los que satisfacen las condiciones de 
contorno en la superficie S son unicos. Mientras las fuentes exteriores a V 
y las condiciones de contorno en S, superficie que engloba a V, no cambien, 
no cambiara la solution (E, H). 

En el problema equivalente, las fuentes exteriores a V no cambian, dado que 
no existen. En cuanto a las condiciones de contorno, se trata de introducir 
los cambios adecuados en el problema inicial para que tampoco cambien. 

Supongamos que en el problema original, los campos en la superficie S son 
E(S) y H(S). En el problema equivalente eliminamos las fuentes originates 
(dentro del volumen V) e introducimos en la superficie S las corrientes apro- 
piadas para que se sigan cumpliendo las condiciones de contorno originates. 

Si los campos internos a S son E\ y Hi, y los externos E y H, las corrientes 
equivalentes se pueden obtener a partir de las condiciones de contorno sobre 
S son: 


j Sl =hx{H- Hi] 

(D.l) 

Ms i = [E- Ei] x h 

(D.2) 


Estas corrientes equivalentes, que se obtienen unicamente a partir de las 
componentes tangenciales de los campos, se pueden emplear para calcular 
los campos fuera de la superficie S. Sin embargo, es necesario conocer los 
campos sobre S para poder determinar dichas corrientes equivalentes. Tam¬ 
poco sabemos todavia como calcular los campos externos a partir de las 
corrientes equivalentes. 

Dado que los campos E\ y H\ son arbitrarios, podemos hacerlos cero por 
simplicidad. Entonces las ecuaciones anteriores se convierten en: 

J s = fix H(S) 


(D.3) 
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M s = E(S ) x h (D.4) 

donde E(S) y H(S) son los campos sobre la superficie S. Esta formulacion 
se conoce como Principio de Equivalencia de Love. 

Dado que los campos en el interior de S son cero, en el problema equivalen- 
te podemos introducir materiales dentro de S. Si colocamos un conductor 
electrico perfecto a lo largo de S, Js desaparecera y unicamente tendremos 
la densidad de corriente magnetica Ms en presencia de un piano conductor 
perfecto generando los campos. De manera analoga, si situamos un conduc¬ 
tor magnetico perfecto, Ms se anulara. 

En cualquier caso, los problemas mencionados son dificiles de resolver siem- 
pre que S sea una superficie arbitraria. Notar que si la antena contiene zonas 
conductoras, Js debe ser igual a la densidad de corriente sobre tal portion, 
que contendra Js y Ms antes de que se introduzca ningun conductor ficticio. 

Si la superficie S es grande en terminos de A y la curvatura de S es pequeiia, 
se puede aplicar localmente la teorfa de las imagenes sobre la superficie 
para encontrar las corrientes que operan en presencia de los conductores 
introducidos. Sin embargo, dado que la superficie S la hemos introducido 
por conveniencia y estamos interesados en problemas de radiation, podemos 
extenderla hasta oo y como S debe ser una superficie cerrada, suponemos 
que S se encuentra en z = 0 y todas las fuentes en z < 0. En este caso, 
podemos aplicar teorfa de imagenes de forma exacta y resolver facilmente 
muchos problemas practicos relacionados con las aperturas planares. 

A partir de las corrientes Js y Ms en z = 0 aplicamos la teorfa de imagenes 
para reducir la formulacion al semiespacio z > 0. Tendremos tres posibili- 
dades: 

■ Espacio libre: Js y Ms 

■ Conductor electrico perfecto: Js = 0 y 2 Ms 

■ Conductor magnetico perfecto: 2 Js y 2 Ms = 0 

Los campos ( E , H) en la region z > 0 se calcularan evaluando previamente 
A y el resultado sera exacto en z > 0. No obstante, si restringimos nuestra 
solution a la zona lejana, el procedimiento se simplifica considerablemente: 

A. Obtener A 

B. Calcular Ea 

C. Obtener F 

D. Calcular Hp 

E. Calcular Ep = TjHp x r 

F. Sumar los dos resultados: E = Ea + Ep 

En el caso en que trabajemos con conductores perfectos, los calculos se 
reducen considerablemente. Resumiendo, las posibilidades son: 


Js y Ms- E = Ea + Ep 
2 Js- E = 2 Ea 
2 Ms'- E — 2Ep 



D.l. Campo a partir de la distribution en una apertura 


Una vez que ha quedado aclarado el procedimiento, falta por ver como 
determinar las corrientes en la apertura. Hasta ahora, no se ha introducido 
ninguna aproximacion en el desarrollo, salfo las relativas al campo lejano. 
Sin embargo, es muy dificil disponer de la distribution de campos E(S) 
y H(S). Lo normal es disponer de information aproximada sobre dichos 
campos en una portion finita del piano de la apertura. Tal aproximacion se 
conoce como Aproximacion de Optica Fisica, en la cual, se asume que los 
campos en la apertura E a y H a son los de la onda incidente y se supone que 
estpos solo existen en una portion finita de la apertura S a siendo cero en 
el resto de S. En la mayoria de los casos, S„ coincide con la apertura fisica 
de la antena. Esta aproximacion mejora a medida que las dimensiones de la 
antena con relation a A aumentan. 

El procedimiento completo se resume a continuation. Supongamos que cono- 
cemos los campos E a y H a que exiten en S a y son tangentes a dicha superfi- 
cie, portion de la superficie infinita S. Estos campo han podido determinar- 
se, p. ej. mediante la aproximacion de optica fisica. Las corrientes lineales 


en las superficie equivalentes son: 

J s = hxH a (D.5) 

Ms = E a x h (D.6) 

sobre S a y cero en el resto. Entonces los potenciales vectores son: 

A = p < —Lfix [[ H a e jkf "ds' (D.7) 

47rr JJs a 

P = ff E a e^ds' (D.8) 

y dado que el campo lejano responde a: 

E = —juA — jwrjF x f (D.9) 


tras ciertas manipulaciones se sigue que: 


E = —jkf x 



rjr x (h x H a )\e ■ jfer ' r ds' 


(D.10) 


analoga a la formula de Goudet. 

Las integrates que aparecen en estas expresiones de Ay F son transformadas 
de Fourier bidimensionales. Si definimos: 


P = JJ E a e^ kf ' r ds' 

Q = /I Hands' 

podemos escribir: 

P x = U E ax {x\y')e jk{x ' sinScos,p+v ' sindsin ^dx'dy' 

P y = If E av {x', y y^ x ' s ^ ecos<l,+v,sinesinrt,) dx'dy' 
Q x = ff H ax (x\y')e ik ( x ' s ™ 9ms * +y ' sin0s ™®dx'dy' 


(D.ll) 

(D.12) 

(D.13) 

(D.14) 

(D.15) 
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Qy = JJ E ay {x',y')e jk ^' sinecosMs ^ es ^^dx'dy' (D.16) 

que junto con h = z conducen a: 

_ g -jkr 

A = ii— - (~Q y x + Q x y) (D.17) 

47t r 


Para aperturas montadas sobre un piano 
conductor se prefiere utilizar un piano 
conductor electrico infinito, ya que asf se 
cumple la condicion de campo electrico 
tangencial nulo en el conductor 


F=-e e -^{-P y x + P x y) (D.18) 

47rr 

Notar que si se conocen exactamente los campos en la apertura, las tres for- 
mulaciones proporcionan el mismo resultado. En la practica, dado que solo 
se dispone de information aproximada sobre los campos en la apertura, las 
tres formulaciones proporcionan resultados distintos, si bien, las diferencias 
no suelen ser significativas. 



